1888     DM              massDM;  in LandauCreateJacobianMatrix()  local
1892     PetscCall(DMClone(ctx->plex[grid], &massDM));  in LandauCreateJacobianMatrix()
1893     PetscCall(DMCopyFields(ctx->plex[grid], PETSC_DETERMINE, PETSC_DETERMINE, massDM));  in LandauCreateJacobianMatrix()
1894     PetscCall(DMCreateDS(massDM));  in LandauCreateJacobianMatrix()
1895     PetscCall(DMGetDS(massDM, &prob));  in LandauCreateJacobianMatrix()
1898     PetscCall(DMCreateMatrix(massDM, &gMat));  in LandauCreateJacobianMatrix()
1903     PetscCall(DMPlexSNESComputeJacobianFEM(massDM, tvec, gMat, gMat, ctx));  in LandauCreateJacobianMatrix()
1905     PetscCall(DMDestroy(&massDM));  in LandauCreateJacobianMatrix()
2559   DM         mass_pack, massDM[LANDAU_MAX_GRIDS];  in DMPlexLandauCreateMassMatrix()  local
2575     PetscCall(DMClone(ctx->plex[grid], &massDM[grid]));  in DMPlexLandauCreateMassMatrix()
2576     PetscCall(DMCopyFields(ctx->plex[grid], PETSC_DETERMINE, PETSC_DETERMINE, massDM[grid]));  in DMPlexLandauCreateMassMatrix()
2577     PetscCall(DMCreateDS(massDM[grid]));  in DMPlexLandauCreateMassMatrix()
2578     PetscCall(DMGetDS(massDM[grid], &prob));  in DMPlexLandauCreateMassMatrix()
2584     PetscCall(DMCompositeAddDM(mass_pack, massDM[grid]));  in DMPlexLandauCreateMassMatrix()
2587       PetscCall(DMCompositeAddDM(mass_pack, massDM[grid]));  in DMPlexLandauCreateMassMatrix()
2590     PetscCall(DMCreateMatrix(massDM[grid], &subM[grid]));  in DMPlexLandauCreateMassMatrix()
2595 …cInt grid = 0; grid < ctx->num_grids; grid++) PetscCall(DMCompositeAddDM(mass_pack, massDM[grid]));  in DMPlexLandauCreateMassMatrix()
2607     DM  plex = massDM[grid];  in DMPlexLandauCreateMassMatrix()
2612     PetscCall(DMDestroy(&massDM[grid]));  in DMPlexLandauCreateMassMatrix()