20 …PetscInt               mwidth, nwidth, pwidth, cnt, mp, np, pp, Ntotal, gl[26], *globals, Ng, *IIi…  in DMDAGetWireBasketInterpolation()  local
237   gl[cnt++] = 0;  in DMDAGetWireBasketInterpolation()
239     gl[cnt++] = 1;  in DMDAGetWireBasketInterpolation()
241   gl[cnt++] = m - istart - 1;  in DMDAGetWireBasketInterpolation()
243     gl[cnt++] = mwidth;  in DMDAGetWireBasketInterpolation()
245       gl[cnt++] = mwidth + 1;  in DMDAGetWireBasketInterpolation()
247     gl[cnt++] = mwidth + m - istart - 1;  in DMDAGetWireBasketInterpolation()
249   gl[cnt++] = mwidth * (n - jstart - 1);  in DMDAGetWireBasketInterpolation()
251     gl[cnt++] = mwidth * (n - jstart - 1) + 1;  in DMDAGetWireBasketInterpolation()
253   gl[cnt++] = mwidth * (n - jstart - 1) + m - istart - 1;  in DMDAGetWireBasketInterpolation()
255     gl[cnt++] = mwidth * nwidth;  in DMDAGetWireBasketInterpolation()
257       gl[cnt++] = mwidth * nwidth + 1;  in DMDAGetWireBasketInterpolation()
259     gl[cnt++] = mwidth * nwidth + m - istart - 1;  in DMDAGetWireBasketInterpolation()
261       gl[cnt++] = mwidth * nwidth + mwidth; /* these are the interior nodes */  in DMDAGetWireBasketInterpolation()
262       gl[cnt++] = mwidth * nwidth + mwidth + m - istart - 1;  in DMDAGetWireBasketInterpolation()
264     gl[cnt++] = mwidth * nwidth + mwidth * (n - jstart - 1);  in DMDAGetWireBasketInterpolation()
266       gl[cnt++] = mwidth * nwidth + mwidth * (n - jstart - 1) + 1;  in DMDAGetWireBasketInterpolation()
268     gl[cnt++] = mwidth * nwidth + mwidth * (n - jstart - 1) + m - istart - 1;  in DMDAGetWireBasketInterpolation()
270   gl[cnt++] = mwidth * nwidth * (p - kstart - 1);  in DMDAGetWireBasketInterpolation()
272     gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + 1;  in DMDAGetWireBasketInterpolation()
274   gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + m - istart - 1;  in DMDAGetWireBasketInterpolation()
276     gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + mwidth;  in DMDAGetWireBasketInterpolation()
278       gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + mwidth + 1;  in DMDAGetWireBasketInterpolation()
280     gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + mwidth + m - istart - 1;  in DMDAGetWireBasketInterpolation()
282   gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + mwidth * (n - jstart - 1);  in DMDAGetWireBasketInterpolation()
284     gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + mwidth * (n - jstart - 1) + 1;  in DMDAGetWireBasketInterpolation()
286   gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + mwidth * (n - jstart - 1) + m - istart - 1;  in DMDAGetWireBasketInterpolation()
290   PetscCall(ISLocalToGlobalMappingApply(ltg, 26, gl, gl));  in DMDAGetWireBasketInterpolation()
293 …PetscCallMPI(MPI_Allgather(gl, 26, MPIU_INT, globals, 26, MPIU_INT, PetscObjectComm((PetscObject)d…  in DMDAGetWireBasketInterpolation()
310     PetscCall(PetscHMapIGetWithDefault(ht, gl[i] + 1, 0, gl + i));  in DMDAGetWireBasketInterpolation()
311     --gl[i];  in DMDAGetWireBasketInterpolation()
325   PetscCall(MatSetValues(*P, Nint, IIint, 26, gl, rxint, INSERT_VALUES));  in DMDAGetWireBasketInterpolation()
328   PetscCall(MatSetValues(*P, Nsurf, IIsurf, 26, gl, rxint, INSERT_VALUES));  in DMDAGetWireBasketInterpolation()
368 …PetscInt               mwidth, nwidth, pwidth, cnt, mp, np, pp, Ntotal, gl[6], *globals, Ng, *IIin…  in DMDAGetFaceInterpolation()  local
559     gl[cnt++] = mwidth + 1;  in DMDAGetFaceInterpolation()
563       gl[cnt++] = mwidth * nwidth + 1;  in DMDAGetFaceInterpolation()
566       gl[cnt++] = mwidth * nwidth + mwidth; /* these are the interior nodes */  in DMDAGetFaceInterpolation()
567       gl[cnt++] = mwidth * nwidth + mwidth + m - istart - 1;  in DMDAGetFaceInterpolation()
570       gl[cnt++] = mwidth * nwidth + mwidth * (n - jstart - 1) + 1;  in DMDAGetFaceInterpolation()
574     gl[cnt++] = mwidth * nwidth * (p - kstart - 1) + mwidth + 1;  in DMDAGetFaceInterpolation()
579   PetscCall(ISLocalToGlobalMappingApply(ltg, 6, gl, gl));  in DMDAGetFaceInterpolation()
582 …PetscCallMPI(MPI_Allgather(gl, 6, MPIU_INT, globals, 6, MPIU_INT, PetscObjectComm((PetscObject)da)…  in DMDAGetFaceInterpolation()
599     PetscCall(PetscHMapIGetWithDefault(ht, gl[i] + 1, 0, gl + i));  in DMDAGetFaceInterpolation()
600     --gl[i];  in DMDAGetFaceInterpolation()
614   PetscCall(MatSetValues(*P, Nint, IIint, 6, gl, rxint, INSERT_VALUES));  in DMDAGetFaceInterpolation()
617   PetscCall(MatSetValues(*P, Nsurf, IIsurf, 6, gl, rxint, INSERT_VALUES));  in DMDAGetFaceInterpolation()