geometry-max-connections
I would like this feature to be set on auto. Currently, I am using Long-range geometry that is using connections that extend up-to 12 sites. For example: Connection SiSj(i,j=i+12) != 0 but SiSj(i,j!=i+12) = 0. Currently, If use geometrymaxconnection=12, DMRG++ will save all operators from site i to i+12. However, I only need to save operators of site i and site i+12 at the center-of-orthogonality i in DMRG++ step. At this center, DMRG++ redundantly save operators that should not be stored taking 11 to 12 times more RAM per step than required.
In fact, all stored operators have to be transformed, so DMRG++ also takes more time than necessary.
It can actually be made even more efficient by being selective which operator. For example in Kitaev model. I can have SxSx(i,i+12) !=0 SySy(i,i+2) != 0 SzSz(i,i+1) != 0
Then, at the center-of-orthogonality at i, I need to save only Sx[i], Sx[i+12], Sy[i], Sy[i+2], Sz[i], Sz[i+1].
This can significantly reduce the memory requirement of the code, especially when operators become dense in highly entangled phases like spin-liquids.