Trying to implement QR decomposition using Givens rotations, I calculate G1 to zero n32 of original matrix A, then G2 to zero n31 of G1 * A, then G2 to zero n21 of G2 * G1 * A. Residual matrix, R = G3 * G2 * G1 * A comes upper triangular as expected, so I believe my code is correct.
Looking...