I've modified the code to work for a square (fully) OBC cluster.
It seems to test well on a 4x4 square with a 2x2 region A, which
according to
http://rgmelko.webfactional.com/index.php?title=TFIM_EEs_from_Lanczos
should give a S_2 = 0.08088534039