Thanks.
Yes, that would be great!
Btw, I've noticed in MILU(0) you had converted to the common CSR format where the number of non-zero values is stored in an additional field of the row pointers array. I actually had planned switching to this representation for a while and finally did the change; now in master. Beside of the simpler loop boundaries, this format is also directly supported by MKL in case we want to extend the native providers to sparse routines in the future.
Thanks,
Christoph
Yes, that would be great!
Btw, I've noticed in MILU(0) you had converted to the common CSR format where the number of non-zero values is stored in an additional field of the row pointers array. I actually had planned switching to this representation for a while and finally did the change; now in master. Beside of the simpler loop boundaries, this format is also directly supported by MKL in case we want to extend the native providers to sparse routines in the future.
Thanks,
Christoph