Lines Matching refs:Nsurf
19 …PetscInt dim, i, j, k, m, n, p, dof, Nint, Nface, Nwire, Nsurf, *Iint, *Isurf, cint … in DMDAGetWireBasketInterpolation() local
61 Nsurf = Nface + Nwire; in DMDAGetWireBasketInterpolation()
63 PetscCall(MatCreateSeqDense(MPI_COMM_SELF, Nsurf, 26, NULL, &Xsurf)); in DMDAGetWireBasketInterpolation()
77 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + Nsurf] = 1; in DMDAGetWireBasketInterpolation()
78 xsurf[cnt++ + 2 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
81 xsurf[cnt++ + 3 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
82 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 4 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
83 xsurf[cnt++ + 5 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
86 xsurf[cnt++ + 6 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
87 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 7 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
88 xsurf[cnt++ + 8 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
91 xsurf[cnt++ + 9 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
92 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 10 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
93 xsurf[cnt++ + 11 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
96 xsurf[cnt++ + 12 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
98 xsurf[cnt++ + 13 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
101 xsurf[cnt++ + 14 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
102 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 15 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
103 xsurf[cnt++ + 16 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
106 xsurf[cnt++ + 17 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
107 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 18 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
108 xsurf[cnt++ + 19 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
111 xsurf[cnt++ + 20 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
112 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 21 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
113 xsurf[cnt++ + 22 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
116 xsurf[cnt++ + 23 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
117 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 24 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
118 xsurf[cnt++ + 25 * Nsurf] = 1; in DMDAGetWireBasketInterpolation()
122 for (i = 0; i < Nsurf; i++) { in DMDAGetWireBasketInterpolation()
124 for (j = 0; j < 26; j++) tmp += xsurf[i + j * Nsurf]; in DMDAGetWireBasketInterpolation()
139 PetscCall(PetscMalloc3(N, &II, Nint, &Iint, Nsurf, &Isurf)); in DMDAGetWireBasketInterpolation()
140 PetscCall(PetscMalloc2(Nint, &IIint, Nsurf, &IIsurf)); in DMDAGetWireBasketInterpolation()
159 PetscCheck(csurf == Nsurf, PETSC_COMM_SELF, PETSC_ERR_PLIB, "csurf != Nsurf"); in DMDAGetWireBasketInterpolation()
163 PetscCall(ISLocalToGlobalMappingApply(ltg, Nsurf, IIsurf, IIsurf)); in DMDAGetWireBasketInterpolation()
167 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, Nsurf, Isurf, PETSC_COPY_VALUES, &issurf)); in DMDAGetWireBasketInterpolation()
319 …PetscObject)da), Ng, PETSC_DECIDE, PETSC_DECIDE, Ntotal, Nint + Nsurf, NULL, Nint + Nsurf, NULL, P… in DMDAGetWireBasketInterpolation()
328 PetscCall(MatSetValues(*P, Nsurf, IIsurf, 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
409 Nsurf = Nface + Nwire; in DMDAGetFaceInterpolation()
411 PetscCall(MatCreateSeqDense(MPI_COMM_SELF, Nsurf, 6, NULL, &Xsurf)); in DMDAGetFaceInterpolation()
424 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 0 * Nsurf] = 1; in DMDAGetFaceInterpolation()
428 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 1 * Nsurf] = 1; in DMDAGetFaceInterpolation()
430 xsurf[cnt++ + 2 * Nsurf] = 1; in DMDAGetFaceInterpolation()
432 xsurf[cnt++ + 3 * Nsurf] = 1; in DMDAGetFaceInterpolation()
434 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 4 * Nsurf] = 1; in DMDAGetFaceInterpolation()
437 for (i = 1; i < m - istart - 1; i++) xsurf[cnt++ + 5 * Nsurf] = 1; in DMDAGetFaceInterpolation()
441 for (i = 0; i < Nsurf; i++) { in DMDAGetFaceInterpolation()
443 for (j = 0; j < 6; j++) tmp += xsurf[i + j * Nsurf]; in DMDAGetFaceInterpolation()
459 PetscCall(PetscMalloc3(N, &II, Nint, &Iint, Nsurf, &Isurf)); in DMDAGetFaceInterpolation()
460 PetscCall(PetscMalloc2(Nint, &IIint, Nsurf, &IIsurf)); in DMDAGetFaceInterpolation()
479 PetscCheck(csurf == Nsurf, PETSC_COMM_SELF, PETSC_ERR_PLIB, "csurf != Nsurf"); in DMDAGetFaceInterpolation()
483 PetscCall(ISLocalToGlobalMappingApply(ltg, Nsurf, IIsurf, IIsurf)); in DMDAGetFaceInterpolation()
487 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, Nsurf, Isurf, PETSC_COPY_VALUES, &issurf)); in DMDAGetFaceInterpolation()
608 …jectComm((PetscObject)da), Ng, PETSC_DECIDE, PETSC_DECIDE, Ntotal, Nint + Nsurf, NULL, Nint, NULL,… in DMDAGetFaceInterpolation()
617 PetscCall(MatSetValues(*P, Nsurf, IIsurf, 6, gl, rxint, INSERT_VALUES)); in DMDAGetFaceInterpolation()