- use ilu with type=ilutp instead of nofill. Contrary to luinc command ilu with nofill option doesn't not allow for partial pivoting