QR decomposition
Let formula_62 be a positive semi-definite Hermitian matrix. Then it can be written as a product of its square root matrix, formula_121. Now QR decomposition can be applied to formula_122, resulting in formula_123