A Riemannian Approach To Large-Scale Constrained Least-Squares With Symmetries