On 7 January 2013 14:40, Giovanni Matteo Fumarola
<[hidden email]> wrote:
> I'm working on an university project about the ILU function. I've developed
> every aspect of this function
Are you hoping to include this function in core Octave? If so, can you
please write it in Octave style instead of Matlab style? Read some
Octave source code to see examples of this style. Also consult this:
thank you very much for considering the inclusion of your
implementation of ILU function in core Octave.
Of course your actual implementation (maybe rewritten in Octave style) is
fine for small matrices, and in fact you are considering to "translate" to
C language to improve the performance.
However, since the code will be used for sparse matrices, it is not just a
translation needed, but a complete rewriting using the CCS format.
Some months ago, I also tried to write the ilu(0) factorization and I
ended up with some files here enclosed. Basically, luincfull.m is your
nofill.m and luincsparse.m is the same done with vectors data, ridx, cidx
of CCS format. If you run testluincsparse.m, you discover that even in
Octave language the sparse version is faster.
To conclude, it is a translation into C of luincsparse.m and not of
luincfull.m that we need. The result would be very similar to subroutine
ilu0 in chapter 10 in Saad's book, except that we should use the CCS
format and not the CSR (if you are interested, I translated luincsparse.m
At that time, I stopped working on it and contacted Saad, who knows and