Excellent, thanks! Yes, I'll include it into the library.
I've noticed that the existing ILU(0) and ILUTP implementations do not actually leverage the sparse storage format at all yet, which may explain at least a major part of their slowness (but of course has no effect on the number of iterations).
Thanks,
Christoph
I've noticed that the existing ILU(0) and ILUTP implementations do not actually leverage the sparse storage format at all yet, which may explain at least a major part of their slowness (but of course has no effect on the number of iterations).
Thanks,
Christoph