Somehow your suggestion made me think about using diag, Thanks :)
>>> mm(u,mm(numpy.diag(s),v))
matrix([[ 4., 2.],
[ 2., 4.]])
Now if anyone knows how to set headers for matrices or arrays?
dpn
On 5/6/07, Matthieu Brucher <[EMAIL PROTECTED]> wrote:
> Don't for
Hi,
i have two questions, both loosely related to SVD.
I've seen this post:
http://thread.gmane.org/gmane.comp.python.numeric.general/4575
>>> u,s,v = numpy.linalg.svd(numpy.array([[4,2],[2,4]])) # symmetric
matrix u == v
>>> u
array([[-0.70710678, -0.70710678],
[-0.70710678, 0.70710678]]
= matrixmultiply(matrixmultiply(u, s), v)
A minor gripe though. I'm more inclined to err on the side of clean
python rather than more operators, especially since aforementioned
GF's will help with many cases
dpn
On 3/26/07, Colin J. Williams <[EMAIL PROTECTED]> wrote:
> Cha
