[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: code for ILU factorization preconditioner
From: |
David Bateman |
Subject: |
Re: code for ILU factorization preconditioner |
Date: |
Sun, 23 Dec 2007 14:06:10 +0100 |
User-agent: |
Thunderbird 1.5.0.7 (X11/20060921) |
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.