| /petsc/src/dm/impls/da/ |
| H A D | dagetarray.c | 45 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecGetArray() local
51 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecGetArray()
57 if (N == xm * ym * zm * dof) { in DMDAVecGetArray()
59 gym = ym; in DMDAVecGetArray()
64 …ompatible with DMDA local sizes %" PetscInt_FMT " %" PetscInt_FMT, N, xm * ym * zm * dof, gxm * gy… in DMDAVecGetArray()
94 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecRestoreArray() local
100 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecRestoreArray()
106 if (N == xm * ym * zm * dof) { in DMDAVecRestoreArray()
108 gym = ym; in DMDAVecRestoreArray()
113 …ompatible with DMDA local sizes %" PetscInt_FMT " %" PetscInt_FMT, N, xm * ym * zm * dof, gxm * gy… in DMDAVecRestoreArray()
[all …]
|
| H A D | grvtk.c | 78 rloc[3] = info.ym; in DMDAVTKWriteAll_VTS()
88 PetscInt xs = -1, xm = -1, ys = -1, ym = -1, zs = -1, zm = -1, nnodes = 0; in DMDAVTKWriteAll_VTS() local
94 ym = grloc[r][3]; in DMDAVTKWriteAll_VTS()
97 nnodes = xm * ym * zm; in DMDAVTKWriteAll_VTS()
100 …FMT " %" PetscInt_FMT " %" PetscInt_FMT "\">\n", xs, xs + xm - 1, ys, ys + ym - 1, zs, zs + zm - 1… in DMDAVTKWriteAll_VTS()
151 PetscInt xs = -1, xm = -1, ys = -1, ym = -1, zs = -1, zm = -1, nnodes = 0; in DMDAVTKWriteAll_VTS() local
157 ym = grloc[r][3]; in DMDAVTKWriteAll_VTS()
160 nnodes = xm * ym * zm; in DMDAVTKWriteAll_VTS()
162 nnodes = info.xm * info.ym * info.zm; in DMDAVTKWriteAll_VTS()
180 for (j = 0; j < ym; j++) { in DMDAVTKWriteAll_VTS()
[all …]
|
| H A D | dadd.c | 43 PetscInt xm = 1, ym = 1, zm = 1; in DMDACreatePatchIS() local
178 ym = ne - ns; in DMDACreatePatchIS()
180 base = ms * ym * zm + ns * M * zm + ps * M * N; in DMDACreatePatchIS()
186 indices[idx] = l + dof * (base + si + xm * sj + xm * ym * sk); in DMDACreatePatchIS()
212 ym = ne - ns; in DMDACreatePatchIS()
214 base = ms * ym * zm + ns * M * zm + ps * M * N; in DMDACreatePatchIS()
244 indices[idx] = l + dof * (base + si + xm * sj + xm * ym * sk); in DMDACreatePatchIS()
273 PetscInt xm = 1, ym = 1, zm = 1; in DMDASubDomainDA_Private() local
290 N = info.ym; in DMDASubDomainDA_Private()
336 ym = N / n + ((N % n) > j); in DMDASubDomainDA_Private()
[all …]
|
| H A D | dalocal.c | 410 PetscInt j, i, xs, ys, xm, ym, zs, zm; in DMDAGetArray() local
432 ym = dd->Ye - dd->Ys; in DMDAGetArray()
449 ym = dd->ye - dd->ys; in DMDAGetArray()
466 … PetscCall(PetscMalloc((ym + 1) * sizeof(void *) + xm * ym * sizeof(PetscScalar), &iarray_start)); in DMDAGetArray()
468 ptr = (void **)(iarray_start + xm * ym * sizeof(PetscScalar) - ys * sizeof(void *)); in DMDAGetArray()
469 … for (j = ys; j < ys + ym; j++) ptr[j] = iarray_start + sizeof(PetscScalar) * (xm * (j - ys) - xs); in DMDAGetArray()
476 …PetscCall(PetscMalloc((zm + 1) * sizeof(void **) + (ym * zm + 1) * sizeof(void *) + xm * ym * zm *… in DMDAGetArray()
478 ptr = (void ***)(iarray_start + xm * ym * zm * sizeof(PetscScalar) - zs * sizeof(void *)); in DMDAGetArray()
479 bptr = (void **)(iarray_start + xm * ym * zm * sizeof(PetscScalar) + zm * sizeof(void **)); in DMDAGetArray()
480 for (i = zs; i < zs + zm; i++) ptr[i] = bptr + ((i - zs) * ym - ys); in DMDAGetArray()
[all …]
|
| /petsc/src/dm/impls/da/ftn-custom/ |
| H A D | zda1f90.c | 68 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray1_() local
71 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray1_()
81 if (N == xm * ym * zm * dof) { in dmdavecgetarray1_()
83 gym = ym; in dmdavecgetarray1_()
109 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray2_() local
112 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray2_()
122 if (N == xm * ym * zm * dof) { in dmdavecgetarray2_()
124 gym = ym; in dmdavecgetarray2_()
156 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in dmdavecgetarray3_() local
159 *ierr = DMDAGetCorners(*da, &xs, &ys, &zs, &xm, &ym, &zm); in dmdavecgetarray3_()
[all …]
|
| /petsc/src/dm/tutorials/ |
| H A D | ex6.c | 40 PetscInt i, start, end, *ifrom, x, y, xm, ym; in main() local
57 PetscCall(DMDAGetCorners(da, &x, &y, NULL, &xm, &ym, NULL)); in main()
79 PetscCall(ISCreateStride(PETSC_COMM_SELF, ym, xm + 2, xm + 2, &to)); in main()
83 PetscCall(PetscMalloc1(ym, &ifrom)); in main()
84 for (i = y; i < y + ym; i++) ifrom[i - y] = i; in main()
86 PetscCall(AOApplicationToPetsc(ao, ym, ifrom)); in main()
88 PetscCall(ISCreateGeneral(PETSC_COMM_WORLD, ym, ifrom, PETSC_OWN_POINTER, &from)); in main()
|
| H A D | ex10.c | 25 PetscScalar ym = -1.0, yp = 1.0; in main() local
34 dy = (yp - ym) / (Ny - 1); in main()
59 …) gauss_ptr[j][i] = PetscExpScalar(-(xm + i * dx) * (xm + i * dx) - (ym + j * dy) * (ym + j * dy)); in main()
|
| /petsc/src/dm/tests/ |
| H A D | ex2k.kokkos.cxx | 36 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in Update1() local
39 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in Update1()
44 for (j = ys; j < ys + ym; j++) { in Update1()
59 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in Update2() local
62 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in Update2()
65 #define Y2(k, j, i) y2[(k - zs) * ym * xm + (j - ys) * xm + (i - xs)] in Update2()
69 for (j = ys; j < ys + ym; j++) { in Update2()
84 PetscInt xm, ym, zm, xs, ys, zs, gxm, gym, gzm, gxs, gys, gzs; in main() local
111 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in main()
146 …xecutionSpace, Rank<3, Iterate::Right, Iterate::Right>>({zs, ys, xs}, {zs + zm, ys + ym, xs + xm}), in main()
|
| H A D | ex10k.kokkos.cxx | 30 PetscInt M = 5, N = 7, xm, ym, xs, ys; in main() local
52 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in main()
60 for (PetscInt j = ys; j < ys + ym; j++) { /* run on host */ in main()
70 …"init 1", MDRangePolicy<Rank<3, Iterate::Right, Iterate::Right>>({ys, xs, 0}, {ys + ym, xs + xm, d… in main()
82 for (PetscInt j = ys; j < ys + ym; j++) { in main()
94 …"stencil 1", MDRangePolicy<Rank<3, Iterate::Right, Iterate::Right>>({ys, xs, 0}, {ys + ym, xs + xm… in main()
125 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in main()
133 for (PetscInt j = ys; j < ys + ym; j++) { in main()
146 …"init 2", MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>({ys, xs}, {ys + ym, xs + xm}), KO… in main()
160 for (PetscInt j = ys; j < ys + ym; j++) { in main()
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex55k.kokkos.cxx | 41 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info->my; in FormFunctionLocalVec() local
61 …RangePolicy<Rank<2, Iterate::Right, Iterate::Right>>(exec, {ys, xs}, {ys + ym, xs + xm}), KOKKOS_L… in FormFunctionLocalVec()
112 PetscCall(PetscLogFlops(11.0 * info->ym * info->xm)); in FormFunctionLocalVec()
118 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info->my; in FormObjectiveLocalVec() local
140 …alVec", MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>(exec, {ys, xs}, {ys + ym, xs + xm}), in FormObjectiveLocalVec()
168 PetscCall(PetscLogFlops(12.0 * info->ym * info->xm)); in FormObjectiveLocalVec()
176 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info-… in FormJacobianLocalVec() local
191 hy = ym > 1 ? PetscRealPart(coords[ys + 1][xs].y) - PetscRealPart(coords[ys][xs].y) : 1.0; in FormJacobianLocalVec()
202 PetscCount ncoo = ((PetscCount)xm) * ((PetscCount)ym) * 5; in FormJacobianLocalVec()
208 for (j = ys; j < ys + ym; j++) { in FormJacobianLocalVec()
[all …]
|
| H A D | ex16.c | 602 PetscInt xm = info->xm, ym = info->ym, zm = info->zm; in FormJacobianLocal() local
622 yee = ys + ym; in FormJacobianLocal()
628 if (ys + ym == my) yee = ys + ym - 1; in FormJacobianLocal()
646 …if (i + ii >= xs && i + ii < xm + xs && j + jj >= ys && j + jj < ys + ym && k + kk >= zs && k + kk… in FormJacobianLocal()
677 for (j = ys; j < ys + ym; j++) { in FormJacobianLocal()
711 PetscInt xm = info->xm, ym = info->ym, zm = info->zm; in FormFunctionLocal() local
723 PetscCall(DMDAGetCorners(info->da, &xs, &ys, &zs, &xm, &ym, &zm)); in FormFunctionLocal()
727 for (j = ys; j < ys + ym; j++) { in FormFunctionLocal()
738 yee = ys + ym; in FormFunctionLocal()
744 if (ys + ym == my) yee = ys + ym - 1; in FormFunctionLocal()
[all …]
|
| /petsc/src/ts/tutorials/autodiff/ |
| H A D | adr_ex5adj_mf.cxx | 184 PetscInt i, j, xs, ys, xm, ym, Mx, My; in InitialConditions() local
202 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitialConditions()
207 for (j = ys; j < ys + ym; j++) { in InitialConditions()
228 PetscInt i, j, Mx, My, xs, ys, xm, ym; in InitializeLambda() local
236 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitializeLambda()
238 if (xs <= i && i < xs + xm && ys <= j && j < ys + ym) { in InitializeLambda()
251 PetscInt i, j, xs, ys, xm, ym; in IFunctionLocalPassive() local
265 ym = info->ym; in IFunctionLocalPassive()
268 for (j = ys; j < ys + ym; j++) { in IFunctionLocalPassive()
280 PetscCall(PetscLogFlops(16.0 * xm * ym)); in IFunctionLocalPassive()
[all …]
|
| H A D | adr_ex5adj.cxx | 296 PetscInt i, j, xs, ys, xm, ym, Mx, My; in InitialConditions() local
314 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitialConditions()
319 for (j = ys; j < ys + ym; j++) { in InitialConditions()
340 PetscInt i, j, Mx, My, xs, ys, xm, ym; in InitializeLambda() local
348 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitializeLambda()
350 if (xs <= i && i < xs + xm && ys <= j && j < ys + ym) { in InitializeLambda()
363 PetscInt i, j, xs, ys, xm, ym; in IFunctionLocalPassive() local
377 ym = info->ym; in IFunctionLocalPassive()
380 for (j = ys; j < ys + ym; j++) { in IFunctionLocalPassive()
392 PetscCall(PetscLogFlops(16.0 * xm * ym)); in IFunctionLocalPassive()
[all …]
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | minsurf2.c | 166 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunction() local
177 PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL)); in FormFunction()
188 for (j = ys; j < ys + ym; j++) { in FormFunction()
235 for (j = ys; j < ys + ym; j++) { in FormFunction()
249 for (j = ys; j < ys + ym; j++) { in FormFunction()
255 if (ys + ym == my) { /* top side */ in FormFunction()
267 if (ys + ym == my && xs + xm == mx) { in FormFunction()
268 d1 = (user->right[ym + 1] - user->right[ym]) * rhy; in FormFunction()
300 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunctionGradient() local
312 PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL)); in FormFunctionGradient()
[all …]
|
| /petsc/src/binding/petsc4py/demo/legacy/wrap-swig/ |
| H A D | Bratu3D.c | 24 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormInitGuess() local
55 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormInitGuess()
62 for (j=ys; j<ys+ym; j++) { in FormInitGuess()
84 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormFunction() local
120 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormFunction()
126 for (j=ys; j<ys+ym; j++) { in FormFunction()
155 PetscCall(PetscLogFlops(11.0*ym*xm)); in FormFunction()
161 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormJacobian() local
197 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormJacobian()
215 for (j=ys; j<ys+ym; j++) { in FormJacobian()
|
| /petsc/src/binding/petsc4py/demo/legacy/wrap-cython/ |
| H A D | Bratu3Dimpl.c | 24 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormInitGuess() local
56 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormInitGuess()
63 for (j=ys; j<ys+ym; j++) { in FormInitGuess()
85 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormFunction() local
124 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormFunction()
130 for (j=ys; j<ys+ym; j++) { in FormFunction()
159 PetscCall(PetscLogFlops(11.0*ym*xm)); in FormFunction()
165 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormJacobian() local
204 PetscCall(DMDAGetCorners(da,&xs,&ys,&zs,&xm,&ym,&zm)); in FormJacobian()
222 for (j=ys; j<ys+ym; j++) { in FormJacobian()
|
| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | reaction_diffusion.c | 31 PetscInt i, j, Mx, My, xs, ys, xm, ym; in RHSFunction() local
64 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in RHSFunction()
69 for (j = ys; j < ys + ym; j++) { in RHSFunction()
81 PetscCall(PetscLogFlops(16.0 * xm * ym)); in RHSFunction()
96 PetscInt i, j, Mx, My, xs, ys, xm, ym; in RHSJacobian() local
131 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in RHSJacobian()
143 for (j = ys; j < ys + ym; j++) { in RHSJacobian()
209 PetscCall(PetscLogFlops(19.0 * xm * ym)); in RHSJacobian()
239 PetscInt i, j, Mx, My, xs, ys, xm, ym; in IFunction() local
273 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in IFunction()
[all …]
|
| H A D | ex5adj_mf.c | 38 PetscInt i, j, Mx, My, xs, ys, xm, ym; in InitializeLambda() local
46 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitializeLambda()
48 if (xs <= i && i < xs + xm && ys <= j && j < ys + ym) { in InitializeLambda()
63 PetscInt i, j, Mx, My, xs, ys, xm, ym; in MyRHSMatMultTranspose() local
93 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in MyRHSMatMultTranspose()
96 for (j = ys; j < ys + ym; j++) { in MyRHSMatMultTranspose()
122 PetscInt i, j, Mx, My, xs, ys, xm, ym; in MyIMatMultTranspose() local
152 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in MyIMatMultTranspose()
155 for (j = ys; j < ys + ym; j++) { in MyIMatMultTranspose()
183 PetscInt i, j, Mx, My, xs, ys, xm, ym; in MyIMatMult() local
[all …]
|
| H A D | ex5opt_ic.c | 82 PetscInt i, j, xs, ys, xm, ym, Mx, My; in InitialConditions() local
94 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitialConditions()
97 for (j = ys; j < ys + ym; j++) { in InitialConditions()
114 PetscInt i, j, xs, ys, xm, ym, Mx, My; in PerturbedInitialConditions() local
126 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in PerturbedInitialConditions()
129 for (j = ys; j < ys + ym; j++) { in PerturbedInitialConditions()
146 PetscInt i, j, xs, ys, xm, ym, Mx, My; in PerturbedInitialConditions2() local
158 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in PerturbedInitialConditions2()
161 for (j = ys; j < ys + ym; j++) { in PerturbedInitialConditions2()
179 PetscInt i, j, xs, ys, xm, ym, Mx, My; in PerturbedInitialConditions3() local
[all …]
|
| H A D | ex5adj.c | 22 PetscInt i, j, Mx, My, xs, ys, xm, ym; in InitializeLambda() local
30 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitializeLambda()
32 if (xs <= i && i < xs + xm && ys <= j && j < ys + ym) { in InitializeLambda()
155 PetscInt i, j, xs, ys, xm, ym, Mx, My; in InitialConditions() local
173 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in InitialConditions()
178 for (j = ys; j < ys + ym; j++) { in InitialConditions()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex66.c | 66 PetscInt i, j, M, N, xm, ym, xs, ys; in ComputeRHS() local
81 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); /* Fine grid */ in ComputeRHS()
83 for (j = ys; j < ys + ym; j++) { in ComputeRHS()
100 PetscInt i, j, M, N, xm, ym, xs, ys; in ComputeJacobian() local
113 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in ComputeJacobian()
114 for (j = ys; j < ys + ym; j++) { in ComputeJacobian()
|
| /petsc/src/dm/impls/da/kokkos/ |
| H A D | dagetov.kokkos.cxx | 8 #define DMDA_VEC_GET_SHAPE(da, vec, xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, d… argument
10 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); \
15 if (N == xm * ym * zm * dof) { \
17 gym = ym; \
22 …ompatible with DMDA local sizes %" PetscInt_FMT " %" PetscInt_FMT, N, xm * ym * zm * dof, gxm * gy…
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()
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()
95 …PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, di… in DMDAVecGetKokkosOffsetView_Private() local
[all …]
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2f.F90 | 68 PetscInt ys, ym
93 PetscCall(DMDAGetCorners(dm, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr))
112 do j = ys, ys + ym - 1
206 do j = ys, ys + ym - 1
222 do j = ys, ys + ym - 1
229 if (ys + ym == my) then
243 if ((ys + ym == my) .and. (xs + xm == mx)) then
244 d1 = (right_v(1 + ym + 1) - right_v(1 + ym))*rhy
264 PetscCall(PetscLogFlops(70.0d0*xm*ym, ierr))
303 PetscInt ys, ym, gys, gym
[all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex13.c | 120 PetscInt i, j, Mx, My, xs, ys, xm, ym; in RHSFunction() local
149 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in RHSFunction()
152 for (j = ys; j < ys + ym; j++) { in RHSFunction()
169 PetscCall(PetscLogFlops(11.0 * ym * xm)); in RHSFunction()
202 for (j = info.ys; j < info.ys + info.ym; j++) { in RHSJacobian()
246 PetscInt i, j, xs, ys, xm, ym, Mx, My; in FormInitialSolution() local
260 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in FormInitialSolution()
263 for (j = ys; j < ys + ym; j++) { in FormInitialSolution()
|
| H A D | ex7.c | 127 PetscInt i, j, Mx, My, xs, ys, xm, ym; in FormFunction() local
160 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in FormFunction()
165 for (j = ys; j < ys + ym; j++) { in FormFunction()
185 PetscCall(PetscLogFlops(11.0 * ym * xm)); in FormFunction()
192 PetscInt i, j, xs, ys, xm, ym, Mx, My; in FormInitialSolution() local
210 PetscCall(DMDAGetCorners(da, &xs, &ys, NULL, &xm, &ym, NULL)); in FormInitialSolution()
215 for (j = ys; j < ys + ym; j++) { in FormInitialSolution()
|