[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