Understood. But, the petsc ILU storage for L and U is non traditional In ILU, all of L is stored and then U is stored. For U, the diagonal element of a row in U is placed at the end of each row (to get better cache performance in the solve).
I'm wondering if anything similar happens in Cholesky where the diagonals are placed at the end of each row. -Paul > Cholesky/ICC only stores upper half (U) matrix, not L. > ILU stores L and U. > > Hong >> -Paul >> >>> Insert it as the upper triangular part in a MATSBAIJ, then everything you >>> want will work. >>
