[SciPy-user] linsolve.factorized was: Re: Using umfpack to calculate a incomplete LU factorisation (ILU)
Thu Mar 8 11:10:30 CST 2007
Robert Cimrman wrote:
> Nils Wagner wrote:
>> Just out of interest
>> Is there something comparable with
>> sigma_solve = linsolve.splu(A - sigma*B).solve
>> available in UMFPACK ?
> Sure, although it is not a oneliner. This is the relevant part of the
> Would you be interested in exposing this in linsolve, as, say,
> def factorized( A ):
> if isUmfpack and useUmfpack:
> # This must be written...
> return splu( A ).solve
Well, I did it since I am going to need this, too :-)
Return a fuction for solving a linear system, with A pre-factorized.
solve = factorized( A ) # Makes LU decomposition.
x1 = solve( rhs1 ) # Uses the LU factors.
x2 = solve( rhs2 ) # Uses again the LU factors.
This uses UMFPACK if available.
More information about the SciPy-user