I wrote some open source code a few years ago that introduces a simple way to convert LU factorizations into a Determinantal Point Process sampler and used large Aztec diamonds as the example

https://gitlab.com/hodge_star/catamari