Implementation of a primal-dual semidefinite optimization program for the variational determination of the 2DM. Block digaonalization is implemented for Hamiltonians which have SU(2) symmetry.