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()
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()
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 …d]) < PETSC_SMALL, PETSC_COMM_WORLD, PETSC_ERR_ARG_INCOMP, "Error vertex v %" PetscInt_FMT " (dim …  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