[NumPy-Tickets] [NumPy] #1655: Poor precision using linalg.solve for certain matrices

NumPy Trac numpy-tickets@scipy....
Fri Oct 29 21:23:21 CDT 2010

#1655: Poor precision using linalg.solve for certain matrices
 Reporter:  zoof          |       Owner:  pv   
     Type:  defect        |      Status:  new  
 Priority:  normal        |   Milestone:  2.0.0
Component:  numpy.linalg  |     Version:  1.5.0
 Keywords:                |  
 There appears to be a bug in either the linear system solver or in the way
 than linalg computes LU decompositions.  On a number of different
 distributions, using somewhat different versions of numpy, I get much
 better precision by inverting the matrix (linalg.inv) and multiplying than
 by either linalg.solve alone or linalg.cholesky in combination with
 linalg.solve.  An example has been posted to: http://pastebin.com/KpXdTVh7

 numpy 1.3.0 on Ubuntu Lucid
 numpy 1.5.0 on Arch
 numpy 1.5.1rc1 on Centos

Ticket URL: <http://projects.scipy.org/numpy/ticket/1655>
NumPy <http://projects.scipy.org/numpy>
My example project

More information about the NumPy-Tickets mailing list