julia> orthogonalize(ψ.AL,:).C
InfiniteMPS
[1] ((dim=20|id=746|"Left,Link,c=1,l=1") <Out>
1: QN("SzParity",0,2) => 10
2: QN("SzParity",1,2) => 10, (dim=20|id=746|"Link,Right,c=1,l=1") <In>
1: QN("SzParity",0,2) => 10
2: QN("SzParity",1,2) => 10)
[2] ((dim=20|id=248|"Left,Link,c=1,l=2") <Out>
1: QN("SzParity",0,2) => 10
2: QN("SzParity",1,2) => 10, (dim=20|id=248|"Link,Right,c=1,l=2") <In>
1: QN("SzParity",0,2) => 10
2: QN("SzParity",1,2) => 10)
The resulting
Cmatrix fromorthogonalize(psi, :)has links that have the same id with different tags. We should instead have different ids so we can contract the tensors properly