19 …PetscInt               dim, i, j, k, m, n, p, dof, Nint, Nface, Nwire, Nsurf, *Iint, *Isurf, cint …  in DMDAGetWireBasketInterpolation()  local
58   Nint  = (m - 2 - istart) * (n - 2 - jstart) * (p - 2 - kstart);  in DMDAGetWireBasketInterpolation()
62   PetscCall(MatCreateSeqDense(MPI_COMM_SELF, Nint, 26, NULL, &Xint));  in DMDAGetWireBasketInterpolation()
139   PetscCall(PetscMalloc3(N, &II, Nint, &Iint, Nsurf, &Isurf));  in DMDAGetWireBasketInterpolation()
140   PetscCall(PetscMalloc2(Nint, &IIint, Nsurf, &IIsurf));  in DMDAGetWireBasketInterpolation()
158   PetscCheck(cint == Nint, PETSC_COMM_SELF, PETSC_ERR_PLIB, "cint != Nint");  in DMDAGetWireBasketInterpolation()
162   PetscCall(ISLocalToGlobalMappingApply(ltg, Nint, IIint, IIint));  in DMDAGetWireBasketInterpolation()
166   PetscCall(ISCreateGeneral(PETSC_COMM_SELF, Nint, Iint, PETSC_COPY_VALUES, &isint));  in DMDAGetWireBasketInterpolation()
198     PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, Nint, NULL, &x));  in DMDAGetWireBasketInterpolation()
200     PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, Nint, NULL, &b));  in DMDAGetWireBasketInterpolation()
203       PetscCall(VecPlaceArray(x, xint + i * Nint));  in DMDAGetWireBasketInterpolation()
204       PetscCall(VecPlaceArray(b, xint_tmp + i * Nint));  in DMDAGetWireBasketInterpolation()
219   for (i = 0; i < Nint; i++) {  in DMDAGetWireBasketInterpolation()
221     for (j = 0; j < 26; j++) tmp += rxint[i + j * Nint];  in DMDAGetWireBasketInterpolation()
319 …tComm((PetscObject)da), Ng, PETSC_DECIDE, PETSC_DECIDE, Ntotal, Nint + Nsurf, NULL, Nint + Nsurf, …  in DMDAGetWireBasketInterpolation()
325   PetscCall(MatSetValues(*P, Nint, IIint, 26, gl, rxint, INSERT_VALUES));  in DMDAGetWireBasketInterpolation()
367 …PetscInt               dim, i, j, k, m, n, p, dof, Nint, Nface, Nwire, Nsurf, *Iint, *Isurf, cint …  in DMDAGetFaceInterpolation()  local
406   Nint  = (m - 2 - istart) * (n - 2 - jstart) * (p - 2 - kstart);  in DMDAGetFaceInterpolation()
410   PetscCall(MatCreateSeqDense(MPI_COMM_SELF, Nint, 6, NULL, &Xint));  in DMDAGetFaceInterpolation()
459   PetscCall(PetscMalloc3(N, &II, Nint, &Iint, Nsurf, &Isurf));  in DMDAGetFaceInterpolation()
460   PetscCall(PetscMalloc2(Nint, &IIint, Nsurf, &IIsurf));  in DMDAGetFaceInterpolation()
478   PetscCheck(cint == Nint, PETSC_COMM_SELF, PETSC_ERR_PLIB, "cint != Nint");  in DMDAGetFaceInterpolation()
482   PetscCall(ISLocalToGlobalMappingApply(ltg, Nint, IIint, IIint));  in DMDAGetFaceInterpolation()
486   PetscCall(ISCreateGeneral(PETSC_COMM_SELF, Nint, Iint, PETSC_COPY_VALUES, &isint));  in DMDAGetFaceInterpolation()
520     PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, Nint, NULL, &x));  in DMDAGetFaceInterpolation()
522     PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF, 1, Nint, NULL, &b));  in DMDAGetFaceInterpolation()
525       PetscCall(VecPlaceArray(x, xint + i * Nint));  in DMDAGetFaceInterpolation()
526       PetscCall(VecPlaceArray(b, xint_tmp + i * Nint));  in DMDAGetFaceInterpolation()
541   for (i = 0; i < Nint; i++) {  in DMDAGetFaceInterpolation()
543     for (j = 0; j < 6; j++) tmp += rxint[i + j * Nint];  in DMDAGetFaceInterpolation()
608 …tComm((PetscObject)da), Ng, PETSC_DECIDE, PETSC_DECIDE, Ntotal, Nint + Nsurf, NULL, Nint, NULL, P)…  in DMDAGetFaceInterpolation()
614   PetscCall(MatSetValues(*P, Nint, IIint, 6, gl, rxint, INSERT_VALUES));  in DMDAGetFaceInterpolation()