Implementation of three different semidefinite programming algorithms for the optimization of the 2DM, tailored to the 2D Hubbard model. All spatial symmetries have been exploited.