[SciPy-user] problem using umfpack.lu()

Robert Cimrman cimrman3@ntc.zcu...
Tue Dec 23 08:44:25 CST 2008


Giampaolo Cuoghi wrote:
> Robert Cimrman <cimrman3 <at> ntc.zcu.cz> writes:
> 
> 
>> Hi Giampaolo,
>>
>> Umfpack has several families of routines for single, double precisions 
>> and real or complex matrices - these can be set upon constuction of the 
>> umfpack context. Try something like:
>>
>> _family = {dtype( 'float64' ) : 'di',
>>             dtype( 'complex128' ) : 'zi'}
>>
>> family = _family[A.dtype]
>> umfpack = UmfpackContext( family = family )
>> umfpack.numeric( self.A )
>> ...
>>
>> r.
>>
> 
> Thanks a lot, Robert... now all work well! I was using the "standard context"
> that use float64 as default, now I have set the family to 'zi' to treat complex
> matrix.

You are welcome,
HTH,
r.


More information about the SciPy-user mailing list