6 Jan
2010
6 Jan
'10
8:55 p.m.
SV-decomposition of matrix A (m x n) produces matrices U (m x n) and V (n x n) and a vector s representing the diagonal matrix S (n x n) so that A = U S V'. The solution of the equation Ax = b, x, is of dimension n that is equal to (dim0 V) rather than (dim0 U) as is written in SV-solve for its optional parameter. Hence the patch. Obviously, bug appears only for non-rectanular A, so it passed the tests.