[SciPy-user] Using umfpack to calculate a incomplete LU factorisation (ILU)
Thu Mar 8 07:19:54 CST 2007
Robert Cimrman wrote:
> Neilen Marais wrote:
>> I'm trying to get an incomplete LU factorisation of a sparse matrix using
>> scipy.linsolve.splu. It's using umfpack as far as I can tell. I'm telling
>> it to drop elements using the drop_tol keyword argument, though it seems
>> to be having no effect:
> Just one correction: linsolve.splu uses SuperLU function *gstrf, so if
> the threshold does not work, the problem is there...
> Otherwise you can get LU factors with UmfpackContext.lu, but those are
> always complete.
> SciPy-user mailing list
Just out of interest
Is there something comparable with
sigma_solve = linsolve.splu(A - sigma*B).solve
available in UMFPACK ?
More information about the SciPy-user