Implementation of a dual-only potential reduction semidefinite optimization program for the variational determination of the 2DM, in which the Hessian is explicitly constructed for the solution of the Newton equations.