Closed
Description
Would be nice to implement in the near future the computation of the 4RDM for DOCI wave-functions. This computation can be added to the function compute_rdms_34 in the source code.
The non-zero elements of the DOCI-4RDM are pretty straightforward to compute, as are the non-zero spin blocks of the 4RDM. The major complication in its implementation is figure out the best way to include all the symmetries the 4RDM has. I think I can do my best to implement this as soon as possible.