2563   Mat        packM, subM[LANDAU_MAX_GRIDS];  in DMPlexLandauCreateMassMatrix()  local
2599   PetscCall(DMCreateMatrix(mass_pack, &packM));  in DMPlexLandauCreateMassMatrix()
2601   PetscCall(MatSetOption(packM, MAT_STRUCTURALLY_SYMMETRIC, PETSC_TRUE));  in DMPlexLandauCreateMassMatrix()
2602   PetscCall(MatSetOption(packM, MAT_IGNORE_ZERO_ENTRIES, PETSC_TRUE));  in DMPlexLandauCreateMassMatrix()
2615   PetscCall(MatGetSize(packM, &N2, NULL));  in DMPlexLandauCreateMassMatrix()
2638         PetscCall(MatSetValues(packM, 1, &row, nzl, colbuf, vals, INSERT_VALUES));  in DMPlexLandauCreateMassMatrix()
2646   PetscCall(MatAssemblyBegin(packM, MAT_FINAL_ASSEMBLY));  in DMPlexLandauCreateMassMatrix()
2647   PetscCall(MatAssemblyEnd(packM, MAT_FINAL_ASSEMBLY));  in DMPlexLandauCreateMassMatrix()
2648   PetscCall(PetscObjectSetName((PetscObject)packM, "mass"));  in DMPlexLandauCreateMassMatrix()
2649   PetscCall(MatViewFromOptions(packM, NULL, "-dm_landau_mass_view"));  in DMPlexLandauCreateMassMatrix()
2650   ctx->M = packM;  in DMPlexLandauCreateMassMatrix()
2651   if (Amat) *Amat = packM;  in DMPlexLandauCreateMassMatrix()