51   Mat         cMat, Prols[5], Rtot;  in MatCoarsenApply_MISK_private()  local
251     PetscCall(MatCreate(comm, &Prols[iterIdx]));  in MatCoarsenApply_MISK_private()
252     PetscCall(MatSetType(Prols[iterIdx], jtype));  in MatCoarsenApply_MISK_private()
253     PetscCall(MatSetSizes(Prols[iterIdx], nloc_inner, nselected, PETSC_DETERMINE, PETSC_DETERMINE));  in MatCoarsenApply_MISK_private()
254     PetscCall(MatSeqAIJSetPreallocation(Prols[iterIdx], 1, NULL));  in MatCoarsenApply_MISK_private()
255     PetscCall(MatMPIAIJSetPreallocation(Prols[iterIdx], 1, NULL, 1, NULL));  in MatCoarsenApply_MISK_private()
260       PetscCall(MatGetOwnershipRangeColumn(Prols[iterIdx], &colIndex, &Iend));  in MatCoarsenApply_MISK_private()
268           PetscCall(MatSetValues(Prols[iterIdx], 1, &fgid, 1, &colIndex, &one, INSERT_VALUES));  in MatCoarsenApply_MISK_private()
274     PetscCall(MatAssemblyBegin(Prols[iterIdx], MAT_FINAL_ASSEMBLY));  in MatCoarsenApply_MISK_private()
275     PetscCall(MatAssemblyEnd(Prols[iterIdx], MAT_FINAL_ASSEMBLY));  in MatCoarsenApply_MISK_private()
279       PetscCall(MatPtAP(cMat, Prols[iterIdx], MAT_INITIAL_MATRIX, PETSC_DETERMINE, &new_mat));  in MatCoarsenApply_MISK_private()
287   Rtot = Prols[misk - 1]; // compose P then transpose to get R  in MatCoarsenApply_MISK_private()
291     PetscCall(MatMatMult(Prols[iterIdx - 1], Rtot, MAT_INITIAL_MATRIX, PETSC_CURRENT, &P));  in MatCoarsenApply_MISK_private()
292     PetscCall(MatDestroy(&Prols[iterIdx - 1]));  in MatCoarsenApply_MISK_private()