variational determination of the reduced density matrix for the 1D Hubbard model using three different semidefinite programming algorithms