The following code:
from scipy import * from scipy import linalg
m = matrix( [ [1,1,0,0], [1,1,0,0], [0,0,1,1], [0,0,1,1] ] )
u,s,v = linalg.svd( m )
fails with the following message:
Traceback (most recent call last): File "boo.py", line 10, in <module> u,s,v = linalg.svd( m ) File "/usr/lib64/python2.6/site-packages/scipy/linalg/decomp.py", line 509, in svd lwork = calc_lwork.gesdd(gesdd.prefix,m,n,compute_uv)[1] RuntimeError: more argument specifiers than keyword list entries (remaining format:'|:calc_lwork.gesdd')
On the other hand, calculating la, v = eig( m ) works just fine.
If I see this correctly, my SciPy version is 0.6.0; running on 64bit Suse 11.
Any thoughts?
Best,
Ph.
Hi Philipp
2009/5/5 Philipp K. Janert python@beyondcode.org:
If I see this correctly, my SciPy version is 0.6.0; running on 64bit Suse 11.
SciPy 0.6 is quite old, and it is likely that the problem was fixed in the mean time.
On SciPy 0.7 I see:
In [31]: u,s,v = linalg.svd(m)
In [32]: u Out[32]: array([[-0.70710678, 0. , -0.70710678, 0. ], [-0.70710678, 0. , 0.70710678, 0. ], [ 0. , -0.70710678, 0. , -0.70710678], [ 0. , -0.70710678, 0. , 0.70710678]])
In [33]: s Out[33]: array([ 2., 2., 0., 0.])
In [34]: v Out[34]: array([[-0.70710678, -0.70710678, -0. , -0. ], [-0. , -0. , -0.70710678, -0.70710678], [-0.70710678, 0.70710678, 0. , 0. ], [ 0. , 0. , -0.70710678, 0.70710678]])
Regards Stéfan
Thanks for the quick reply. I'll try upgrading.
Best,
Ph.
On Monday 04 May 2009 04:21:29 pm Stéfan van der Walt wrote:
Hi Philipp
2009/5/5 Philipp K. Janert python@beyondcode.org:
If I see this correctly, my SciPy version is 0.6.0; running on 64bit Suse 11.
SciPy 0.6 is quite old, and it is likely that the problem was fixed in the mean time.
On SciPy 0.7 I see:
In [31]: u,s,v = linalg.svd(m)
In [32]: u Out[32]: array([[-0.70710678, 0. , -0.70710678, 0. ], [-0.70710678, 0. , 0.70710678, 0. ], [ 0. , -0.70710678, 0. , -0.70710678], [ 0. , -0.70710678, 0. , 0.70710678]])
In [33]: s Out[33]: array([ 2., 2., 0., 0.])
In [34]: v Out[34]: array([[-0.70710678, -0.70710678, -0. , -0. ], [-0. , -0. , -0.70710678, -0.70710678], [-0.70710678, 0.70710678, 0. , 0. ], [ 0. , 0. , -0.70710678, 0.70710678]])
Regards Stéfan _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
participants (2)
-
Philipp K. Janert -
Stéfan van der Walt