In Matlab, an ILU with threshold and pivoting (ILUTP) can be passed by default as:

setup.type = 'ilutp';
[L, U] = ilu(A, setup);

Looking for an equivalent in scipy, ex:

sparse.linalg.spilu(A, options = ??)

Tried to look up how matlab does the pivot and threshold but didn't seem to find much information.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.