[SciPy-user] Trouble with linsolve

Robert Cimrman cimrman3 at ntc.zcu.cz
Wed Nov 29 03:29:11 CST 2006


Nils Wagner wrote:
> With UMFPACK 4.4 and a dense RHS I get
> 
>     x = spsolve(K_dyn, f)
>   File "/usr/lib64/python2.4/site-packages/scipy/linsolve/linsolve.py",
> line 65, in spsolve
>     return umf.linsolve( umfpack.UMFPACK_A, mat, b, autoTranspose = True )
>   File
> "/usr/lib64/python2.4/site-packages/scipy/linsolve/umfpack/umfpack.py",
> line 566, in linsolve
>     sol = self.solve( sys, mtx, rhs, autoTranspose )
>   File
> "/usr/lib64/python2.4/site-packages/scipy/linsolve/umfpack/umfpack.py",
> line 508, in solve
>     self._numeric, self.control, self.info )
>   File
> "/usr/lib64/python2.4/site-packages/scipy/linsolve/umfpack/_umfpack.py",
> line 214, in umfpack_di_solve
>     return __umfpack.umfpack_di_solve(*args)
> ValueError: object too deep for desired array
>>>> K_dyn
> <71987x71987 sparse matrix of type '<type 'numpy.float64'>'
>         with 3083884 stored elements (space for 3083884)
>         in Compressed Sparse Column format>
> 
> Any idea ?

I use umfpack on a 32 bit system regularly without problems, so it may 
be a 64-bit issue. I have yet to set up a 64-bit installation of scipy 
to test it, though, so I cannot tell you more immediately.

Did you also check your inputs to spsolve?

r.


More information about the SciPy-user mailing list