Fix a typo in documentation of fs_svd()
Fix reordering of singular vectors of a non-square matrix

so that in the SVD, m == u * sigma * vh, sigma(i,j) == 0 if i != j.
Generalize documentation of SVD to rectangular matrices
Add diagonalization procedures that estimate error bounds