# The tensors calculated by idmrg do not obey the right orthogonality condition.

Hi Miles,

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
for(auto ni=1;ni<=N;++ni)
{
auto pa=psi.A(ni);
auto indc = commonIndex(pa,psi.A(ni-1));
PrintDat(pa*prime(dag(pa),indc));
}


where right-orthogonality is satisfied except for the center tensor, i.e., ni=5. The result for that site is not an identity matrix.

Best,
Su Wei

psi.A(ni-1)