def rq(A): from scipy.linalg import qr Q,R = qr(flipud(A).T) R = flipud(R.T) Q = Q.T return R[:,::-1],Q[::-1,:]