QR decomposition

Given a matrix A, its QR decomposition refers to a factorisation of A

A = QR,

where Q is a orthogonal matrix (QQ^T = I), and R is a upper triangular matrix.