[SciPy-user] Trouble with linsolve
Nils Wagner
nwagner at iam.uni-stuttgart.de
Wed Nov 29 05:06:55 EST 2006
Robert Cimrman wrote:
> Nils Wagner wrote:
>
>> 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, it was easy :) Both tickets can be resolved by:
>
> b = io.mmread('sherman2_rhs1.mtx').squeeze()
>
> in your test scripts. The rhs vector must be 1d array...
>
> I will update linsolve to do it automatically, but for now just squeeze
> the extra dimensions.
>
> cheers,
> r.
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
BTW, is there a way to store the result in a file using io.mmwrite ?
>>> io.mmwrite('sol',x)
Traceback (most recent call last):
File "<stdin>", line 1, in ?
File "/usr/lib64/python2.4/site-packages/scipy/io/mmio.py", line 244,
in mmwrite
raise ValueError, 'expected matrix'
ValueError: expected matrix
Nils
More information about the SciPy-User
mailing list