Projection

If $$A$$ is skinny and full rank, then the projection operation is $$\mathcal P_{\mathcal R(A)}(y) = A(A^TA)^{-1}A^T y$$

TBC