[SciPy-dev] Possibly bug in Sparse.py

[Robert Kern]

Nils Wagner wrote:

> Traceback (most recent call last):
>  File "meerbergen.py", line 35, in ?
>    y, info = linalg.gmres(omega[i]*B-A,v)
>  File "/usr/lib/python2.3/site-packages/scipy/linalg/iterative.py", line 
> 583, in gmres
>    work = sb.zeros((6+restrt)*n,typ)
> MemoryError: can't allocate memory for array

Well, since most everything seems to operate with CSC format matrices, 
you could try explicitly converting your arrays to CSC format at the 
beginning and delete the COO instances.

-- 
Robert Kern
rkern at ucsd.edu

"In the fields of hell where the grass grows high
  Are the graves of dreams allowed to die."
   -- Richard Harter