[SciPy-user] linsolve.factorized was: Re: Using umfpack to calculate a incomplete LU factorisation (ILU)
Thu Mar 8 12:17:32 CST 2007
On Thu, 08 Mar 2007 18:10:30 +0100
Robert Cimrman <email@example.com> wrote:
> 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,
>> 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 :-)
> In :scipy.linsolve.factorized?
> Definition: scipy.linsolve.factorized(A)
> Return a fuction for solving a linear system, with A
> 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.
> SciPy-user mailing list
Great ! Thank you very much !
More information about the SciPy-user