Project vector a onto vector b: proj_b(a) = (a·b / b·b)·b. Also computes the orthogonal component and verifies that projection + orthogonal = a.