19 …rrorCode CreateDM(MPI_Comm comm, ProblemData problem, MatType mat_type, VecType vec_type, DM *dm) {  in CreateDM()  argument
22   PetscCall(DMCreate(comm, dm));  in CreateDM()
23   PetscCall(DMSetType(*dm, DMPLEX));  in CreateDM()
27     PetscCall(DMSetMatrixPreallocateSkip(*dm, skip));  in CreateDM()
29   PetscCall(DMSetMatType(*dm, mat_type));  in CreateDM()
30   PetscCall(DMSetVecType(*dm, vec_type));  in CreateDM()
39   PetscCall(DMSetFromOptions(*dm));  in CreateDM()
40   PetscCall(DMViewFromOptions(*dm, NULL, "-dm_view"));  in CreateDM()
45 PetscErrorCode SetUpDM(DM dm, ProblemData problem, PetscInt degree, PetscInt q_extra, SimpleBC bc, …  in SetUpDM()  argument
49 …MSetupByOrderBegin_FEM(PETSC_TRUE, PETSC_TRUE, degree, PETSC_DECIDE, q_extra, 1, &num_comp_q, dm));  in SetUpDM()
53     PetscCall(DMGetLabel(dm, "Face Sets", &label));  in SetUpDM()
54     PetscCall(DMPlexLabelComplete(dm, label));  in SetUpDM()
65 …PetscCall(DMAddBoundary(dm, DM_BC_ESSENTIAL, name, label, num_label_values, label_values, 0, num_e…  in SetUpDM()
72       if (use_strongstg) PetscCall(SetupStrongStg(dm, bc, problem, phys));  in SetUpDM()
76   PetscCall(DMSetupByOrderEnd_FEM(PETSC_TRUE, dm));  in SetUpDM()
80   PetscCall(DMGetLocalSection(dm, &section));  in SetUpDM()
111 PetscErrorCode VizRefineDM(DM dm, User user, ProblemData problem, SimpleBC bc, Physics phys) {  in VizRefineDM()  argument
116   PetscCall(DMPlexSetRefinementUniform(dm, PETSC_TRUE));  in VizRefineDM()
118   dm_hierarchy[0] = dm;  in VizRefineDM()
128     PetscCall(DMGetVecType(dm, &vec_type));  in VizRefineDM()