Fix SVD::solve with matrix r.h.s. B for non-square A (fixed in V. 3.02)


Bill Press
12-25-2007, 12:22 PM
As written, SVD::solve, in the overloaded version for matrix right-hand sides b (book page 71), requires (and tests) that the matrix A be square. This is not exactly a bug, since the subsection is named "SVD of a Square Matrix". However, it is clear in later subsections that the intent was that this routine also work for non-square A, which it doesn't.

User take2rec has kindly contributed a fixed version, at http://www.nr.com/forum/showthread.php?t=1355.

We will use this fixed version in reprintings of the book and CDROM.