general purpose boundary point method for the optimization of the 2DM. No symmetries implemented