16   DM             celldm, swarm;  in pic_advect()  local
38 … DMDA_STENCIL_BOX, nxy, nxy, PETSC_DECIDE, PETSC_DECIDE, dof, stencil_width, NULL, NULL, &celldm));  in pic_advect()
40     PetscCall(DMDASetElementType(celldm, DMDA_ELEMENT_Q1));  in pic_advect()
42     PetscCall(DMSetFromOptions(celldm));  in pic_advect()
44     PetscCall(DMSetUp(celldm));  in pic_advect()
46     PetscCall(DMDASetUniformCoordinates(celldm, 0.0, 1.0, 0.0, 1.0, 0.0, 1.5));  in pic_advect()
64 …BoxMesh(PETSC_COMM_WORLD, dim, PETSC_TRUE, faces, NULL, NULL, PETSC_TRUE, 0, PETSC_TRUE, &celldm));  in pic_advect()
67     PetscCall(DMPlexDistribute(celldm, 0, NULL, &distributedMesh));  in pic_advect()
69       PetscCall(DMDestroy(&celldm));  in pic_advect()
70       celldm = distributedMesh;  in pic_advect()
73     PetscCall(DMSetFromOptions(celldm));  in pic_advect()
75 …PetscCall(DMPlexCreateSection(celldm, NULL, numComp, numDof, numBC, NULL, NULL, NULL, NULL, &secti…  in pic_advect()
76     PetscCall(DMSetLocalSection(celldm, section));  in pic_advect()
78     PetscCall(DMSetUp(celldm));  in pic_advect()
81     PetscCall(DMPlexComputeGeometryFVM(celldm, &cellgeom, &facegeom));  in pic_advect()
82     PetscCall(DMPlexGetMinRadius(celldm, &minradius));  in pic_advect()
96   PetscCall(DMSwarmSetCellDM(swarm, celldm));  in pic_advect()
133   PetscCall(DMCreateGlobalVector(celldm, &pfields[0]));  in pic_advect()
150       PetscCall(DMView(celldm, viewer));  in pic_advect()
156   PetscCall(DMView(celldm, PETSC_VIEWER_STDOUT_WORLD));  in pic_advect()
228         PetscCall(DMView(celldm, viewer));  in pic_advect()
235   PetscCall(DMDestroy(&celldm));  in pic_advect()