xref: /petsc/src/dm/impls/plex/tests/ex1.c (revision 2f613bf53f46f9356e00a2ca2bd69453be72fc31)
1 static char help[] = "Tests various DMPlex routines to construct, refine and distribute a mesh.\n\n";
2 
3 #include <petscdmplex.h>
4 #include <petscdmplextransform.h>
5 #include <petscsf.h>
6 
7 enum {STAGE_LOAD, STAGE_DISTRIBUTE, STAGE_REFINE, STAGE_OVERLAP};
8 
9 typedef struct {
10   PetscLogEvent createMeshEvent;
11   PetscLogStage stages[4];
12   /* Domain and mesh definition */
13   PetscInt      dim;                             /* The topological mesh dimension */
14   PetscInt      overlap;                         /* The cell overlap to use during partitioning */
15   PetscBool     testp4est[2];
16   PetscBool     redistribute;
17   PetscBool     final_ref;                       /* Run refinement at the end */
18   PetscBool     final_diagnostics;               /* Run diagnostics on the final mesh */
19 } AppCtx;
20 
21 PetscErrorCode ProcessOptions(MPI_Comm comm, AppCtx *options)
22 {
23   PetscErrorCode ierr;
24 
25   PetscFunctionBegin;
26   options->dim               = 2;
27   options->overlap           = 0;
28   options->testp4est[0]      = PETSC_FALSE;
29   options->testp4est[1]      = PETSC_FALSE;
30   options->redistribute      = PETSC_FALSE;
31   options->final_ref         = PETSC_FALSE;
32   options->final_diagnostics = PETSC_TRUE;
33 
34   ierr = PetscOptionsBegin(comm, "", "Meshing Problem Options", "DMPLEX");CHKERRQ(ierr);
35   ierr = PetscOptionsRangeInt("-dim", "The topological mesh dimension", "ex1.c", options->dim, &options->dim, NULL,1,3);CHKERRQ(ierr);
36   ierr = PetscOptionsBoundedInt("-overlap", "The cell overlap for partitioning", "ex1.c", options->overlap, &options->overlap, NULL,0);CHKERRQ(ierr);
37   ierr = PetscOptionsBool("-test_p4est_seq", "Test p4est with sequential base DM", "ex1.c", options->testp4est[0], &options->testp4est[0], NULL);CHKERRQ(ierr);
38   ierr = PetscOptionsBool("-test_p4est_par", "Test p4est with parallel base DM", "ex1.c", options->testp4est[1], &options->testp4est[1], NULL);CHKERRQ(ierr);
39   ierr = PetscOptionsBool("-test_redistribute", "Test redistribution", "ex1.c", options->redistribute, &options->redistribute, NULL);CHKERRQ(ierr);
40   ierr = PetscOptionsBool("-final_ref", "Run uniform refinement on the final mesh", "ex1.c", options->final_ref, &options->final_ref, NULL);CHKERRQ(ierr);
41   ierr = PetscOptionsBool("-final_diagnostics", "Run diagnostics on the final mesh", "ex1.c", options->final_diagnostics, &options->final_diagnostics, NULL);CHKERRQ(ierr);
42   ierr = PetscOptionsEnd();CHKERRQ(ierr);
43 
44   ierr = PetscLogEventRegister("CreateMesh", DM_CLASSID, &options->createMeshEvent);CHKERRQ(ierr);
45   ierr = PetscLogStageRegister("MeshLoad",       &options->stages[STAGE_LOAD]);CHKERRQ(ierr);
46   ierr = PetscLogStageRegister("MeshDistribute", &options->stages[STAGE_DISTRIBUTE]);CHKERRQ(ierr);
47   ierr = PetscLogStageRegister("MeshRefine",     &options->stages[STAGE_REFINE]);CHKERRQ(ierr);
48   ierr = PetscLogStageRegister("MeshOverlap",    &options->stages[STAGE_OVERLAP]);CHKERRQ(ierr);
49   PetscFunctionReturn(0);
50 }
51 
52 PetscErrorCode CreateMesh(MPI_Comm comm, AppCtx *user, DM *dm)
53 {
54   PetscInt       dim           = user->dim;
55   PetscBool      testp4est_seq = user->testp4est[0];
56   PetscBool      testp4est_par = user->testp4est[1];
57   PetscMPIInt    rank, size;
58   PetscBool      periodic;
59   PetscErrorCode ierr;
60 
61   PetscFunctionBegin;
62   ierr = PetscLogEventBegin(user->createMeshEvent,0,0,0,0);CHKERRQ(ierr);
63   ierr = MPI_Comm_rank(comm, &rank);CHKERRMPI(ierr);
64   ierr = MPI_Comm_size(comm, &size);CHKERRMPI(ierr);
65   ierr = PetscLogStagePush(user->stages[STAGE_LOAD]);CHKERRQ(ierr);
66   ierr = DMCreate(comm, dm);CHKERRQ(ierr);
67   ierr = DMSetType(*dm, DMPLEX);CHKERRQ(ierr);
68   ierr = DMSetFromOptions(*dm);CHKERRQ(ierr);
69 
70   /* For topologically periodic meshes, we first localize coordinates,
71      and then remove any information related with the
72      automatic computation of localized vertices.
73      This way, refinement operations and conversions to p4est
74      will preserve the shape of the domain in physical space */
75   ierr = DMLocalizeCoordinates(*dm);CHKERRQ(ierr);
76   ierr = DMGetPeriodicity(*dm, &periodic, NULL, NULL, NULL);CHKERRQ(ierr);
77   if (periodic) {ierr = DMSetPeriodicity(*dm, PETSC_TRUE, NULL, NULL, NULL);CHKERRQ(ierr);}
78 
79   ierr = DMViewFromOptions(*dm,NULL,"-init_dm_view");CHKERRQ(ierr);
80   ierr = DMGetDimension(*dm, &dim);CHKERRQ(ierr);
81 
82   if (testp4est_seq) {
83 #if defined(PETSC_HAVE_P4EST)
84     DM dmConv = NULL;
85 
86     ierr = DMPlexCheckSymmetry(*dm);CHKERRQ(ierr);
87     ierr = DMPlexCheckSkeleton(*dm, 0);CHKERRQ(ierr);
88     ierr = DMPlexCheckFaces(*dm, 0);CHKERRQ(ierr);
89     ierr = DMPlexCheckGeometry(*dm);CHKERRQ(ierr);
90     ierr = DMPlexCheckPointSF(*dm);CHKERRQ(ierr);
91     ierr = DMPlexCheckInterfaceCones(*dm);CHKERRQ(ierr);
92     ierr = DMPlexSetRefinementUniform(*dm, PETSC_TRUE);CHKERRQ(ierr);
93     ierr = DMPlexSetTransformType(*dm, DMPLEXREFINETOBOX);CHKERRQ(ierr);
94     ierr = DMRefine(*dm, PETSC_COMM_WORLD, &dmConv);CHKERRQ(ierr);
95     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, NULL);CHKERRQ(ierr);
96     if (dmConv) {
97       ierr = DMDestroy(dm);CHKERRQ(ierr);
98       *dm  = dmConv;
99     }
100     ierr = DMViewFromOptions(*dm,NULL,"-initref_dm_view");CHKERRQ(ierr);
101     ierr = DMPlexCheckSymmetry(*dm);CHKERRQ(ierr);
102     ierr = DMPlexCheckSkeleton(*dm, 0);CHKERRQ(ierr);
103     ierr = DMPlexCheckFaces(*dm, 0);CHKERRQ(ierr);
104     ierr = DMPlexCheckGeometry(*dm);CHKERRQ(ierr);
105     ierr = DMPlexCheckPointSF(*dm);CHKERRQ(ierr);
106     ierr = DMPlexCheckInterfaceCones(*dm);CHKERRQ(ierr);
107 
108     ierr = DMConvert(*dm,dim == 2 ? DMP4EST : DMP8EST,&dmConv);CHKERRQ(ierr);
109     if (dmConv) {
110       ierr = PetscObjectSetOptionsPrefix((PetscObject) dmConv, "conv_seq_1_");CHKERRQ(ierr);
111       ierr = DMSetFromOptions(dmConv);CHKERRQ(ierr);
112       ierr = DMDestroy(dm);CHKERRQ(ierr);
113       *dm  = dmConv;
114     }
115     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, "conv_seq_1_");CHKERRQ(ierr);
116     ierr = DMSetUp(*dm);CHKERRQ(ierr);
117     ierr = DMViewFromOptions(*dm, NULL, "-dm_view");CHKERRQ(ierr);
118     ierr = DMConvert(*dm,DMPLEX,&dmConv);CHKERRQ(ierr);
119     if (dmConv) {
120       ierr = PetscObjectSetOptionsPrefix((PetscObject) dmConv, "conv_seq_2_");CHKERRQ(ierr);
121       ierr = DMSetFromOptions(dmConv);CHKERRQ(ierr);
122       ierr = DMDestroy(dm);CHKERRQ(ierr);
123       *dm  = dmConv;
124     }
125     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, "conv_seq_2_");CHKERRQ(ierr);
126     ierr = DMViewFromOptions(*dm, NULL, "-dm_view");CHKERRQ(ierr);
127     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, NULL);CHKERRQ(ierr);
128 #else
129     SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_SUP,"Recompile with --download-p4est");
130 #endif
131   }
132 
133   ierr = PetscLogStagePop();CHKERRQ(ierr);
134   if (!testp4est_seq) {
135     ierr = PetscLogStagePush(user->stages[STAGE_DISTRIBUTE]);CHKERRQ(ierr);
136     ierr = DMViewFromOptions(*dm, NULL, "-dm_pre_dist_view");CHKERRQ(ierr);
137     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, "dist_");CHKERRQ(ierr);
138     ierr = DMSetFromOptions(*dm);CHKERRQ(ierr);
139     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, NULL);CHKERRQ(ierr);
140     ierr = PetscLogStagePop();CHKERRQ(ierr);
141     ierr = DMViewFromOptions(*dm, NULL, "-distributed_dm_view");CHKERRQ(ierr);
142   }
143   ierr = PetscLogStagePush(user->stages[STAGE_REFINE]);CHKERRQ(ierr);
144   ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, "ref_");CHKERRQ(ierr);
145   ierr = DMSetFromOptions(*dm);CHKERRQ(ierr);
146   ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, NULL);CHKERRQ(ierr);
147   ierr = PetscLogStagePop();CHKERRQ(ierr);
148 
149   if (testp4est_par) {
150 #if defined(PETSC_HAVE_P4EST)
151     DM dmConv = NULL;
152 
153     ierr = DMViewFromOptions(*dm, NULL, "-dm_tobox_view");CHKERRQ(ierr);
154     ierr = DMPlexSetRefinementUniform(*dm, PETSC_TRUE);CHKERRQ(ierr);
155     ierr = DMPlexSetTransformType(*dm, DMPLEXREFINETOBOX);CHKERRQ(ierr);
156     ierr = DMRefine(*dm, PETSC_COMM_WORLD, &dmConv);CHKERRQ(ierr);
157     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, NULL);CHKERRQ(ierr);
158     if (dmConv) {
159       ierr = DMDestroy(dm);CHKERRQ(ierr);
160       *dm  = dmConv;
161     }
162     ierr = DMViewFromOptions(*dm, NULL, "-dm_tobox_view");CHKERRQ(ierr);
163     ierr = DMPlexCheckSymmetry(*dm);CHKERRQ(ierr);
164     ierr = DMPlexCheckSkeleton(*dm, 0);CHKERRQ(ierr);
165     ierr = DMPlexCheckFaces(*dm, 0);CHKERRQ(ierr);
166     ierr = DMPlexCheckGeometry(*dm);CHKERRQ(ierr);
167     ierr = DMPlexCheckPointSF(*dm);CHKERRQ(ierr);
168     ierr = DMPlexCheckInterfaceCones(*dm);CHKERRQ(ierr);
169 
170     ierr = DMConvert(*dm,dim == 2 ? DMP4EST : DMP8EST,&dmConv);CHKERRQ(ierr);
171     if (dmConv) {
172       ierr = PetscObjectSetOptionsPrefix((PetscObject) dmConv, "conv_par_1_");CHKERRQ(ierr);
173       ierr = DMSetFromOptions(dmConv);CHKERRQ(ierr);
174       ierr = DMDestroy(dm);CHKERRQ(ierr);
175       *dm  = dmConv;
176     }
177     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, "conv_par_1_");CHKERRQ(ierr);
178     ierr = DMSetUp(*dm);CHKERRQ(ierr);
179     ierr = DMViewFromOptions(*dm, NULL, "-dm_view");CHKERRQ(ierr);
180     ierr = DMConvert(*dm, DMPLEX, &dmConv);CHKERRQ(ierr);
181     if (dmConv) {
182       ierr = PetscObjectSetOptionsPrefix((PetscObject) dmConv, "conv_par_2_");CHKERRQ(ierr);
183       ierr = DMSetFromOptions(dmConv);CHKERRQ(ierr);
184       ierr = DMDestroy(dm);CHKERRQ(ierr);
185       *dm  = dmConv;
186     }
187     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, "conv_par_2_");CHKERRQ(ierr);
188     ierr = DMViewFromOptions(*dm, NULL, "-dm_view");CHKERRQ(ierr);
189     ierr = PetscObjectSetOptionsPrefix((PetscObject) *dm, NULL);CHKERRQ(ierr);
190 #else
191     SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_SUP,"Recompile with --download-p4est");
192 #endif
193   }
194 
195   /* test redistribution of an already distributed mesh */
196   if (user->redistribute) {
197     DM       distributedMesh;
198     PetscSF  sf;
199     PetscInt nranks;
200 
201     ierr = DMViewFromOptions(*dm, NULL, "-dm_pre_redist_view");CHKERRQ(ierr);
202     ierr = DMPlexDistribute(*dm, 0, NULL, &distributedMesh);CHKERRQ(ierr);
203     if (distributedMesh) {
204       ierr = DMGetPointSF(distributedMesh, &sf);CHKERRQ(ierr);
205       ierr = PetscSFSetUp(sf);CHKERRQ(ierr);
206       ierr = DMGetNeighbors(distributedMesh, &nranks, NULL);CHKERRQ(ierr);
207       ierr = MPI_Allreduce(MPI_IN_PLACE, &nranks, 1, MPIU_INT, MPI_MIN, PetscObjectComm((PetscObject)*dm));CHKERRMPI(ierr);
208       ierr = PetscViewerASCIIPrintf(PETSC_VIEWER_STDOUT_(PetscObjectComm((PetscObject)*dm)), "Minimum number of neighbors: %D\n", nranks);CHKERRQ(ierr);
209       ierr = DMDestroy(dm);CHKERRQ(ierr);
210       *dm  = distributedMesh;
211     }
212     ierr = DMViewFromOptions(*dm, NULL, "-dm_post_redist_view");CHKERRQ(ierr);
213   }
214 
215   if (user->overlap) {
216     DM overlapMesh = NULL;
217 
218     /* Add the overlap to refined mesh */
219     ierr = PetscLogStagePush(user->stages[STAGE_OVERLAP]);CHKERRQ(ierr);
220     ierr = DMViewFromOptions(*dm, NULL, "-dm_pre_overlap_view");CHKERRQ(ierr);
221     ierr = DMPlexDistributeOverlap(*dm, user->overlap, NULL, &overlapMesh);CHKERRQ(ierr);
222     if (overlapMesh) {
223       PetscInt overlap;
224       ierr = DMPlexGetOverlap(overlapMesh, &overlap);CHKERRQ(ierr);
225       ierr = PetscViewerASCIIPrintf(PETSC_VIEWER_STDOUT_WORLD, "Overlap: %D\n", overlap);CHKERRQ(ierr);
226       ierr = DMDestroy(dm);CHKERRQ(ierr);
227       *dm = overlapMesh;
228     }
229     ierr = DMViewFromOptions(*dm, NULL, "-dm_post_overlap_view");CHKERRQ(ierr);
230     ierr = PetscLogStagePop();CHKERRQ(ierr);
231   }
232   if (user->final_ref) {
233     DM refinedMesh = NULL;
234 
235     ierr = DMPlexSetRefinementUniform(*dm, PETSC_TRUE);CHKERRQ(ierr);
236     ierr = DMRefine(*dm, comm, &refinedMesh);CHKERRQ(ierr);
237     if (refinedMesh) {
238       ierr = DMDestroy(dm);CHKERRQ(ierr);
239       *dm  = refinedMesh;
240     }
241   }
242 
243   ierr = PetscObjectSetName((PetscObject) *dm, "Generated Mesh");CHKERRQ(ierr);
244   ierr = DMViewFromOptions(*dm, NULL, "-dm_view");CHKERRQ(ierr);
245   if (user->final_diagnostics) {
246     DMPlexInterpolatedFlag interpolated;
247     PetscInt  dim, depth;
248 
249     ierr = DMGetDimension(*dm, &dim);CHKERRQ(ierr);
250     ierr = DMPlexGetDepth(*dm, &depth);CHKERRQ(ierr);
251     ierr = DMPlexIsInterpolatedCollective(*dm, &interpolated);CHKERRQ(ierr);
252 
253     ierr = DMPlexCheckSymmetry(*dm);CHKERRQ(ierr);
254     if (interpolated == DMPLEX_INTERPOLATED_FULL) {
255       ierr = DMPlexCheckFaces(*dm, 0);CHKERRQ(ierr);
256     }
257     ierr = DMPlexCheckSkeleton(*dm, 0);CHKERRQ(ierr);
258     ierr = DMPlexCheckGeometry(*dm);CHKERRQ(ierr);
259   }
260   ierr = PetscLogEventEnd(user->createMeshEvent,0,0,0,0);CHKERRQ(ierr);
261   PetscFunctionReturn(0);
262 }
263 
264 int main(int argc, char **argv)
265 {
266   DM             dm;
267   AppCtx         user;
268   PetscErrorCode ierr;
269 
270   ierr = PetscInitialize(&argc, &argv, NULL, help);if (ierr) return ierr;
271   ierr = ProcessOptions(PETSC_COMM_WORLD, &user);CHKERRQ(ierr);
272   ierr = CreateMesh(PETSC_COMM_WORLD, &user, &dm);CHKERRQ(ierr);
273   ierr = DMDestroy(&dm);CHKERRQ(ierr);
274   ierr = PetscFinalize();
275   return ierr;
276 }
277 
278 /*TEST
279 
280   # CTetGen 0-1
281   test:
282     suffix: 0
283     requires: ctetgen
284     args: -dm_coord_space 0 -dm_plex_dim 3 -dim 3 -dm_plex_interpolate 0 -ctetgen_verbose 4 -dm_view ascii::ascii_info_detail -info :~sys
285   test:
286     suffix: 1
287     requires: ctetgen
288     args: -dm_coord_space 0 -dm_plex_dim 3 -dim 3 -dm_plex_interpolate 0 -ctetgen_verbose 4 -dm_refine_volume_limit_pre 0.0625 -dm_view ascii::ascii_info_detail -info :~sys
289 
290   # 2D LaTex and ASCII output 2-9
291   test:
292     suffix: 2
293     requires: triangle
294     args: -dm_plex_interpolate 0 -dm_view ascii::ascii_latex
295   test:
296     suffix: 3
297     requires: triangle
298     args: -ref_dm_refine 1 -dm_view ascii::ascii_info_detail
299   test:
300     suffix: 4
301     requires: triangle
302     nsize: 2
303     args: -dm_coord_space 0 -ref_dm_refine 1 -dist_dm_distribute -petscpartitioner_type simple -dm_view ascii::ascii_info_detail
304   test:
305     suffix: 5
306     requires: triangle
307     nsize: 2
308     args: -dm_coord_space 0 -ref_dm_refine 1 -dist_dm_distribute -petscpartitioner_type simple -dm_view ascii::ascii_latex
309   test:
310     suffix: 6
311     args: -dm_coord_space 0 -dm_plex_simplex 0 -dm_view ascii::ascii_info_detail
312   test:
313     suffix: 7
314     args: -dm_coord_space 0 -dm_plex_simplex 0 -ref_dm_refine 1 -dm_view ascii::ascii_info_detail
315   test:
316     suffix: 8
317     nsize: 2
318     args: -dm_plex_simplex 0 -ref_dm_refine 1 -dist_dm_distribute -petscpartitioner_type simple -dm_view ascii::ascii_latex
319 
320   # 1D ASCII output
321   testset:
322     args: -dm_coord_space 0 -dm_plex_dim 1 -dm_view ascii::ascii_info_detail -dm_plex_check_all
323     test:
324       suffix: 1d_0
325       args:
326     test:
327       suffix: 1d_1
328       args: -ref_dm_refine 2
329     test:
330       suffix: 1d_2
331       args: -dm_plex_box_faces 5 -dm_plex_box_bd periodic
332 
333   # Parallel refinement tests with overlap
334   test:
335     suffix: refine_overlap_1d
336     nsize: 2
337     args: -dm_plex_dim 1 -dim 1 -dm_plex_box_faces 4 -dm_plex_box_faces 4 -ref_dm_refine 1 -overlap {{0 1 2}separate output} -dist_dm_distribute -petscpartitioner_type simple -dm_view ascii::ascii_info
338   test:
339     suffix: refine_overlap_2d
340     requires: triangle
341     nsize: {{2 8}separate output}
342     args: -dm_coord_space 0 -ref_dm_refine 1 -dist_dm_distribute -petscpartitioner_type simple -overlap {{0 1 2}separate output} -dm_view ascii::ascii_info
343 
344   # Parallel extrusion tests
345   test:
346     suffix: spheresurface_extruded
347     nsize : 4
348     args: -dm_coord_space 0 -dm_plex_shape sphere -dm_extrude 3 -dist_dm_distribute -petscpartitioner_type simple \
349           -dm_plex_check_all -dm_view ::ascii_info_detail -dm_plex_view_coord_system spherical
350 
351   test:
352     suffix: spheresurface_extruded_symmetric
353     nsize : 4
354     args: -dm_coord_space 0 -dm_plex_shape sphere -dm_extrude 3 -dm_plex_transform_extrude_symmetric -dist_dm_distribute -petscpartitioner_type simple \
355           -dm_plex_check_all -dm_view ::ascii_info_detail -dm_plex_view_coord_system spherical
356 
357   # Parallel simple partitioner tests
358   test:
359     suffix: part_simple_0
360     requires: triangle
361     nsize: 2
362     args: -dm_coord_space 0 -dm_plex_interpolate 0 -dist_dm_distribute -petscpartitioner_type simple -dist_partition_view -dm_view ascii::ascii_info_detail
363   test:
364     suffix: part_simple_1
365     requires: triangle
366     nsize: 8
367     args: -dm_coord_space 0 -ref_dm_refine 1 -dist_dm_distribute -petscpartitioner_type simple -dist_partition_view -dm_view ascii::ascii_info_detail
368 
369   # Parallel partitioner tests
370   test:
371     suffix: part_parmetis_0
372     requires: parmetis
373     nsize: 2
374     args: -dm_plex_simplex 0 -ref_dm_refine 1 -dist_dm_distribute -petscpartitioner_type parmetis -dm_view -petscpartitioner_view -test_redistribute -dm_plex_csr_via_mat {{0 1}} -dm_pre_redist_view ::load_balance -dm_post_redist_view ::load_balance -petscpartitioner_view_graph
375   test:
376     suffix: part_ptscotch_0
377     requires: ptscotch
378     nsize: 2
379     args: -dm_plex_simplex 0 -dist_dm_distribute -petscpartitioner_type ptscotch -petscpartitioner_view -petscpartitioner_ptscotch_strategy quality -test_redistribute -dm_plex_csr_via_mat {{0 1}} -dm_pre_redist_view ::load_balance -dm_post_redist_view ::load_balance -petscpartitioner_view_graph
380   test:
381     suffix: part_ptscotch_1
382     requires: ptscotch
383     nsize: 8
384     args: -dm_plex_simplex 0 -ref_dm_refine 1 -dist_dm_distribute -petscpartitioner_type ptscotch -petscpartitioner_view -petscpartitioner_ptscotch_imbalance 0.1
385 
386   # CGNS reader tests 10-11 (need to find smaller test meshes)
387   test:
388     suffix: cgns_0
389     requires: cgns
390     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/tut21.cgns -dm_view
391 
392   # Gmsh mesh reader tests
393   testset:
394     args: -dm_coord_space 0 -dm_view
395 
396     test:
397       suffix: gmsh_0
398       requires: !single
399       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/doublet-tet.msh
400     test:
401       suffix: gmsh_1
402       requires: !single
403       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square.msh
404     test:
405       suffix: gmsh_2
406       requires: !single
407       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_bin.msh
408     test:
409       suffix: gmsh_3
410       nsize: 3
411       requires: !single
412       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square.msh -dist_dm_distribute -petscpartitioner_type simple
413     test:
414       suffix: gmsh_4
415       nsize: 3
416       requires: !single
417       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_bin.msh -dist_dm_distribute -petscpartitioner_type simple
418     test:
419       suffix: gmsh_5
420       requires: !single
421       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_quad.msh
422     # TODO: it seems the mesh is not a valid gmsh (inverted cell)
423     test:
424       suffix: gmsh_6
425       requires: !single
426       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_bin_physnames.msh -final_diagnostics 0
427     test:
428       suffix: gmsh_7
429       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/mesh-3d-box-innersphere_bin.msh -dm_view ::ascii_info_detail -dm_plex_check_all
430     test:
431       suffix: gmsh_8
432       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/mesh-3d-box-innersphere.msh -dm_view ::ascii_info_detail -dm_plex_check_all
433   testset:
434     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic_bin.msh -dm_view ::ascii_info_detail -dm_plex_check_all
435     test:
436       suffix: gmsh_9
437     test:
438       suffix: gmsh_9_periodic_0
439       args: -dm_plex_gmsh_periodic 0
440   testset:
441     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_view ::ascii_info_detail -dm_plex_check_all
442     test:
443       suffix: gmsh_10
444     test:
445       suffix: gmsh_10_periodic_0
446       args: -dm_plex_gmsh_periodic 0
447   testset:
448     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_view ::ascii_info_detail -dm_plex_check_all -ref_dm_refine 1
449     test:
450       suffix: gmsh_11
451     test:
452       suffix: gmsh_11_periodic_0
453       args: -dm_plex_gmsh_periodic 0
454   # TODO: it seems the mesh is not a valid gmsh (inverted cell)
455   test:
456     suffix: gmsh_12
457     nsize: 4
458     requires: !single mpiio
459     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_bin_physnames.msh -viewer_binary_mpiio -dist_dm_distribute -petscpartitioner_type simple -dm_view -final_diagnostics 0
460   test:
461     suffix: gmsh_13_hybs2t
462     nsize: 4
463     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_triquad.msh -dist_dm_distribute -petscpartitioner_type simple -dm_view -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox -dm_plex_check_all
464   test:
465     suffix: gmsh_14_ext
466     requires: !single
467     args: -dm_coord_space 0 -dm_extrude 2 -dm_plex_transform_extrude_thickness 1.5 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_bin.msh -dm_view -dm_plex_check_all
468   test:
469     suffix: gmsh_14_ext_s2t
470     requires: !single
471     args: -dm_coord_space 0 -dm_extrude 2 -dm_plex_transform_extrude_thickness 1.5 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_bin.msh -dm_view -dm_plex_check_all -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox
472   test:
473     suffix: gmsh_15_hyb3d
474     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_tetwedge.msh -dm_view -dm_plex_check_all
475   test:
476     suffix: gmsh_15_hyb3d_vtk
477     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_tetwedge.msh -dm_view vtk: -dm_plex_gmsh_hybrid -dm_plex_check_all
478   test:
479     suffix: gmsh_15_hyb3d_s2t
480     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_tetwedge.msh -dm_view -dm_plex_check_all -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox
481   test:
482     suffix: gmsh_16_spheresurface
483     nsize : 4
484     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3 -dm_plex_check_all -dm_view -dist_dm_distribute -petscpartitioner_type simple
485   test:
486     suffix: gmsh_16_spheresurface_s2t
487     nsize : 4
488     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3 -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox -dm_plex_check_all -dm_view -dist_dm_distribute -petscpartitioner_type simple
489   test:
490     suffix: gmsh_16_spheresurface_extruded
491     nsize : 4
492     args: -dm_coord_space 0 -dm_extrude 3 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3 -dm_plex_check_all -dm_view -dist_dm_distribute -petscpartitioner_type simple
493   test:
494     suffix: gmsh_16_spheresurface_extruded_s2t
495     nsize : 4
496     args: -dm_coord_space 0 -dm_extrude 3 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3 -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox -dm_plex_check_all -dm_view -dist_dm_distribute -petscpartitioner_type simple
497   test:
498     suffix: gmsh_17_hyb3d_interp_ascii
499     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_hexwedge.msh -dm_view -dm_plex_check_all
500   test:
501     suffix: exodus_17_hyb3d_interp_ascii
502     requires: exodusii
503     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_hexwedge.exo -dm_view -dm_plex_check_all
504 
505   # Legacy Gmsh v22/v40 ascii/binary reader tests
506   testset:
507     output_file: output/ex1_gmsh_3d_legacy.out
508     args: -dm_coord_space 0 -dm_view ::ascii_info_detail -dm_plex_check_all
509     test:
510       suffix: gmsh_3d_ascii_v22
511       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-ascii.msh2
512     test:
513       suffix: gmsh_3d_ascii_v40
514       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-ascii.msh4
515     test:
516       suffix: gmsh_3d_binary_v22
517       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary.msh2
518     test:
519       suffix: gmsh_3d_binary_v40
520       requires: long64
521       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary.msh4
522 
523   # Gmsh v41 ascii/binary reader tests
524   testset: # 32bit mesh, sequential
525     args: -dm_coord_space 0 -dm_view ::ascii_info_detail -dm_plex_check_all
526     output_file: output/ex1_gmsh_3d_32.out
527     test:
528       suffix: gmsh_3d_ascii_v41_32
529       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-ascii-32.msh
530     test:
531       suffix: gmsh_3d_binary_v41_32
532       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary-32.msh
533     test:
534       suffix: gmsh_3d_binary_v41_32_mpiio
535       requires: defined(PETSC_HAVE_MPIIO)
536       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary-32.msh -viewer_binary_mpiio
537   testset:  # 32bit mesh, parallel
538     args: -dm_coord_space 0 -dist_dm_distribute -petscpartitioner_type simple -dm_view ::ascii_info_detail -dm_plex_check_all
539     nsize: 2
540     output_file: output/ex1_gmsh_3d_32_np2.out
541     test:
542       suffix: gmsh_3d_ascii_v41_32_np2
543       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-ascii-32.msh
544     test:
545       suffix: gmsh_3d_binary_v41_32_np2
546       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary-32.msh
547     test:
548       suffix: gmsh_3d_binary_v41_32_np2_mpiio
549       requires: defined(PETSC_HAVE_MPIIO)
550       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary-32.msh -viewer_binary_mpiio
551   testset: # 64bit mesh, sequential
552     args: -dm_coord_space 0 -dm_view ::ascii_info_detail -dm_plex_check_all
553     output_file: output/ex1_gmsh_3d_64.out
554     test:
555       suffix: gmsh_3d_ascii_v41_64
556       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-ascii-64.msh
557     test:
558       suffix: gmsh_3d_binary_v41_64
559       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary-64.msh
560     test:
561       suffix: gmsh_3d_binary_v41_64_mpiio
562       requires: defined(PETSC_HAVE_MPIIO)
563       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary-64.msh -viewer_binary_mpiio
564   testset:  # 64bit mesh, parallel
565     args: -dm_coord_space 0 -dist_dm_distribute -petscpartitioner_type simple -dm_view ::ascii_info_detail -dm_plex_check_all
566     nsize: 2
567     output_file: output/ex1_gmsh_3d_64_np2.out
568     test:
569       suffix: gmsh_3d_ascii_v41_64_np2
570       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-ascii-64.msh
571     test:
572       suffix: gmsh_3d_binary_v41_64_np2
573       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary-64.msh
574     test:
575       suffix: gmsh_3d_binary_v41_64_np2_mpiio
576       requires: defined(PETSC_HAVE_MPIIO)
577       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/gmsh-3d-binary-64.msh -viewer_binary_mpiio
578 
579   # Fluent mesh reader tests
580   # TODO: Geometry checks fail
581   test:
582     suffix: fluent_0
583     requires: !complex
584     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square.cas -dm_view -final_diagnostics 0
585   test:
586     suffix: fluent_1
587     nsize: 3
588     requires: !complex
589     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square.cas -dist_dm_distribute -petscpartitioner_type simple -dm_view -final_diagnostics 0
590   test:
591     suffix: fluent_2
592     requires: !complex
593     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/cube_5tets_ascii.cas -dm_view -final_diagnostics 0
594   test:
595     suffix: fluent_3
596     requires: !complex
597     TODO: Fails on non-linux: fseek(), fileno() ? https://gitlab.com/petsc/petsc/merge_requests/2206#note_238166382
598     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/cube_5tets.cas -dm_view -final_diagnostics 0
599 
600   # Med mesh reader tests, including parallel file reads
601   test:
602     suffix: med_0
603     requires: med
604     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square.med -dm_view
605   test:
606     suffix: med_1
607     requires: med
608     nsize: 3
609     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square.med -dist_dm_distribute -petscpartitioner_type simple -dm_view
610   test:
611     suffix: med_2
612     requires: med
613     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/cylinder.med -dm_view
614   test:
615     suffix: med_3
616     requires: med
617     TODO: MED
618     nsize: 3
619     args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/cylinder.med -dist_dm_distribute -petscpartitioner_type simple -dm_view
620 
621   # Test shape quality
622   test:
623     suffix: test_shape
624     requires: ctetgen
625     args: -dm_plex_dim 3 -dim 3 -dm_refine_hierarchy 3 -dm_plex_check_all -dm_plex_check_cell_shape
626 
627   # Test simplex to tensor conversion
628   test:
629     suffix: s2t2
630     requires: triangle
631     args: -dm_coord_space 0 -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox -dm_refine_volume_limit_pre 0.0625 -dm_view ascii::ascii_info_detail
632 
633   test:
634     suffix: s2t3
635     requires: ctetgen
636     args: -dm_coord_space 0 -dm_plex_dim 3 -dim 3 -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox -dm_refine_volume_limit_pre 0.0625 -dm_view ascii::ascii_info_detail
637 
638   # Test cylinder
639   testset:
640     args: -dm_plex_shape cylinder -dm_plex_check_all -dm_view
641     test:
642       suffix: cylinder
643       args: -ref_dm_refine 1
644     test:
645       suffix: cylinder_per
646       args: -dm_plex_cylinder_bd periodic -ref_dm_refine 1 -ref_dm_refine_remap 0
647     test:
648       suffix: cylinder_wedge
649       args: -dm_coord_space 0 -dm_plex_interpolate 0 -dm_plex_cell tensor_triangular_prism -dm_view vtk:
650     test:
651       suffix: cylinder_wedge_int
652       output_file: output/ex1_cylinder_wedge.out
653       args: -dm_coord_space 0 -dm_plex_cell tensor_triangular_prism -dm_view vtk:
654 
655   test:
656     suffix: box_2d
657     args: -dm_plex_simplex 0 -ref_dm_refine 2 -dm_plex_check_all -dm_view
658 
659   test:
660     suffix: box_2d_per
661     args: -dm_plex_simplex 0 -ref_dm_refine 2 -dm_plex_check_all -dm_view
662 
663   test:
664     suffix: box_2d_per_unint
665     args: -dm_coord_space 0 -dm_plex_simplex 0 -dm_plex_interpolate 0 -dm_plex_box_faces 3,3 -dm_plex_box_faces 3,3 -dm_plex_check_all -dm_view ::ascii_info_detail
666 
667   test:
668     suffix: box_3d
669     args: -dm_plex_dim 3 -dim 3 -dm_plex_simplex 0 -ref_dm_refine 3 -dm_plex_check_all -dm_view
670 
671   test:
672     requires: triangle
673     suffix: box_wedge
674     args: -dm_coord_space 0 -dm_plex_dim 3 -dim 3 -dm_plex_simplex 0 -dm_plex_cell tensor_triangular_prism -dm_view vtk: -dm_plex_check_all
675 
676   testset:
677     requires: triangle
678     args: -dm_coord_space 0 -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_cell tensor_triangular_prism -dm_plex_box_faces 2,3,1 -dm_view -dm_plex_check_all -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox
679     test:
680       suffix: box_wedge_s2t
681     test:
682       nsize: 3
683       args: -dist_dm_distribute -petscpartitioner_type simple
684       suffix: box_wedge_s2t_parallel
685 
686   # Test GLVis output
687   testset:
688     args: -dm_coord_space 0 -dm_plex_interpolate 0
689     test:
690       suffix: glvis_2d_tet
691       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_plex_gmsh_periodic 0 -dm_view glvis:
692     test:
693       suffix: glvis_2d_tet_per
694       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary 0
695     test:
696       suffix: glvis_3d_tet
697       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/mesh-3d-box-innersphere_bin.msh -dm_plex_gmsh_periodic 0 -dm_view glvis:
698   testset:
699     args: -dm_coord_space 0
700     test:
701       suffix: glvis_2d_tet_per_mfem
702       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -viewer_glvis_dm_plex_enable_boundary -viewer_glvis_dm_plex_enable_mfem -dm_view glvis:
703     test:
704       suffix: glvis_2d_quad
705       args: -dm_plex_simplex 0 -dm_plex_box_faces 3,3 -dm_view glvis:
706     test:
707       suffix: glvis_2d_quad_per
708       args: -dm_plex_simplex 0 -dm_plex_box_faces 3,3 -dm_plex_box_bd periodic,periodic -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary
709     test:
710       suffix: glvis_2d_quad_per_mfem
711       args: -dm_plex_simplex 0 -dm_plex_box_faces 3,3 -dm_plex_box_bd periodic,periodic -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary -viewer_glvis_dm_plex_enable_mfem
712     test:
713       suffix: glvis_3d_tet_per
714       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/mesh-3d-box-innersphere_bin.msh -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary
715     test:
716       suffix: glvis_3d_tet_per_mfem
717       TODO: broken
718       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/mesh-3d-box-innersphere_bin.msh -viewer_glvis_dm_plex_enable_mfem -dm_view glvis:
719     test:
720       suffix: glvis_3d_hex
721       args: -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_box_faces 3,3,3 -dm_view glvis:
722     test:
723       suffix: glvis_3d_hex_per
724       args: -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_box_faces 3,3,3 -dm_plex_box_bd periodic,periodic,periodic -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary 0
725     test:
726       suffix: glvis_3d_hex_per_mfem
727       args: -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_box_faces 3,3,3 -dm_plex_box_bd periodic,periodic,periodic -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary -viewer_glvis_dm_plex_enable_mfem
728     test:
729       suffix: glvis_2d_hyb
730       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_triquad.msh -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary
731     test:
732       suffix: glvis_3d_hyb
733       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_tetwedge.msh -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary
734     test:
735       suffix: glvis_3d_hyb_s2t
736       args: -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_3d_cube.msh -dm_view glvis: -viewer_glvis_dm_plex_enable_boundary -ref_dm_refine 1 -ref_dm_plex_transform_type refine_tobox -dm_plex_check_all
737 
738   # Test P4EST
739   testset:
740     requires: p4est
741     args: -dm_coord_space 0 -dm_view -test_p4est_seq -conv_seq_2_dm_plex_check_all -conv_seq_1_dm_forest_minimum_refinement 1
742     test:
743       suffix: p4est_periodic
744       args: -dm_plex_simplex 0 -dm_plex_box_bd periodic,periodic -dm_plex_box_faces 3,5 -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 2 -conv_seq_1_dm_p4est_refine_pattern hash
745     test:
746       suffix: p4est_periodic_3d
747       args: -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_box_bd periodic,periodic,none -dm_plex_box_faces 3,5,4 -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 2 -conv_seq_1_dm_p4est_refine_pattern hash
748     test:
749       suffix: p4est_gmsh_periodic
750       args: -dm_coord_space 0 -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 1 -conv_seq_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh
751     test:
752       suffix: p4est_gmsh_surface
753       args: -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 1 -conv_seq_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3
754     test:
755       suffix: p4est_gmsh_surface_parallel
756       nsize: 2
757       args: -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 1 -conv_seq_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3 -petscpartitioner_type simple -dm_view ::load_balance
758     test:
759       suffix: p4est_hyb_2d
760       args: -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 1 -conv_seq_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_triquad.msh
761     test:
762       suffix: p4est_hyb_3d
763       args: -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 1 -conv_seq_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_tetwedge.msh
764     test:
765       requires: ctetgen
766       suffix: p4est_s2t_bugfaces_3d
767       args: -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 0 -dm_plex_dim 3 -dm_plex_box_faces 1,1
768     test:
769       suffix: p4est_bug_overlapsf
770       nsize: 3
771       args: -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_box_faces 2,2,1 -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 1 -conv_seq_1_dm_p4est_refine_pattern hash -petscpartitioner_type simple
772     test:
773       suffix: p4est_redistribute
774       nsize: 3
775       args: -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_box_faces 2,2,1 -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 1 -conv_seq_1_dm_p4est_refine_pattern hash -petscpartitioner_type simple -test_redistribute -dm_plex_csr_via_mat {{0 1}} -dm_view ::load_balance
776     test:
777       suffix: p4est_gmsh_s2t_3d
778       args: -conv_seq_1_dm_forest_initial_refinement 1 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/doublet-tet.msh
779     test:
780       suffix: p4est_gmsh_s2t_3d_hash
781       args: -conv_seq_1_dm_forest_initial_refinement 1 -conv_seq_1_dm_forest_maximum_refinement 2 -conv_seq_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/doublet-tet.msh
782     test:
783       requires: long_runtime
784       suffix: p4est_gmsh_periodic_3d
785       args: -dm_coord_space 0 -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 1 -conv_seq_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/mesh-3d-box-innersphere.msh
786 
787   testset:
788     requires: p4est
789     nsize: 6
790     args: -dm_coord_space 0 -test_p4est_par -conv_par_2_dm_plex_check_all -conv_par_1_dm_forest_minimum_refinement 1 -conv_par_1_dm_forest_partition_overlap 0 -dist_dm_distribute
791     test:
792       TODO: interface cones do not conform
793       suffix: p4est_par_periodic
794       args: -dm_plex_simplex 0 -dm_plex_box_bd periodic,periodic -dm_plex_box_faces 3,5 -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 2 -conv_par_1_dm_p4est_refine_pattern hash
795     test:
796       TODO: interface cones do not conform
797       suffix: p4est_par_periodic_3d
798       args: -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_box_bd periodic,periodic,periodic -dm_plex_box_faces 3,5,4 -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 2 -conv_par_1_dm_p4est_refine_pattern hash
799     test:
800       TODO: interface cones do not conform
801       suffix: p4est_par_gmsh_periodic
802       args: -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 1 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh
803     test:
804       suffix: p4est_par_gmsh_surface
805       args: -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 1 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3
806     test:
807       suffix: p4est_par_gmsh_s2t_3d
808       args: -conv_par_1_dm_forest_initial_refinement 1 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/doublet-tet.msh
809     test:
810       TODO: interface cones do not conform
811       suffix: p4est_par_gmsh_s2t_3d_hash
812       args: -conv_par_1_dm_forest_initial_refinement 1 -conv_par_1_dm_forest_maximum_refinement 2 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/doublet-tet.msh
813     test:
814       requires: long_runtime
815       suffix: p4est_par_gmsh_periodic_3d
816       args: -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 1 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/mesh-3d-box-innersphere.msh
817 
818   testset:
819     requires: p4est
820     nsize: 6
821     args: -dm_coord_space 0 -test_p4est_par -conv_par_2_dm_plex_check_all -conv_par_1_dm_forest_minimum_refinement 1 -conv_par_1_dm_forest_partition_overlap 1 -dist_dm_distribute -petscpartitioner_type simple
822     test:
823       suffix: p4est_par_ovl_periodic
824       args: -dm_plex_simplex 0 -dm_plex_box_bd periodic,periodic -dm_plex_box_faces 3,5 -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 2 -conv_par_1_dm_p4est_refine_pattern hash
825     #TODO Mesh cell 201 is inverted, vol = 0. (FVM Volume. Is it correct? -> Diagnostics disabled)
826     test:
827       suffix: p4est_par_ovl_periodic_3d
828       args: -dm_plex_dim 3 -dm_plex_simplex 0 -dm_plex_box_bd periodic,periodic,none -dm_plex_box_faces 3,5,4 -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 2 -conv_par_1_dm_p4est_refine_pattern hash -final_diagnostics 0
829     test:
830       suffix: p4est_par_ovl_gmsh_periodic
831       args: -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 1 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh
832     test:
833       suffix: p4est_par_ovl_gmsh_surface
834       args: -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 1 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3
835     test:
836       suffix: p4est_par_ovl_gmsh_s2t_3d
837       args: -conv_par_1_dm_forest_initial_refinement 1 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/doublet-tet.msh
838     test:
839       suffix: p4est_par_ovl_gmsh_s2t_3d_hash
840       args: -conv_par_1_dm_forest_initial_refinement 1 -conv_par_1_dm_forest_maximum_refinement 2 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/doublet-tet.msh
841     test:
842       requires: long_runtime
843       suffix: p4est_par_ovl_gmsh_periodic_3d
844       args: -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 1 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/mesh-3d-box-innersphere.msh
845     test:
846       suffix: p4est_par_ovl_hyb_2d
847       args: -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 1 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_triquad.msh
848     test:
849       suffix: p4est_par_ovl_hyb_3d
850       args: -conv_par_1_dm_forest_initial_refinement 0 -conv_par_1_dm_forest_maximum_refinement 1 -conv_par_1_dm_p4est_refine_pattern hash -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_tetwedge.msh
851 
852   test:
853     TODO: broken
854     requires: p4est
855     nsize: 2
856     suffix: p4est_bug_labels_noovl
857     args: -test_p4est_seq -dm_plex_check_all -dm_forest_minimum_refinement 0 -dm_forest_partition_overlap 1 -dm_plex_simplex 0 -dm_plex_box_faces 3,3 -dm_forest_initial_refinement 0 -dm_forest_maximum_refinement 2 -dm_p4est_refine_pattern hash -dist_dm_distribute -petscpartitioner_type simple -dm_forest_print_label_error
858 
859   test:
860     requires: p4est
861     nsize: 2
862     suffix: p4est_bug_distribute_overlap
863     args: -dm_coord_space 0 -test_p4est_seq -conv_seq_2_dm_plex_check_all -conv_seq_1_dm_forest_minimum_refinement 0 -conv_seq_1_dm_forest_partition_overlap 0 -dm_plex_simplex 0 -dm_plex_box_faces 3,3 -conv_seq_1_dm_forest_initial_refinement 0 -conv_seq_1_dm_forest_maximum_refinement 2 -conv_seq_1_dm_p4est_refine_pattern hash -petscpartitioner_type simple -overlap 1 -dm_view ::load_balance
864     args: -dm_post_overlap_view
865 
866   test:
867     suffix: ref_alfeld2d_0
868     requires: triangle
869     args: -dm_plex_box_faces 5,3 -dm_view -dm_plex_check_all -ref_dm_refine 1 -ref_dm_plex_transform_type refine_alfeld -final_diagnostics
870   test:
871     suffix: ref_alfeld3d_0
872     requires: ctetgen
873     args: -dm_plex_dim 3 -dm_plex_box_faces 5,1,1 -dm_view -dm_plex_check_all -ref_dm_refine 1 -ref_dm_plex_transform_type refine_alfeld -final_diagnostics
874 
875   # Boundary layer refiners
876   test:
877     suffix: ref_bl_1
878     args: -dm_plex_dim 1 -dm_plex_simplex 0 -dm_plex_box_faces 5,1 -dm_view -dm_plex_check_all 0 -ref_dm_refine 1 -ref_dm_plex_transform_type refine_boundary_layer -dm_extrude 2 -final_diagnostics -ref_dm_plex_transform_bl_splits 3
879   test:
880     suffix: ref_bl_2_tri
881     requires: triangle
882     args: -dm_plex_box_faces 5,3 -dm_view -dm_plex_check_all 0 -ref_dm_refine 1 -ref_dm_plex_transform_type refine_boundary_layer -dm_extrude 3 -final_diagnostics -ref_dm_plex_transform_bl_splits 4
883   test:
884     suffix: ref_bl_3_quad
885     args: -dm_plex_simplex 0 -dm_plex_box_faces 5,1 -dm_view -dm_plex_check_all 0 -ref_dm_refine 1 -ref_dm_plex_transform_type refine_boundary_layer -dm_extrude 3 -final_diagnostics -ref_dm_plex_transform_bl_splits 4
886   test:
887     suffix: ref_bl_spheresurface_extruded
888     nsize : 4
889     args: -dm_extrude 3 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/surfacesphere_bin.msh -dm_plex_gmsh_spacedim 3 -dm_plex_check_all -dm_view -dist_dm_distribute -petscpartitioner_type simple -final_diagnostics -ref_dm_refine 1 -ref_dm_plex_transform_type refine_boundary_layer -ref_dm_plex_transform_bl_splits 2
890   test:
891     suffix: ref_bl_3d_hyb
892     nsize : 4
893     args: -dm_coord_space 0 -dm_plex_filename ${wPETSC_DIR}/share/petsc/datafiles/meshes/hybrid_3d_cube.msh -dm_plex_check_all -dm_view -dist_dm_distribute -petscpartitioner_type simple -final_diagnostics -ref_dm_refine 1 -ref_dm_plex_transform_type refine_boundary_layer -ref_dm_plex_transform_bl_splits 4 -ref_dm_plex_transform_bl_height_factor 3.1
894 
895   testset:
896     args: -dm_plex_shape sphere -dm_plex_check_all -dm_view
897     test:
898       suffix: sphere_0
899       args:
900     test:
901       suffix: sphere_1
902       args: -ref_dm_refine 2
903     test:
904       suffix: sphere_2
905       args: -dm_plex_simplex 0
906     test:
907       suffix: sphere_3
908       args: -dm_plex_simplex 0 -ref_dm_refine 2
909 
910   test:
911     suffix: ball_0
912     requires: ctetgen
913     args: -dm_plex_dim 3 -dm_plex_shape ball -dm_plex_check_all -dm_view
914 
915   test:
916     suffix: ball_1
917     requires: ctetgen
918     args: -dm_plex_dim 3 -dm_plex_shape ball -bd_dm_refine 2 -dm_plex_check_all -dm_view
919 
920 TEST*/
921