code for ILU factorization preconditioner

[David Bateman]

Kamaraju S Kusumanchi wrote:
> David Bateman wrote:
>> Octave doesn't have the ilu function yet. It does have a very limited
>> luinc function.
>>
> 
> Looks like it is very limited. Even ILU(0) is not implemented.
> 
> [M N O] = luinc(sparse(unknown), '0');
> error: luinc: zero-level factorization not implemented
> 
> raju

Yup, it only does a drop tolerance for the partial LU factorization, as
that was easy to implement as UMFPACK supports a drop tolerance. I have
no need of an iterative solver and so never implemented this stuff. I'd
be happy to see someone else do it though :-)

D.