SVD for an arbitrary matrix a:
A[aLo1..aHi1, aLo2..aHi2], m = A.Dim1, n = A.Dim2, nm = min(n,m)
fullmat = false (default)
A = U * Q * VT UT * A * V = Q
m x n m x nm nm x nm nm x n nm x m m x n n x nm nm x nm
U[aLo1..aHi1, 1..nm], V[aLo2..aHi2, 1..nm2], Q[1..nm]
fullmat = true
A = U * Q * VT UT * A * V = Q
m x n m x m m x n n x n m x m m x n n x n m x n
U[aLo1..aHi1, 1..m], V[aLo2..aHi2, 1..n]
Q[1, nm] - first nm values