25 Jul
2006
25 Jul
'06
4:45 p.m.
Hi, there was a thread about this before, diag() is currently only partly useful if you work with numpy-matrices, because the 1d->2d direction doesn't work, as there are no 1d-numpy-matrices. This is unfortunate because a numpy-matrix with shape (n,1) or (1,m) should be naturally treated as a vector, imho. So it would be nice if this could be fixed. It's probably not the most efficient solution, but what I want for numpy-matrix input x is to get: mat(diag(x.A.squeeze)) where diag is the current implementation. This means that if x is not a vector ("truly 2d"), then nothing is changed. But if one of the dimensions of x is ==1, then it's turned into a 1d-array, and diag works as it should. Does that sound reasonable? Thanks, Sven
6445
Age (days ago)
6457
Last active (days ago)
13 comments
6 participants
participants (6)
-
Alan G Isaac
-
Charles R Harris
-
Christopher Barker
-
Keith Goodman
-
Robert Kern
-
Sven Schreiber