Classical Gram-Schmidt

Classical Gram-Schmidt computes all orthogonal projection coefficients $⟨v_j,q_i⟩$ from the original input $v_j$, then subtracts their sum as in the displayed Gram-Schmidt process formula.