QR decomposition
perform two iterations of the Givens rotation (note that the Givens rotation algorithm used here differs slightly from above) to yield an upper triangular matrix in order to compute the QR decomposition.