Used for solving systems and iterative eigenvalue algorithms. A=QRA = QRA=QR Here: QQQ orthonormal (via Gram–Schmidt) RRR upper triangular.