Aren't there a couple of MatDenseRestoreArray missing in CPARDISO MatMatSolve? https://www.mcs.anl.gov/petsc/petsc-dev/src/mat/impls/aij/mpi/mkl_cpardiso/mkl_cpardiso.c.html#line309 Thanks, Pierre