15 …cErrorCode CreateDM(Honee honee, ProblemData problem, MatType mat_type, VecType vec_type, DM *dm) {  in CreateDM()  argument
26     PetscCall(DMPlexCreateFromFile(comm, filename, "HONEE", PETSC_TRUE, dm));  in CreateDM()
29     PetscCall(DMCreate(comm, dm));  in CreateDM()
30     PetscCall(DMSetType(*dm, DMPLEX));  in CreateDM()
36     PetscCall(DMSetMatrixPreallocateSkip(*dm, skip));  in CreateDM()
38   PetscCall(DMSetMatType(*dm, mat_type));  in CreateDM()
39   PetscCall(DMSetVecType(*dm, vec_type));  in CreateDM()
45     PetscCall(DMPlexDistributeGetDefault(*dm, &distribute));  in CreateDM()
49   PetscCall(DMSetSparseLocalize(*dm, PETSC_FALSE));  in CreateDM()
50   PetscCall(DMSetBlockingType(*dm, DM_BLOCKING_FIELD_NODE));  in CreateDM()
51   PetscCall(DMPlexSetPartitionBalance(*dm, PETSC_TRUE));  in CreateDM()
54   PetscCall(DMSetFromOptions(*dm));  in CreateDM()
55   PetscCall(DMViewFromOptions(*dm, NULL, "-dm_view"));  in CreateDM()
59     PetscCall(DMGetNumFields(*dm, &num_fields));  in CreateDM()
62       PetscCall(DMGetGlobalSection(*dm, &dummy));  in CreateDM()
66   PetscCall(DMPlexSetScale(*dm, PETSC_UNIT_LENGTH, honee->units->meter));  in CreateDM()
67   PetscCall(DMPlexSetScale(*dm, PETSC_UNIT_MASS, honee->units->kilogram));  in CreateDM()
68   PetscCall(DMPlexSetScale(*dm, PETSC_UNIT_TIME, honee->units->second));  in CreateDM()
69   PetscCall(DMPlexSetScale(*dm, PETSC_UNIT_TEMPERATURE, honee->units->Kelvin));  in CreateDM()
74 PetscErrorCode SetUpDM(DM dm, ProblemData problem, PetscInt degree, PetscInt q_extra, Physics phys)…  in SetUpDM()  argument
78   PetscCall(DMClearFields(dm));  in SetUpDM()
79   PetscCall(DMSetLocalSection(dm, NULL));  in SetUpDM()
80 …MSetupByOrderBegin_FEM(PETSC_TRUE, PETSC_TRUE, degree, PETSC_DECIDE, q_extra, 1, &num_comp_q, dm));  in SetUpDM()
85     if (use_strongstg) PetscCall(SetupStrongStg(dm, problem, phys));  in SetUpDM()
90     PetscCall(DMGetLabel(dm, "Face Sets", &label));  in SetUpDM()
91 …PetscCheck(label, PetscObjectComm((PetscObject)dm), PETSC_ERR_SUP, "DM does not have 'Face Sets' l…  in SetUpDM()
92     PetscCall(DMPlexLabelComplete(dm, label));  in SetUpDM()
100       PetscCall(BCDefinitionSetDM(problem->bc_defs[i], dm));  in SetUpDM()
104 …PetscCall(DMAddBoundary(dm, DM_BC_ESSENTIAL, name, label, num_label_values, label_values, 0, num_e…  in SetUpDM()
110   PetscCall(DMSetupByOrderEnd_FEM(PETSC_TRUE, dm));  in SetUpDM()
114   PetscCall(DMGetLocalSection(dm, &section));  in SetUpDM()
123 PetscErrorCode VizRefineDM(DM dm, Honee honee, ProblemData problem, Physics phys) {  in VizRefineDM()  argument
128   PetscCall(DMPlexSetRefinementUniform(dm, PETSC_TRUE));  in VizRefineDM()
130   dm_hierarchy[0] = dm;  in VizRefineDM()
140     PetscCall(DMGetVecType(dm, &vec_type));  in VizRefineDM()