Home
last modified time | relevance | path

Searched refs:gxs (Results 1 – 25 of 38) sorted by relevance

12

/petsc/src/dm/impls/da/
H A Ddagetarray.c45 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecGetArray() local
52 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, &gzs, &gxm, &gym, &gzm)); in DMDAVecGetArray()
61 gxs = xs; in DMDAVecGetArray()
67 PetscCall(VecGetArray1d(vec, gxm * dof, gxs * dof, (PetscScalar **)array)); in DMDAVecGetArray()
69 PetscCall(VecGetArray2d(vec, gym, gxm * dof, gys, gxs * dof, (PetscScalar ***)array)); in DMDAVecGetArray()
71 … PetscCall(VecGetArray3d(vec, gzm, gym, gxm * dof, gzs, gys, gxs * dof, (PetscScalar ****)array)); in DMDAVecGetArray()
94 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecRestoreArray() local
101 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, &gzs, &gxm, &gym, &gzm)); in DMDAVecRestoreArray()
110 gxs = xs; in DMDAVecRestoreArray()
116 PetscCall(VecRestoreArray1d(vec, gxm * dof, gxs * dof, (PetscScalar **)array)); in DMDAVecRestoreArray()
[all …]
H A Dfdda.c261 PetscInt xs, ys, nx, ny, i, j, ii, gxs, gys, gnx, gny, m, n, M, N, dim, s, k, nc, col; in DMCreateColoring_DA_2d_MPIAIJ() local
278 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, NULL, &gnx, &gny, NULL)); in DMCreateColoring_DA_2d_MPIAIJ()
303 for (i = gxs; i < gxs + gnx; i++) { in DMCreateColoring_DA_2d_MPIAIJ()
325 …PetscInt xs, ys, nx, ny, i, j, gxs, gys, gnx, gny, m, n, p, dim, s, k, nc, col, zs, gzs, i… in DMCreateColoring_DA_3d_MPIAIJ() local
345 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, &gzs, &gnx, &gny, &gnz)); in DMCreateColoring_DA_3d_MPIAIJ()
370 for (i = gxs; i < gxs + gnx; i++) { in DMCreateColoring_DA_3d_MPIAIJ()
390 PetscInt xs, nx, i, i1, gxs, gnx, l, m, M, dim, s, nc, col; in DMCreateColoring_DA_1d_MPIAIJ() local
406 PetscCall(DMDAGetGhostCorners(da, &gxs, NULL, NULL, &gnx, NULL, NULL)); in DMCreateColoring_DA_1d_MPIAIJ()
441 for (i = gxs; i < gxs + gnx; i++) { in DMCreateColoring_DA_1d_MPIAIJ()
459 PetscInt xs, ys, nx, ny, i, j, ii, gxs, gys, gnx, gny, m, n, dim, s, k, nc; in DMCreateColoring_DA_2d_5pt_MPIAIJ() local
[all …]
/petsc/src/vec/vec/tutorials/
H A Dex9.c25 Vec lx, gx, gxs; in main() local
74 …Call(VecCreateGhostWithArray(PETSC_COMM_WORLD, nlocal, PETSC_DECIDE, nghost, ifrom, tarray, &gxs)); in main()
76 PetscCall(VecCreate(PETSC_COMM_WORLD, &gxs)); in main()
77 PetscCall(VecSetType(gxs, VECMPI)); in main()
78 PetscCall(VecSetSizes(gxs, nlocal, PETSC_DECIDE)); in main()
79 PetscCall(VecMPISetGhost(gxs, nghost, ifrom)); in main()
81 PetscCall(VecCreateGhost(PETSC_COMM_WORLD, nlocal, PETSC_DECIDE, nghost, ifrom, &gxs)); in main()
82 PetscCall(VecSet(gxs, 1.0)); in main()
83 if (rank == 1) PetscCall(VecSetValueLocal(gxs, 0, 2.0, INSERT_VALUES)); in main()
84 PetscCall(VecAssemblyBegin(gxs)); in main()
[all …]
H A Dex14f.F9023 Vec lx, gx, gxs
71 …eGhostBlockWithArray(PETSC_COMM_WORLD, bs, nlocal, PETSC_DECIDE, nghost, ifrom, tarray, gxs, ierr))
73 …PetscCallA(VecCreateGhostBlock(PETSC_COMM_WORLD, bs, nlocal, PETSC_DECIDE, nghost, ifrom, gxs, ier…
78 PetscCallA(VecDuplicate(gxs, gx, ierr))
79 PetscCallA(VecDestroy(gxs, ierr))
H A Dex9f.F9023 Vec lx, gx, gxs
69 …(VecCreateGhostWithArray(PETSC_COMM_WORLD, nlocal, PETSC_DECIDE, nghost, ifrom, tarray, gxs, ierr))
71 PetscCallA(VecCreateGhost(PETSC_COMM_WORLD, nlocal, PETSC_DECIDE, nghost, ifrom, gxs, ierr))
76 PetscCallA(VecDuplicate(gxs, gx, ierr))
77 PetscCallA(VecDestroy(gxs, ierr))
/petsc/src/dm/impls/da/ftn-custom/
H A Dzda1f90.c68 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray1_() local
73 *ierr = DMDAGetGhostCorners(*da, &gxs, &gys, &gzs, &gxm, &gym, &gzm); in dmdavecgetarray1_()
85 gxs = xs; in dmdavecgetarray1_()
94 *ierr = F90Array1dCreate(aa, MPIU_SCALAR, gxs, gxm, a PETSC_F90_2PTR_PARAM(ptrd)); in dmdavecgetarray1_()
109 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray2_() local
114 *ierr = DMDAGetGhostCorners(*da, &gxs, &gys, &gzs, &gxm, &gym, &gzm); in dmdavecgetarray2_()
126 gxs = xs; in dmdavecgetarray2_()
134 gys = gxs; in dmdavecgetarray2_()
136 gxs = 0; in dmdavecgetarray2_()
141 *ierr = F90Array2dCreate(aa, MPIU_SCALAR, gxs, gxm, gys, gym, a PETSC_F90_2PTR_PARAM(ptrd)); in dmdavecgetarray2_()
[all …]
/petsc/src/ts/tutorials/
H A Dex8.py12 gxs, gxm = self.da.getGhostCorners()
13 return xs[0], xm[0], gxs[0], gxm[0]
16 (xs,xm,gxs,gxm) = self.getCorners()
25 if lxe == mx: f[mx-1] = xx[mx-1-gxs]/hx; lxe-=1
27 f[i] = xxdot[i-gxs] + (2.0*xx[i-gxs] - xx[i-1-gxs] - xx[i+1-gxs])/hx - hx*math.exp(xx[i-gxs])
31 (xs,xm,gxs,gxm) = self.getCorners()
41 P.setValues([i],[i-1,i,i+1],[-1.0/hx,2.0/hx-hx*math.exp(xx[i-gxs])+shift,-1.0/hx])
H A Dex22f_mf.F9053 PetscInt xs, xe, gxs, gxe, dof, gdof
105 PetscCallA(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
107 gdof = i2*(gxe - gxs + 1)
147 PetscCallA(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
149 gdof = i2*(gxe - gxs + 1)
165 subroutine GetLayout(da, mx, xs, xe, gxs, gxe, ierr) argument
170 PetscInt mx, xs, xe, gxs, gxe
175 …PetscCall(DMDAGetGhostCorners(da, gxs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, gxm, PETSC_NULL_INT…
177 gxs = gxs + 1
179 gxe = gxs + gxm - 1
[all …]
H A Dex22f.F90137 subroutine GetLayout(da, mx, xs, xe, gxs, gxe, ierr) argument
142 PetscInt mx, xs, xe, gxs, gxe
147 …PetscCall(DMDAGetGhostCorners(da, gxs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, gxm, PETSC_NULL_INT…
149 gxs = gxs + 1
151 gxe = gxs + gxm - 1
154 subroutine FormIFunctionLocal(mx, xs, xe, gxs, gxe, x, xdot, f, a, k, s, ierr) argument
156 PetscInt mx, xs, xe, gxs, gxe
180 PetscInt mx, xs, xe, gxs, gxe
184 PetscCall(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
191 PetscCall(FormIFunctionLocal(mx, xs, xe, gxs, gxe, xx, xxdot, ff, ctx%a, ctx%k, ctx%s, ierr))
[all …]
/petsc/src/dm/impls/da/kokkos/
H A Ddagetov.kokkos.cxx8 #define DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, d… argument
11 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, &gzs, &gxm, &gym, &gzm)); \
19 gxs = xs; \
29 …PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, di… in DMDAVecGetKokkosOffsetView_Private() local
36 DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof); in DMDAVecGetKokkosOffsetView_Private()
41 *ov = PetscScalarKokkosOffsetView1DType<MemorySpace>(kv.data(), {gxs * dof}, {(gxs + gxm) * dof}); in DMDAVecGetKokkosOffsetView_Private()
63 …PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, … in DMDAVecGetKokkosOffsetView() local
70 DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof); in DMDAVecGetKokkosOffsetView()
73 …*ov = ConstPetscScalarKokkosOffsetView1DType<MemorySpace>(kv.data(), {gxs * dof}, {(gxs + gxm) * d… in DMDAVecGetKokkosOffsetView()
95 …PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, di… in DMDAVecGetKokkosOffsetView_Private() local
[all …]
/petsc/src/ts/tutorials/autodiff/adolc-utils/
H A Dinit.cxx87 PetscInt gxs; in GiveGhostPoints1d() local
90 PetscCall(DMDAGetGhostCorners(da, &gxs, NULL, NULL, NULL, NULL, NULL)); in GiveGhostPoints1d()
91 *a1d -= gxs; in GiveGhostPoints1d()
111 PetscInt gxs, gys, gxm, gym; in GiveGhostPoints2d() local
114 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, NULL, &gxm, &gym, NULL)); in GiveGhostPoints2d()
115 for (PetscInt j = 0; j < gym; j++) (*a2d)[j] = cgs + j * gxm - gxs; in GiveGhostPoints2d()
H A Dmatfree.cxx61 for (i = info.gxs; i < info.gxs + info.gxm; i++) { in PetscAdolcIJacobianVectorProduct()
77 for (i = info.gxs; i < info.gxs + info.gxm; i++) { in PetscAdolcIJacobianVectorProduct()
143 for (i = info.gxs; i < info.gxs + info.gxm; i++) { in PetscAdolcIJacobianVectorProductIDMass()
211 for (i = info.gxs; i < info.gxs + info.gxm; i++) { in PetscAdolcIJacobianTransposeVectorProduct()
231 for (i = info.gxs; i < info.gxs + info.gxm; i++) { in PetscAdolcIJacobianTransposeVectorProduct()
297 for (i = info.gxs; i < info.gxs + info.gxm; i++) { in PetscAdolcIJacobianTransposeVectorProductIDMass()
/petsc/src/tao/complementarity/tutorials/
H A Dblackscholes.c126 PetscInt xs, xm, gxs, gxm; in main() local
169 PetscCall(DMDAGetGhostCorners(user.dm, &gxs, NULL, NULL, &gxm, NULL, NULL)); in main()
188 sval = (gxs + i) * user.ds; in main()
193 if (gxs + gxm == user.ms) user.Vt1[gxm - 1] = 0; in main()
217 for (i = 0; i < xm; i++) x_array[i] = user.Vt1[i - gxs + xs]; in main()
311 PetscInt i, xs, xm, gxs, gxm; in FormConstraints() local
321 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, NULL, NULL, &gxm, NULL, NULL)); in FormConstraints()
333 if (gxs == 0) { in FormConstraints()
344 if (gxs + gxm == ms) { in FormConstraints()
386 PetscInt gxs, gxm; in FormJacobian() local
[all …]
/petsc/src/tao/unconstrained/tutorials/
H A Deptorsion2.c159 PetscInt xs, ys, xm, ym, gxm, gym, gxs, gys, xe, ye; in FormInitialGuess() local
165 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormInitialGuess()
173 k = (j - gys) * gxm + i - gxs; in FormInitialGuess()
201 PetscInt xs, ys, xm, ym, gxm, gym, gxs, gys; in FormFunctionGradient() local
229 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormFunctionGradient()
246 k = (j - gys) * gxm + i - gxs; in FormFunctionGradient()
278 k = (j - gys) * gxm + i - gxs; in FormFunctionGradient()
331 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormHessian() local
343 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormHessian()
347 row = (j - gys) * gxm + (i - gxs); in FormHessian()
[all …]
H A Deptorsion2f.F9058 PetscInt gxm, gym, gxs, gys
67 …PetscCall(DMDAGetGhostCorners(dm, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr…
75 k = (j - gys)*gxm + i - gxs
124 PetscInt gxs, gys, gxm, gym
146 …PetscCall(DMDAGetGhostCorners(dm, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr…
179 k = (j - gys)*gxm + i - gxs
212 k = (j - gys)*gxm + i - gxs
267 PetscInt xs, xm, gxs, gxm
276 …PetscCall(DMDAGetGhostCorners(dm, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr…
280 row = (j - gys)*gxm + (i - gxs)
[all …]
/petsc/src/snes/tutorials/
H A Dex5f.F9039 PetscInt xs, xe, xm, gxs, gxe, gxm
164 PetscScalar x(gxs:gxe, gys:gye)
250 PetscScalar x(gxs:gxe, gys:gye)
288 row = (j - gys)*gxm + xs - gxs - 1
446 call DMDAGetGhostCorners(da, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr)
454 gxs = gxs + 1
460 gxe = gxs + gxm - 1
H A Dex5f90t.F9040 PetscInt xs, xe, xm, gxs, gxe, gxm
224 PetscScalar x(ctx%gxs:ctx%gxe, ctx%gys:ctx%gye)
389 PetscScalar x(ctx%gxs:ctx%gxe, ctx%gys:ctx%gye)
423 row = (j - ctx%gys)*ctx%gxm + ctx%xs - ctx%gxs - 1
525 …PetscCallA(DMDAGetGhostCorners(ctx%da, ctx%gxs, ctx%gys, PETSC_NULL_INTEGER, ctx%gxm, ctx%gym, PET…
531 ctx%gxs = ctx%gxs + 1
537 ctx%gxe = ctx%gxs + ctx%gxm - 1
H A Dex5f90.F9040 PetscInt xs, xe, xm, gxs, gxe, gxm
234 PetscScalar x(ctx%gxs:ctx%gxe, ctx%gys:ctx%gye)
400 PetscScalar x(ctx%gxs:ctx%gxe, ctx%gys:ctx%gye)
434 row = (j - ctx%gys)*ctx%gxm + ctx%xs - ctx%gxs - 1
537 …PetscCallA(DMDAGetGhostCorners(da, ctx%gxs, ctx%gys, PETSC_NULL_INTEGER, ctx%gxm, ctx%gym, PETSC_N…
543 ctx%gxs = ctx%gxs + 1
549 ctx%gxe = ctx%gxs + ctx%gxm - 1
/petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/
H A DBratu2D.F9028 PetscInt mx, xs, xe, xm, gxs, gxe, gxm
42 …PetscCall(DMDAGetGhostCorners(da, grd%gxs, grd%gys, PETSC_NULL_INTEGER, grd%gxm, grd%gym, PETSC_NU…
46 grd%gxs = grd%gxs + 1
52 grd%gxe = grd%gxs + grd%gxm - 1
90 PetscScalar x(grd%gxs:grd%gxe, grd%gys:grd%gye)
128 PetscScalar x(grd%gxs:grd%gxe, grd%gys:grd%gye)
148 row = (j - grd%gys)*grd%gxm + grd%xs - grd%gxs - 1
/petsc/src/dm/tests/
H A Dex2k.kokkos.cxx36 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in Update1() local
40 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, &gzs, &gxm, &gym, &gzm)); in Update1()
59 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in Update2() local
63 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, &gzs, &gxm, &gym, &gzm)); in Update2()
64 #define X2(k, j, i) x2[(k - gzs) * gym * gxm + (j - gys) * gxm + (i - gxs)] in Update2()
84 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in main() local
112 PetscCall(DMDAGetGhostCorners(da, &gxs, &gys, &gzs, &gxm, &gym, &gzm)); in main()
/petsc/src/tao/bound/tutorials/
H A Dplate2f.F9067 PetscInt gxs, gxm
94 …PetscCall(DMDAGetGhostCorners(dm, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr…
114 row = (j - gys)*gxm + (i - gxs)
137 if (i + 1 == gxs + gxm) then !right side
151 if ((i > gxs) .and. (j + 1 < gys + gym)) then
155 if ((j > gys) .and. (i + 1 < gxs + gxm)) then
216 d3 = (bottom_v(1 + i - xs + 1) - x_v(1 + i - gxs))*rhy
231 d1 = (x_v(1 + (gym - 1)*gxm + i - gxs) - top_v(1 + i - xs + 1))*rhy
302 PetscInt xs, xm, gxs, gxm
325 …PetscCall(DMDAGetGhostCorners(dm, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr…
[all …]
H A Djbearing2.c175 PetscInt nx, ny, xs, xm, gxs, gxm, ys, ym, gys, gym; in ComputeB() local
192 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in ComputeB()
212 PetscInt col[5], row, nx, ny, xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunctionGradient() local
241 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormFunctionGradient()
262 row = (j - gys) * gxm + (i - gxs); in FormFunctionGradient()
276 if (i > gxs) { in FormFunctionGradient()
286 if (i + 1 < gxs + gxm) { in FormFunctionGradient()
327 PetscInt col[5], row, nx, ny, xs, xm, gxs, gxm, ys, ym, gys, gym; in FormHessian() local
348 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormHessian()
373 row = (j - gys) * gxm + (i - gxs); in FormHessian()
[all …]
H A Dplate2.c202 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunctionGradient() local
215 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormFunctionGradient()
235 row = (j - gys) * gxm + (i - gxs); in FormFunctionGradient()
254 if (i + 1 == gxs + gxm) { /* right side */ in FormFunctionGradient()
268 if (i > gxs && j + 1 < gys + gym) xlt = x[row - 1 + gxm]; in FormFunctionGradient()
269 if (j > gys && i + 1 < gxs + gxm) xrb = x[row + 1 - gxm]; in FormFunctionGradient()
327 d3 = (bottom[i - xs + 1] - x[i - gxs]) * rhy; in FormFunctionGradient()
341 d1 = (x[(gym - 1) * gxm + i - gxs] - top[i - xs + 1]) * rhy; in FormFunctionGradient()
420 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym, col[7]; in FormHessian() local
440 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormHessian()
[all …]
/petsc/src/ksp/ksp/tests/
H A Dex14.c279 PetscInt i, j, row, mx, my, xs, ys, xm, ym, gxm, gym, gxs, gys; in FormInitialGuess() local
307 PetscCall(DMDAGetGhostCorners(user->da, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormInitialGuess()
315 row = i - gxs + (j - gys) * gxm; in FormInitialGuess()
343 PetscInt i, j, row, mx, my, xs, ys, xm, ym, gxs, gys, gxm, gym; in ComputeFunction() local
376 PetscCall(DMDAGetGhostCorners(user->da, &gxs, &gys, NULL, &gxm, &gym, NULL)); in ComputeFunction()
382 row = (j - gys) * gxm + xs - gxs - 1; in ComputeFunction()
427 PetscInt xs, ys, xm, ym, gxs, gys, gxm, gym, grow; in ComputeJacobian() local
458 PetscCall(DMDAGetGhostCorners(user->da, &gxs, &gys, NULL, &gxm, &gym, NULL)); in ComputeJacobian()
479 row = (j - gys) * gxm + xs - gxs - 1; in ComputeJacobian()
/petsc/src/snes/utils/ftn-custom/
H A Dzdmdasnesf.c30 PetscCallFortranVoidFunction((*func)(info, &in[info->dof * info->gxs], &A, &m, ctx, &ierr)); in sourlj1d()
42 …PetscCallFortranVoidFunction((*func)(info, &in[info->gys][info->dof * info->gxs], &A, &m, ctx, &ie… in sourlj2d()
54 …PetscCallFortranVoidFunction((*func)(info, &in[info->gzs][info->gys][info->dof * info->gxs], &A, &… in sourlj3d()
92 …PetscCallFortranVoidFunction((*func)(info, &in[info->dof * info->gxs], &out[info->dof * info->xs],… in sourlf1d()
104 …PetscCallFortranVoidFunction((*func)(info, &in[info->gys][info->dof * info->gxs], &out[info->ys][i… in sourlf2d()
116 …PetscCallFortranVoidFunction((*func)(info, &in[info->gzs][info->gys][info->dof * info->gxs], &out[… in sourlf3d()

12