Carsten wrote a checkerboard decomposition before I started working with this code. Since then BlockDiagonal and Complex StructArrays have been introduced, but checkerboard decompositions haven't been touched. It would be good to move the implementation to a matrix wrapper type and generalize it, if necessary, to arbitrary lattices.
Carsten wrote a checkerboard decomposition before I started working with this code. Since then BlockDiagonal and Complex StructArrays have been introduced, but checkerboard decompositions haven't been touched. It would be good to move the implementation to a matrix wrapper type and generalize it, if necessary, to arbitrary lattices.
This might be a good reference: (Appendix A1) https://opus.bibliothek.uni-wuerzburg.de/opus4-wuerzburg/frontdoor/deliver/index/docId/19122/file/Beyl_Stefan_Thesis.pdf