how to do the non-economic QR decompostion for non-square matrix
using lapack?
for example, in matlab, for a mxn full rank matrix where m>n
[Q,R] = qr(A) gives a mxm unitary matrix Q and mxn upper-triangular
matrix R. But xgeqrf returns mxn matrix Q and nxn upper-triangular matrix.
In matlab, they call this result the economic QR decomposition

