[SciPy-user] Trouble with linsolve
Nils Wagner
nwagner at iam.uni-stuttgart.de
Wed Nov 29 04:38:02 EST 2006
Robert Cimrman wrote:
> 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.
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
Robert,
Just now I have filed two tickets. I hope someone can reproduce my results.
http://projects.scipy.org/scipy/scipy/ticket/318
http://projects.scipy.org/scipy/scipy/ticket/319
I have tried SuperLU on 32 and 64 bit machines - the segfault persists.
I have installed UMFPACK on a 64 bit machine.
Any pointer would be appreciated.
Thanks in advance.
Nils
More information about the SciPy-User
mailing list