New Post: MILU(0) preconditioner

Finally the fix for CSR storage! Really like what you're doing for v3 ...

But you have introduced a bug in the MILU0 code. Arrays for MSR storage have to be (nonzeros + 1) long! So

_alu = new double[ia.Length];
_jlu = new int[ia.Length];

has to be

_alu = new double[ia[n] + 1];
_jlu = new int[ia[n] + 1];

I will upload the other stuff I have later on this week.