This works. Thanks. Does the Cholesky/ICC factorization have a specialized storage i.e. like the ILU storage? -Paul > Insert it as the upper triangular part in a MATSBAIJ, then everything > you want will work.