Printing Symmetric Hamiltonians with identity

[VicVanderLinden]

If you create a symmetric Hamiltonians containing the identity

lat = FiniteChain(N);
idd = id(U1Space( 1 => 1,0 => 1,-1 => 1)) 
H_heisenberg = @mpoham sum(J1*((S_exchange(ComplexF64,U1Irrep,spin=1))⊗idd){lat[i],lat[i+1],lat[i+2]} for i in 1:N-2)

or

idd = id(ComplexF64,SU2Space(1 => 1)) 
H_heisenberg = @mpoham sum(J1*((S_exchange(ComplexF64,SU2Irrep,spin=1))⊗idd){lat[i],lat[i+1],lat[i+2]} for i in 1:N-2)

you encounter issues when printing the Hamiltonian (It still creates it perfectly).

TypeError: in typeassert, expected DataType, got Type{Union{BraidingTensor{ComplexF64, GradedSpace{SU2Irrep, TensorKit.SortedVectorDict{SU2Irrep, Int64}}}, TensorMap{ComplexF64, GradedSpace{SU2Irrep, TensorKit.SortedVectorDict{SU2Irrep, Int64}}, 2, 2, Vector{ComplexF64}}}}

[lkdvos]

Thanks for reporting this!
I think this was an issue in TensorKit, that has been solved by QuantumKitHub/TensorKit.jl@db09a2f, but unfortunately has not been released yet.

At least for me, I can locally resolve this by using pkg> add TensorKit#main, although I'm not sure if this breaks other things. In any case, this will be resolved as soon as a new TensorKit version is released and we update MPSKit to be compatible with that.

I'll leave this open as a reminder for that, thanks again for reporting!