QR decomposition with modified Gram-Schmidt orthonormalization
[Q,R] = ddgqr(A)
DD matrix
returns a (m,n) orthogonal matrix Q and a (n,n) upper triangular matrix R.
A = ddrand(10,10) [Q,R] = ddgqr(A)