24 …PetscCall(PetscOptionsFList("-dm_type", "DM implementation on which to define field", "ex20.c", DM…  in main()
27   PetscCall(PetscOptionsGetBool(NULL, NULL, "-ex20_print", &print, NULL));  in main()
28   PetscCall(PetscOptionsGetInt(NULL, NULL, "-dm_plex_dim", &dim, NULL));  in main()
31   PetscCall(PetscOptionsGetIntArray(NULL, NULL, "-dm_plex_box_faces", faces, &n, NULL));  in main()
34   PetscCall(PetscOptionsGetRealArray(NULL, NULL, "-dm_plex_box_lower", lower, &n, NULL));  in main()
37   PetscCall(PetscOptionsGetRealArray(NULL, NULL, "-dm_plex_box_upper", upper, &n, NULL));  in main()
59       for (PetscInt d = 0; d < dim; d++) {  in main()  local
66         PetscCall(DMCreateGlobalVector(dm, &Vecxyz[d]));  in main()
67 …PetscCall(PetscObjectSetName((PetscObject)Vecxyz[d], d == 0 ? "x_coordinate" : (d == 1 ? "y_coordi…  in main()
81           xArray[sOff] = coordArray[cOff + d];  in main()
86         PetscCall(DMLocalToGlobal(dm, loc_vec, INSERT_VALUES, Vecxyz[d]));  in main()
88         PetscCall(VecGetSize(Vecxyz[d], &n_state_global));  in main()
92 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Created DMPlex in %" PetscInt_FMT "D with faces (%" Petsc…  in main()
110       for (PetscInt d = 0; d < dim; d++) {  in main()  local
111         PetscCall(DMCreateGlobalVector(dm, &Vecxyz[d]));  in main()
112         PetscCall(VecSetFromOptions(Vecxyz[d]));  in main()
113 …PetscCall(PetscObjectSetName((PetscObject)Vecxyz[d], d == 0 ? "x_coordinate" : (d == 1 ? "y_coordi…  in main()
114         PetscCall(VecStrideGather(coord, d, Vecxyz[d], INSERT_VALUES));  in main()
115         PetscCall(VecGetSize(Vecxyz[d], &n_state_global));  in main()
118 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Created DMDA of type %s in %" PetscInt_FMT "D with faces …  in main()
120   PetscCall(DMViewFromOptions(dm, NULL, "-ex20_dm_view")); // PetscSleep(10);  in main()
124   for (PetscInt d = 0; d < dim && d < 2; d++) nvertexobs *= 3;  in main()  local
137     for (PetscInt d = 0; d < dim; d++) PetscCall(VecGetArrayRead(Vecxyz[d], &coords[d]));  in main()  local
138     for (PetscInt d = 0; d < dim; d++) gridSpacing[d] = (upper[d] - lower[d]) / faces[d];  in main()  local
144       for (PetscInt d = 0; d < dim; d++) c[d] = PetscRealPart(coords[d][v]);  in main()  local
146       for (PetscInt d = 0; d < dim; d++) {  in main()  local
147         PetscReal relCoord = c[d] - lower[d];  in main()
148         PetscInt  gridIdx  = (PetscInt)PetscFloorReal(relCoord / gridSpacing[d] + 0.5);  in main()
149 …- gridIdx * gridSpacing[d]) < PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_ARG_INCOMP, "Error vertex v…  in main()
157     for (PetscInt d = 0; d < dim; d++) PetscCall(VecRestoreArrayRead(Vecxyz[d], &coords[d]));  in main()  local
195     for (PetscInt d = 0; d < dim; d++) PetscCall(VecGetArray(Vecxyz[d], &coords[d]));  in main()  local
198       for (PetscInt d = 0; d < dim; d++) {  in main()  local
199         PetscReal noise, gridSpacing = (upper[d] - lower[d]) / faces[d];  in main()
203         coords[d][v] += (noise - 0.5) * 0.001 * gridSpacing;  in main()
206     for (PetscInt d = 0; d < dim; d++) PetscCall(VecRestoreArray(Vecxyz[d], &coords[d]));  in main()  local
218   for (PetscInt d = 0; d < dim; d++) PetscCall(VecDestroy(&Vecxyz[d]));  in main()  local
231     diff_args: -j
232 …-dm_plex_dim 1 -dm_plex_box_faces 16 -dm_plex_simplex 0 -dm_plex_box_bd periodic -dm_plex_box_uppe…
237     diff_args: -j
238 …-dm_plex_dim 2 -dm_plex_box_faces 7,7 -dm_plex_simplex 0 -dm_plex_box_bd periodic,none -dm_plex_bo…
243     diff_args: -j
244 …args: -dm_type da -dm_plex_dim 2 -dm_plex_box_faces 7,7 -dm_plex_box_upper 5,5 -ex20_print -ex20_d…
249     diff_args: -j
250 …-dm_plex_dim 3 -dm_plex_box_faces 5,5,5 -dm_plex_simplex 0 -dm_plex_box_bd periodic,none,none -dm_…