[Numpy-discussion] Re: Bug in scipy.linalg.svd?
Robert Kern
robert.kern at gmail.com
Mon Mar 27 15:20:06 CST 2006
Webb Sprague wrote:
> Hi all,
>
> scipy.__version__: 0.3.2
>
> numpy.__version__: 0.9.4
>
> OS: Gentoo.
>
> Maybe I am missing something key here, but my Numerical Recipes book says:
>
> "Any M x N matrix A whose number of rows M is greater than or equal to
> its number of columns N, can be written as the product of of an *M x
> N* column orthogonal matrix U, an N x N diagonal matrix W, ... and the
> transpose of an N x N orthogh matrix V."
This is only a convention. It is frequently called computing the "thin SVD."
Another convention is to return the fully expanded U matrix. scipy.linalg.svd
does the latter.
--
Robert Kern
robert.kern at gmail.com
"I have come to believe that the whole world is an enigma, a harmless enigma
that is made terrible by our own mad attempt to interpret it as though it had
an underlying truth."
-- Umberto Eco
More information about the Numpy-discussion
mailing list