Fua: Constrained least squares