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

[Giampaolo Cuoghi]

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