I'm trying to calculate observables from idmrg. The scheme used in the sample "idmrg.cc" is base on the right-orthogonal gauge. But when I use the following code to check the orthogonality, it turns that this condition is not satisfied in the middle of the unicell
int Nuc = 4;
int N = 2*Nuc;
// Other codes
auto indc = commonIndex(pa,psi.A(ni-1));
where right-orthogonality is satisfied except for the center tensor, i.e., ni=5. The result for that site is not an identity matrix.