returns pair of orthogonal matrix (Q) and upper right triangular matrix (R) factoring full column rank matrix mat.