Call PETSc matrix operator MatDiagonalSet (…)