I have added some operators like DES and CRE_DES_DESCOMP which are required when the bra nd ket states are different
I am adding new operators so we can calculate Hamiltonian matrix elements between two arbitrary mps states