Description: Define orthogonal projection onto a subspace. This is just a wrapping of df-pj1 , but we restrict the domain of this function to only total projection functions. (Contributed by Mario Carneiro, 16-Oct-2015)