Home
last modified time | relevance | path

Searched refs:gxm (Results 1 – 25 of 32) 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()
58 gxm = xm; in DMDAVecGetArray()
64 …gxm * gym * gzm * dof, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Vector local size %" PetscInt_FMT "… 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()
107 gxm = xm; in DMDAVecRestoreArray()
[all …]
/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_()
82 gxm = xm; in dmdavecgetarray1_()
88 } else if (N != gxm * gym * gzm * dof) { 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_()
123 gxm = xm; in dmdavecgetarray2_()
129 } else if (N != gxm * gym * gzm * dof) { in dmdavecgetarray2_()
135 gym = gxm; in dmdavecgetarray2_()
[all …]
/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()
179 PetscCall(PetscMalloc1(gxm, &user.Vt1)); in main()
180 PetscCall(PetscMalloc1(gxm, &user.c)); in main()
181 PetscCall(PetscMalloc1(gxm, &user.d)); in main()
187 for (i = 0; i < gxm; i++) { in main()
193 if (gxs + gxm == user.ms) user.Vt1[gxm - 1] = 0; in main()
235 for (j = 0; j < gxm; j++) user.Vt1[j] = x_array[j]; 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()
[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()
252 if (i > -1 && j < my - 1) vt = x[k + gxm]; in FormFunctionGradient()
266 ind = k + gxm; in FormFunctionGradient()
278 k = (j - gys) * gxm + i - gxs; in FormFunctionGradient()
282 if (i < mx && j > 0) vb = x[k - gxm]; in FormFunctionGradient()
[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
185 if (i > -1 .and. j < my - 1) vt = lx_v(k + 1 + gxm)
199 ind = k + gxm
212 k = (j - gys)*gxm + i - gxs
216 if (i < mx .and. j > 0) vb = lx_v(k + 1 - gxm)
[all …]
H A Dminsurf2.c166 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunction() local
178 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormFunction()
204 if (i + 1 == gxs + gxm) { /* right side */ in FormFunction()
300 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunctionGradient() local
313 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in FormFunctionGradient()
343 if (i + 1 == gxs + gxm) { /* right side */ in FormFunctionGradient()
358 if (j > gys && i + 1 < gxs + gxm) xrb = x[j - 1][i + 1]; in FormFunctionGradient()
506 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in QuadraticH() local
520 PetscCall(DMDAGetGhostCorners(user->dm, &gxs, &gys, NULL, &gxm, &gym, NULL)); in QuadraticH()
690 PetscInt xs, ys, xm, ym, gxs, gys, gxm, gym; in MSA_BoundaryConditions() local
[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)); \
16 gxm = xm; \
22 …gxm * gym * gzm * dof, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Vector local size %" PetscInt_FMT "…
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()
[all …]
/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)
134 xb = x_v(1 + row - gxm)
137 if (i + 1 == gxs + gxm) then !right side
148 xt = x_v(1 + row + gxm)
152 xlt = x_v(1 + row - 1 + gxm)
155 if ((j > gys) .and. (i + 1 < gxs + gxm)) then
156 xrb = x_v(1 + row + 1 - gxm)
208 d2 = (left_v(1 + j - ys + 1) - x_v(1 + (j - gys)*gxm))*rhx
[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()
251 xb = x[row - gxm]; in FormFunctionGradient()
254 if (i + 1 == gxs + gxm) { /* right side */ in FormFunctionGradient()
265 xt = x[row + gxm]; 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()
320 d2 = (left[j - ys + 1] - x[(j - gys) * gxm]) * rhx; in FormFunctionGradient()
334 d1 = (x[(j + 1 - gys) * gxm - 1] - right[j - ys + 1]) * rhx; in FormFunctionGradient()
[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()
272 col[k] = row - gxm; in FormFunctionGradient()
286 if (i + 1 < gxs + gxm) { in FormFunctionGradient()
294 col[k] = row + 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()
[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()
391 uyy = (two * u - x[row - gxm] - x[row + gxm]) * hxdhy; 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()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex14f.F90147 PetscInt gys, gxm, gym
182 …PetscCall(DMDAGetGhostCorners(da, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr…
188 row = (j - gys)*gxm + xs - gxs
199 uyy = (two*u - xx(row - gxm) - xx(row + gxm))*hxdhy
233 PetscInt gxs, gys, gxm, gym
270 …PetscCall(DMDAGetGhostCorners(da, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr…
289 row = (j - gys)*gxm + xs - gxs
298 col(1) = ltog(row - gxm)
306 col(5) = ltog(row + gxm)
/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()
H A Dex43.c27 …ank %d portion (with ghosts, %" PetscInt_FMT " x %" PetscInt_FMT ")\n", rank, info.gxm, info.gym)); in PrintVecWithGhosts()
29 for (i = info.gxs; i < info.gxs + info.gxm; i++) { in PrintVecWithGhosts()
/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()
31 (xs,xm,gxs,gxm) = self.getCorners()
/petsc/src/ts/tutorials/autodiff/adolc-utils/
H A Dinit.cxx111 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/ts/tutorials/autodiff/
H A Dadr_ex5adj_mf.cxx52 PetscInt gxm, gym; in main() local
109 PetscCall(DMDAGetGhostCorners(da, NULL, NULL, NULL, &gxm, &gym, NULL)); in main()
110 matctx.m = 2 * gxm * gym; in main()
111 matctx.n = 2 * gxm * gym; /* Number of dependent and independent variables */ in main()
291 PetscInt i, j, xs, ys, xm, ym, gxs, gys, gxm, gym; in IFunctionActive() local
308 gxm = info.gxm; in IFunctionActive()
337 u_c = new AField[info.gxm * info.gym]; in IFunctionActive()
338 f_c = new AField[info.gxm * info.gym]; in IFunctionActive()
357 for (i = gxs; i < gxs + gxm; i++) { in IFunctionActive()
384 for (i = gxs; i < gxs + gxm; i++) { in IFunctionActive()
H A Dadr_ex5adj.cxx81 PetscInt gxm, gym, i, dofs = 2, ctrl[3] = {0, 0, 0}; in main() local
135 PetscCall(DMDAGetGhostCorners(da, NULL, NULL, NULL, &gxm, &gym, NULL)); in main()
136 adctx->m = dofs * gxm * gym; in main()
137 adctx->n = dofs * gxm * gym; /* Number of dependent and independent variables */ in main()
403 PetscInt i, j, xs, ys, xm, ym, gxs, gys, gxm, gym; in IFunctionActive() local
420 gxm = info.gxm; in IFunctionActive()
449 u_c = new AField[info.gxm * info.gym]; in IFunctionActive()
450 f_c = new AField[info.gxm * info.gym]; in IFunctionActive()
469 for (i = gxs; i < gxs + gxm; i++) { in IFunctionActive()
496 for (i = gxs; i < gxs + gxm; i++) { in IFunctionActive()
[all …]
/petsc/src/snes/tutorials/
H A Dex5f.F9039 PetscInt xs, xe, xm, gxs, gxe, gxm
288 row = (j - gys)*gxm + xs - gxs - 1
305 col(1) = row - gxm
309 col(5) = row + gxm
446 call DMDAGetGhostCorners(da, gxs, gys, PETSC_NULL_INTEGER, gxm, gym, PETSC_NULL_INTEGER, ierr)
460 gxe = gxs + gxm - 1
H A Dex5f90t.F9040 PetscInt xs, xe, xm, gxs, gxe, gxm
423 row = (j - ctx%gys)*ctx%gxm + ctx%xs - ctx%gxs - 1
438 col(1) = row - ctx%gxm
442 col(5) = row + ctx%gxm
525 …PetscCallA(DMDAGetGhostCorners(ctx%da, ctx%gxs, ctx%gys, PETSC_NULL_INTEGER, ctx%gxm, ctx%gym, PET…
537 ctx%gxe = ctx%gxs + ctx%gxm - 1
H A Dex5f90.F9040 PetscInt xs, xe, xm, gxs, gxe, gxm
434 row = (j - ctx%gys)*ctx%gxm + ctx%xs - ctx%gxs - 1
449 col(1) = row - ctx%gxm
453 col(5) = row + ctx%gxm
537 …PetscCallA(DMDAGetGhostCorners(da, ctx%gxs, ctx%gys, PETSC_NULL_INTEGER, ctx%gxm, ctx%gym, PETSC_N…
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…
52 grd%gxe = grd%gxs + grd%gxm - 1
148 row = (j - grd%gys)*grd%gxm + grd%xs - grd%gxs - 1
163 col(1) = row(1) - grd%gxm
167 col(5) = row(1) + grd%gxm
/petsc/include/
H A Dpetscdmdatypes.h65 PetscInt gxm, gym, gzm; /* number of grid points on this processor including ghosts */ member
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A Dpetscdmda.pxi217 cdef PetscInt gxm=0, gym=0, gzm=0
220 &gxm, &gym, &gzm))
228 elif (n == gxm*gym*gzm*dof):
230 xm, ym, zm = gxm, gym, gzm

12