Hi there,

I am using the cisdvec_to_amplitudes function in ucisd.py. I found that when transforming the myci.ci to the amplitudes, the 4-fold degeneracy is not considered for the double excitations within the same spins. Therefore, if we directly take the value returned: c0, c1, c2 = myci.cisdvec_to_amplitudes(myci.ci) and construct the CISD wavefunction, the wavefunction is incorrect. Therefore, I need to manually divide c2[0] and c2[2] by 4.

I wonder if this is the convention or something to be fixed.

Thank you for your time!