[SciPy-user] bug in diagsvd
Carlo Mattoni
carlo at jsbach.harvard.edu
Tue Jun 1 10:34:34 EDT 2004
I'm not sure if this is the best place to report a bug, but I have
found one in the function diagsvd in decomp.py.
Python 2.3 (#1, Sep 13 2003, 00:49:11)
[GCC 3.3 20030304 (Apple Computer, Inc. build 1495)] on darwin
Type "help", "copyright", "credits" or "license" for more information.
>>> import scipy
>>> scipy.linalg.diagsvd( scipy.array( [ 1, 0, 0 ] ), 3, 3 )
Traceback (most recent call last):
File "<stdin>", line 1, in ?
File "/System/Library/Frameworks/Python.framework/Versions/2.3/lib/python2.3/site-packages/scipy/linalg/decomp.py",
line 297, in diagsvd
part = diag(s)
NameError: global name 'diag' is not defined
It appears that the problem is that the definition of diagsvd does not
properly prepend scipy function calls with scip_base as all of the other
functions in decomp.py do.
I believe the diagsvd function (starting on line 295) should be changed
from:
def diagsvd(s,M,N):
"""Return sigma from singular values and original size M,N."""
part = diag(s)
typ = part.typecode()
MorN = len(s)
if MorN == M:
return c_[part,zeros((M,N-M),typ)]
elif MorN == N:
return r_[part,zeros((M-N,N),typ)]
else:
raise ValueError, "Length of s must be M or N."
to:
def diagsvd(s,M,N):
"""Return sigma from singular values and original size M,N."""
part = scipy_base.diag(s)
typ = part.typecode()
MorN = len(s)
if MorN == M:
return c_[part,scipy_base.zeros((M,N-M),typ)]
elif MorN == N:
return r_[part,scipy_base.zeros((M-N,N),typ)]
else:
raise ValueError, "Length of s must be M or N."
Once this change is made, the original test works as expected:
Python 2.3 (#1, Sep 13 2003, 00:49:11)
[GCC 3.3 20030304 (Apple Computer, Inc. build 1495)] on darwin
Type "help", "copyright", "credits" or "license" for more information.
>>> import scipy
>>> scipy.linalg.diagsvd( scipy.array( [ 1, 0, 0 ] ), 3, 3 )
array([[ 1., 0., 0.],
[ 0., 0., 0.],
[ 0., 0., 0.]])
More information about the SciPy-User
mailing list