[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