[SciPy-dev] Ticket #552 - no errors with recent versions

josef.pktd at gmail.com josef.pktd at gmail.com
Wed Aug 6 13:21:07 EDT 2008


I don't have any problem with newer versions of numpy with ticket 552 on
WindowsXP


Ticket #552 (new defect): linalg.svd fails for some matrices [Win32, SciPy
0.6.0, NumPy 1.0.4]


Python 2.4.3 (#69, Mar 29 2006, 17:35:34) [MSC v.1310 32 bit (Intel)] on
win32
Type "help", "copyright", "credits" or "license" for more information.
>>> import scipy
>>> execfile(r'cpos_failure_mintestcase.py')
H:
[[ -8.19200144e+03   9.58422857e+00  -2.18886768e+00]
 [  2.70299079e+00   2.32063587e+03   5.67734088e+00]
 [  0.00000000e+00   0.00000000e+00   0.00000000e+00]]
Now computing SVD of H
Computed SVD of H. Result:
(matrix([[ -1.00000000e+00,  -1.39611407e-06,   0.00000000e+00],
        [ -1.39611407e-06,   1.00000000e+00,   0.00000000e+00],
        [  0.00000000e+00,   0.00000000e+00,   1.00000000e+00]]), array([
8192.0
0734196,  2320.64438555,     0.        ]), matrix([[  9.99999279e-01,
-1.170344
22e-03,   2.67194555e-04],
        [  1.16968708e-03,   9.99996323e-01,   2.44645150e-03],
        [ -2.70056763e-04,  -2.44613720e-03,   9.99996972e-01]]))

>>> scipy.version.version
'0.6.0'
>>> import numpy
>>> numpy.version.version
'1.1.0'
>>>


Python 2.5.2 (r252:60911, Feb 21 2008, 13:11:45) [MSC v.1310 32 bit (Intel)]
on
win32
Type "help", "copyright", "credits" or "license" for more information.
>>> import scipy
>>> execfile(r'cpos_failure_mintestcase.py')
H:
[[ -8.192e+03   9.584e+00  -2.189e+00]
 [  2.703e+00   2.321e+03   5.677e+00]
 [  0.000e+00   0.000e+00   0.000e+00]]
Now computing SVD of H
Computed SVD of H. Result:
(matrix([[ -1.000e+00,  -1.396e-06,   0.000e+00],
        [ -1.396e-06,   1.000e+00,   0.000e+00],
        [  0.000e+00,   0.000e+00,   1.000e+00]]), array([ 8192.007,
2320.644,
    0.   ]), matrix([[  1.000e+00,  -1.170e-03,   2.672e-04],
        [  1.170e-03,   1.000e+00,   2.446e-03],
        [ -2.701e-04,  -2.446e-03,   1.000e+00]]))
>>> scipy.version.version
'0.7.0.dev'
>>> import numpy
>>> numpy.version.version
'1.2.0.dev5608'
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20080806/f34e3900/attachment.html>


More information about the SciPy-Dev mailing list