Bug: different reduced density matrices calculated before and after center_orthogonalization

[ifsheldon]

length = 8
physical_dim = 3
virtual_dim = 4
dtype = torch.complex128
device = torch.device("cpu")
center = 3

para = {"length": length, "d": physical_dim, "chi": virtual_dim, "dtype": dtype}

psi = MPS_basic(para=para)
print(f"psi.center: {psi.center}")
reduced_density_matrix_psi_before = psi.one_body_RDM(center)
psi.center_orthogonalization(center, "svd")
print(f"psi.center: {psi.center}")
reduced_density_matrix_psi_after = psi.one_body_RDM(center)

assert torch.allclose(reduced_density_matrix_psi_before, reduced_density_matrix_psi_after), f"reduced_density_matrix_psi_before: {reduced_density_matrix_psi_before}\n\nreduced_density_matrix_psi_after: {reduced_density_matrix_psi_after}"

This assertion will fail.