[Scipy-tickets] [SciPy] #857: spsolve doesn't work for large dimensions

SciPy scipy-tickets@scipy....
Thu Jan 29 04:54:24 CST 2009

#857: spsolve doesn't work for large dimensions
 Reporter:  sascha.baumanns      |       Owner:  wnbell            
     Type:  defect               |      Status:  new               
 Priority:  normal               |   Milestone:                    
Component:  scipy.sparse.linalg  |     Version:  devel             
 Severity:  normal               |    Keywords:  spsolve, dimension
 the sparse solver for linear systems of equations doesn't work for large
 a short example:

 from numpy import ones

 from scipy.sparse.linalg import spsolve
 from scipy.sparse import spdiags

 N = 1000000

 L = spdiags([-ones(N),2*ones(N),-ones(N)], [-1,0,1], N, N, format='csc')

 b = ones(N)
 y1 = spsolve(L,b)
 print 1

Ticket URL: <http://scipy.org/scipy/scipy/ticket/857>
SciPy <http://www.scipy.org/>
SciPy is open-source software for mathematics, science, and engineering.

More information about the Scipy-tickets mailing list