Construct an orthonormal basis for the range of A using SVD
A : array, shape (M, N)
Q : array, shape (M, K)
Orthonormal basis for the range of A. K = effective rank of A, as determined by automatic cutoff
See also
scipy.linalg.diagsvd
scipy.linalg.cholesky