[SciPy-user] Using umfpack to calculate a incomplete LU factorisation (ILU)

Robert Cimrman cimrman3 at ntc.zcu.cz
Wed Mar 7 04:45:01 EST 2007


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.

cheers,
r.



More information about the SciPy-User mailing list