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

Giampaolo Cuoghi jpaul74@gmail....
Tue Dec 23 08:40:02 CST 2008


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.

Giampaolo





More information about the SciPy-user mailing list