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