Orthogonal projection on a subspace