Home
last modified time | relevance | path

Searched refs:xs (Results 1 – 25 of 162) sorted by relevance

1234567

/petsc/src/ts/tutorials/
H A Dex22f_mf.F9053 PetscInt xs, xe, gxs, gxe, dof, gdof
105 PetscCallA(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
106 dof = i2*(xe - xs + 1)
147 PetscCallA(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
148 dof = i2*(xe - xs + 1)
150 call SaveSolutionToDisk(da, X, gdof, xs, xe)
165 subroutine GetLayout(da, mx, xs, xe, gxs, gxe, ierr) argument
170 PetscInt mx, xs, xe, gxs, gxe
174 …PetscCall(DMDAGetCorners(da, xs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, xm, PETSC_NULL_INTEGER, P…
176 xs = xs + 1
[all …]
H A Dex22f.F90137 subroutine GetLayout(da, mx, xs, xe, gxs, gxe, ierr) argument
142 PetscInt mx, xs, xe, gxs, gxe
146 …PetscCall(DMDAGetCorners(da, xs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, xm, PETSC_NULL_INTEGER, P…
148 xs = xs + 1
150 xe = xs + xm - 1
154 subroutine FormIFunctionLocal(mx, xs, xe, gxs, gxe, x, xdot, f, a, k, s, ierr) argument
156 PetscInt mx, xs, xe, gxs, gxe
157 PetscScalar x(2, xs:xe)
158 PetscScalar xdot(2, xs:xe)
159 PetscScalar f(2, xs:xe)
[all …]
H A Dex8.py11 xs, xm = self.da.getCorners()
13 return xs[0], xm[0], gxs[0], gxm[0]
16 (xs,xm,gxs,gxm) = self.getCorners()
23 lxs, lxe = xs, xs+xm
24 if xs == 0: f[0] = xx[0]/hx; lxs+=1;
31 (xs,xm,gxs,gxm) = self.getCorners()
37 lxs, lxe = xs, xs+xm
38 if xs == 0: P.setValues([0],[0],1.0/hx); lxs+=1
H A Dex32.c41 PetscInt lsize, M, xs, xm, i; in InitialConditions() local
53 PetscCall(DMDAGetCorners(da, &xs, 0, 0, &xm, 0, 0)); in InitialConditions()
55 for (i = xs; i < xs + xm; i++) { in InitialConditions()
60 app->sw[i - xs] = 1; in InitialConditions()
103 PetscInt M, xs, xm, i; in IFunction() local
111 PetscCall(DMDAGetCorners(da, &xs, 0, 0, &xm, 0, 0)); in IFunction()
123 for (i = xs; i < xs + xm; i++) { in IFunction()
129 f[i] = (u[i + 1] - 2 * u[i] + u[i - 1]) / h2 + app->sw[i - xs] * (-app->A * u[i]) - udot[i]; in IFunction()
149 PetscInt M, xs, xm, i; in IJacobian() local
156 PetscCall(DMDAGetCorners(da, &xs, 0, 0, &xm, 0, 0)); in IJacobian()
[all …]
H A Dex50.c86 PetscInt i, xs, xm, ind, j, lenglob; in main() local
143 PetscCall(DMDAGetCorners(appctx.da, &xs, NULL, NULL, &xm, NULL, NULL)); in main()
149 xs = xs / (appctx.param.N - 1); in main()
156 for (i = xs; i < xs + xm; i++) { in main()
260 PetscInt i, xs, xn; in TrueSolution() local
265 PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL)); in TrueSolution()
266 for (i = xs; i < xs + xn; i++) { in TrueSolution()
325 PetscInt i, j, xs, xn; in MatMult_Laplacian() local
346 PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL)); in MatMult_Laplacian()
348 …for (j = xs; j < xs + xn; j += appctx->param.N - 1) PetscCallBLAS("BLASgemv", BLASgemv_("N", &n, &… in MatMult_Laplacian()
[all …]
/petsc/src/dm/impls/da/
H A Dda1.c56 …(viewer, "X range of indices: %" PetscInt_FMT " %" PetscInt_FMT "\n", info.xs, info.xs + info.xm)); in DMView_DA_1d()
95 xmin = dd->xs / dd->w; in DMView_DA_1d()
138 PetscInt i, *idx, nn, left, xs, xe, x, Xs, Xe, start, m, IXs, IXe; in DMSetUp_DA_1D() local
165 xs = rank * M / m; in DMSetUp_DA_1D()
166 x = (rank + 1) * M / m - xs; in DMSetUp_DA_1D()
169 if (M / m == x) xs = rank * x; in DMSetUp_DA_1D()
170 else xs = rank * (x - 1) + (M + rank) % (x * m); in DMSetUp_DA_1D()
174 if (rank >= (M % m)) xs = (rank * (M / m) + M % m); in DMSetUp_DA_1D()
175 else xs = rank * (M / m) + rank; in DMSetUp_DA_1D()
177 PetscCallMPI(MPI_Allgather(&xs, 1, MPIU_INT, dd->lx, 1, MPIU_INT, comm)); in DMSetUp_DA_1D()
[all …]
H A Ddagetelem.c6 PetscInt i, xs, xe, Xs, Xe; in DMDAGetElements_1D() local
14 PetscCall(DMDAGetCorners(dm, &xs, NULL, NULL, &xe, NULL, NULL)); in DMDAGetElements_1D()
16 xe += xs; in DMDAGetElements_1D()
18 if (xs != Xs) xs -= 1; in DMDAGetElements_1D()
19 da->ne = 1 * (xe - xs - 1); in DMDAGetElements_1D()
21 for (i = xs; i < xe - 1; i++) { in DMDAGetElements_1D()
27 corners[0] = (xs - Xs); in DMDAGetElements_1D()
40 PetscInt i, xs, xe, Xs, Xe; in DMDAGetElements_2D() local
65 PetscCall(DMDAGetCorners(dm, &xs, &ys, NULL, &xe, &ye, NULL)); in DMDAGetElements_2D()
67 xe += xs; in DMDAGetElements_2D()
[all …]
H A Ddaltol.c32 left = dd->xs - dd->Xs; in DMLocalToLocalCreate_DA()
33 PetscCall(PetscMalloc1(dd->xe - dd->xs, &idx)); in DMLocalToLocalCreate_DA()
34 for (j = 0; j < dd->xe - dd->xs; j++) idx[j] = left + j; in DMLocalToLocalCreate_DA()
36 left = dd->xs - dd->Xs; in DMLocalToLocalCreate_DA()
39 PetscCall(PetscMalloc1((dd->xe - dd->xs) * (up - down), &idx)); in DMLocalToLocalCreate_DA()
42 for (j = 0; j < dd->xe - dd->xs; j++) idx[count++] = left + i * (dd->Xe - dd->Xs) + j; in DMLocalToLocalCreate_DA()
45 left = dd->xs - dd->Xs; in DMLocalToLocalCreate_DA()
50 count = (dd->xe - dd->xs) * (top - bottom) * (up - down); in DMLocalToLocalCreate_DA()
55 …for (k = 0; k < dd->xe - dd->xs; k++) idx[count++] = (left + j * (dd->Xe - dd->Xs)) + i * (dd->Xe … in DMLocalToLocalCreate_DA()
H A Ddagetarray.c45 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()
61 gxs = xs; 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()
110 gxs = xs; in DMDAVecRestoreArray()
169 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecGetArrayWrite() local
179 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xm, &ym, &zm)); in DMDAVecGetArrayWrite()
189 gxs = xs; in DMDAVecGetArrayWrite()
221 PetscInt xs, ys, zs, xm, ym, zm, gxs, gys, gzs, gxm, gym, gzm, N, dim, dof; in DMDAVecRestoreArrayWrite() local
[all …]
H A Ddageometry.c41 PetscInt xs, xe, Xs, Xe; in DMGetLocalBoundingBox_DA() local
52 PetscCall(DMDAGetCorners(da, &xs, &ys, &zs, &xe, &ye, &ze)); in DMGetLocalBoundingBox_DA()
56 xe += xs; in DMGetLocalBoundingBox_DA()
63 if (xs != Xs && Xs >= 0) xs -= 1; in DMGetLocalBoundingBox_DA()
79 c0 = (xs - Xs); in DMGetLocalBoundingBox_DA()
83 c0 = (xs - Xs) + (ys - Ys) * (Xe - Xs); in DMGetLocalBoundingBox_DA()
87 c0 = (xs - Xs) + (ys - Ys) * (Xe - Xs) + (zs - Zs) * (Xe - Xs) * (Ye - Ys); in DMGetLocalBoundingBox_DA()
118 cs[0] = xs; in DMGetLocalBoundingBox_DA()
134 PetscInt xs, xe, mxlocal; in private_DMDALocatePointsIS_2D_Regular() local
143 xs = cs[0]; in private_DMDALocatePointsIS_2D_Regular()
[all …]
/petsc/src/dm/tutorials/
H A Dex11f90.F9013 PetscInt m, n, p, dof, s, i, j, k, xs, xl
38 …PetscCallA(DMDAGetCorners(ada, xs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, xl, PETSC_NULL_INTEGER,…
40 do i = xs, xs + xl - 1
53 PetscCallA(DMDAGetCorners(ada, xs, ys, PETSC_NULL_INTEGER, xl, yl, PETSC_NULL_INTEGER, ierr))
55 do i = xs, xs + xl - 1
76 PetscCallA(DMDAGetCorners(ada, xs, ys, zs, xl, yl, zl, ierr))
78 do i = xs, xs + xl - 1
99 …PetscCallA(DMDAGetCorners(ada, xs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, xl, PETSC_NULL_INTEGER,…
101 do i = xs, xs + xl - 1
126 PetscCallA(DMDAGetCorners(ada, xs, ys, PETSC_NULL_INTEGER, xl, yl, PETSC_NULL_INTEGER, ierr))
[all …]
/petsc/src/tao/unconstrained/tutorials/
H A Dminsurf2.c166 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()
189 for (i = xs; i < xs + xm; i++) { in FormFunction()
199 xb = user->bottom[i - xs + 1]; in FormFunction()
211 xt = user->top[i - xs + 1]; in FormFunction()
234 if (xs == 0) { /* left side */ in FormFunction()
242 for (i = xs; i < xs + xm; i++) { in FormFunction()
243 d2 = (user->bottom[i + 1 - xs] - user->bottom[i - xs + 2]) * rhx; in FormFunction()
244 d3 = (user->bottom[i - xs + 1] - x[0][i]) * rhy; in FormFunction()
248 if (xs + xm == mx) { /* right side */ in FormFunction()
[all …]
H A Dburgers_spectral.c100 PetscInt i, xs, xm, ind, j, lenglob; in main() local
161 PetscCall(DMDAGetCorners(appctx.da, &xs, NULL, NULL, &xm, NULL, NULL)); in main()
167 xs = xs / (appctx.param.N - 1); in main()
174 for (i = xs; i < xs + xm; i++) { in main()
298 PetscInt i, xs, xn; in InitialConditions() local
303 PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL)); in InitialConditions()
304 …for (i = xs; i < xs + xn; i++) s[i] = 2.0 * appctx->param.mu * PETSC_PI * PetscSinScalar(PETSC_PI … in InitialConditions()
326 PetscInt i, xs, xn; in TrueSolution() local
331 PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL)); in TrueSolution()
332 …for (i = xs; i < xs + xn; i++) s[i] = 2.0 * appctx->param.mu * PETSC_PI * PetscSinScalar(PETSC_PI … in TrueSolution()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex69.c24 PetscInt j, xs, xn; in ComputeSolution() local
30 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xn, NULL, NULL)); in ComputeSolution()
35 for (j = xs; j < xs + xn; j++) { in ComputeSolution()
50 PetscInt i, j, xs, xn, n = gll->n; in ComputeRhs() local
56 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xn, NULL, NULL)); in ComputeRhs()
57 xs = xs / (n - 1); in ComputeRhs()
65 for (j = xs; j < xs + xn; j++) { in ComputeRhs()
91 PetscInt N = 80, n, q = 8, xs, xn, j, l; in main() local
157 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xn, NULL, NULL)); in main()
158 xs = xs / (n - 1); in main()
[all …]
H A Dex44f.F9049 PetscInt xs, xm, i, mx
53 …PetscCall(DMDAGetCorners(da, xs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, xm, PETSC_NULL_INTEGER, P…
56 do i = xs, xs + xm - 1
70 PetscInt xs, xm, i, mx
75 …PetscCall(DMDAGetCorners(da, xs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, xm, PETSC_NULL_INTEGER, P…
77 do i = xs, xs + xm - 1
/petsc/src/ts/tutorials/phasefield/
H A Dheat.c153 PetscInt i, Mx, xs, xm; in FormFunction() local
185 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xm, NULL, NULL)); in FormFunction()
190 for (i = xs; i < xs + xm; i++) { in FormFunction()
207 PetscInt i, xs, xm, Mx, scale = 1, N; in FormInitialSolution() local
228 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xm, NULL, NULL)); in FormInitialSolution()
247 for (i = xs; i < xs + xm; i++) { in FormInitialSolution()
279 PetscInt Mx, i, xs, xm, cnt; in MyMonitor() local
296 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xm, NULL, NULL)); in MyMonitor()
315 xs = xx[0] / hx; in MyMonitor()
324 x = hx * xs; in MyMonitor()
[all …]
/petsc/src/snes/tutorials/
H A Dex3.c86 PetscInt its, N = 5, i, maxit, maxf, xs, xm; in main() local
239 PetscCall(DMDAGetCorners(ctx.da, &xs, NULL, NULL, &xm, NULL, NULL)); in main()
250 xp = ctx.h * xs; in main()
251 for (i = xs; i < xs + xm; i++) { in main()
353 PetscInt i, M, xs, xm; in FormFunction() local
381 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xm, NULL, NULL)); in FormFunction()
389 if (xs == 0) { /* left boundary */ in FormFunction()
391 xs++; in FormFunction()
394 if (xs + xm == M) { /* right boundary */ in FormFunction()
395 ff[xs + xm - 1] = xx[xs + xm - 1] - 1.0; in FormFunction()
[all …]
H A Dex3k.kokkos.cxx52 PetscInt i, M, xs, xm; in CpuFunction() local
62 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xm, NULL, NULL)); in CpuFunction()
65 if (xs == 0) { /* left boundary */ in CpuFunction()
67 xs++; in CpuFunction()
70 if (xs + xm == M) { /* right boundary */ in CpuFunction()
71 R[xs + xm - 1] = X[xs + xm - 1] - 1.0; in CpuFunction()
75 for (i = xs; i < xs + xm; i++) R[i] = d * (X[i - 1] - 2.0 * X[i] + X[i + 1]) + X[i] * X[i] - F[i]; in CpuFunction()
155 PetscInt i, j[3], M, xs, xm; in FormJacobian() local
163 PetscCall(DMDAGetCorners(da, &xs, NULL, NULL, &xm, NULL, NULL)); in FormJacobian()
175 if (xs == 0) { /* left boundary */ in FormJacobian()
[all …]
H A Dex55k.kokkos.cxx41 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info->my; in FormFunctionLocalVec() local
61 …, MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>(exec, {ys, xs}, {ys + ym, xs + xm}), KOKK… 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()
176 …PetscInt xs = info->xs, ys = info->ys, xm = info->xm, ym = info->ym, mx = info->mx, my = info-… in FormJacobianLocalVec() local
190 hx = xm > 1 ? PetscRealPart(coords[ys][xs + 1].x) - PetscRealPart(coords[ys][xs].x) : 1.0; in FormJacobianLocalVec()
191 hy = ym > 1 ? PetscRealPart(coords[ys + 1][xs].y) - PetscRealPart(coords[ys][xs].y) : 1.0; in FormJacobianLocalVec()
209 for (i = xs; i < xs + xm; i++) { in FormJacobianLocalVec()
264 …, MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>(exec, {ys, xs}, {ys + ym, xs + xm}), KOKK… in FormJacobianLocalVec()
265 PetscInt p = ((j - ys) * xm + (i - xs)) * 5; in FormJacobianLocalVec()
/petsc/src/tao/bound/tutorials/
H A Dplate2f.F9066 PetscInt xs, xm
93 PetscCall(DMDAGetCorners(dm, xs, ys, PETSC_NULL_INTEGER, xm, ym, PETSC_NULL_INTEGER, ierr))
113 do i = xs, xs + xm - 1
131 xb = bottom_v(1 + i - xs + 1)
132 xrb = bottom_v(1 + i - xs + 2)
145 xt = top_v(1 + i - xs + 1)
146 xlt = top_v(1 + i - xs)
205 if (xs == 0) then ! left side
214 do i = xs, xs + xm - 1
215 d2 = (bottom_v(1 + i + 1 - xs) - bottom_v(1 + i - xs + 2))*rhx
[all …]
H A Dplate2.c202 PetscInt xs, xm, gxs, gxm, ys, ym, gys, gym; in FormFunctionGradient() local
214 PetscCall(DMDAGetCorners(user->dm, &xs, &ys, NULL, &xm, &ym, NULL)); in FormFunctionGradient()
234 for (i = xs; i < xs + xm; i++) { in FormFunctionGradient()
248 xb = bottom[i - xs + 1]; in FormFunctionGradient()
249 xrb = bottom[i - xs + 2]; in FormFunctionGradient()
262 xt = top[i - xs + 1]; in FormFunctionGradient()
263 xlt = top[i - xs]; in FormFunctionGradient()
317 if (xs == 0) { /* left side */ in FormFunctionGradient()
325 for (i = xs; i < xs + xm; i++) { in FormFunctionGradient()
326 d2 = (bottom[i + 1 - xs] - bottom[i - xs + 2]) * rhx; in FormFunctionGradient()
[all …]
/petsc/src/dm/tests/
H A Dex10k.kokkos.cxx30 PetscInt M = 5, N = 7, xm, ym, xs, ys; in main() local
52 PetscCall(DMDAGetCorners(da, &xs, &ys, 0, &xm, &ym, 0)); in main()
61 for (PetscInt i = xs; i < xs + xm; i++) { in main()
70 …"init 1", MDRangePolicy<Rank<3, Iterate::Right, Iterate::Right>>({ys, xs, 0}, {ys + ym, xs + xm, d… in main()
83 for (PetscInt i = xs; i < xs + xm; i++) { 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()
134 for (PetscInt i = xs; i < xs + xm; i++) { in main()
146 …"init 2", MDRangePolicy<Rank<2, Iterate::Right, Iterate::Right>>({ys, xs}, {ys + ym, xs + xm}), KO… in main()
161 for (PetscInt i = xs; i < xs + xm; i++) { in main()
[all …]
H A Dex2k.kokkos.cxx36 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()
45 …for (i = xs; i < xs + xm; i++) y1[k][j][i] = 6 * x1[k][j][i] - x1[k - 1][j][i] - x1[k][j - 1][i] -… 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()
70 …for (i = xs; i < xs + xm; i++) Y2(k, j, i) = 6 * X2(k, j, i) - X2(k - 1, j, i) - X2(k, j - 1, i) -… 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()
/petsc/src/ts/tutorials/multirate/
H A Dex8.c126 PetscInt i, j, k, dof, xs, xm, Mx; in FVSample_3WaySplit() local
133 PetscCall(DMDAGetCorners(da, &xs, 0, 0, &xm, 0, 0)); in FVSample_3WaySplit()
140 for (i = xs; i < xs + xm; i++) { in FVSample_3WaySplit()
223 PetscInt i, j, k, Mx, dof, xs, xm, sm = ctx->sm, mf = ctx->mf, fm = ctx->fm, ms = ctx->ms; in FVRHSFunction_3WaySplit() local
245 PetscCall(DMDAGetCorners(da, &xs, 0, 0, &xm, 0, 0)); in FVRHSFunction_3WaySplit()
248 for (i = xs - 2; i < 0; i++) { in FVRHSFunction_3WaySplit()
251 for (i = Mx; i < xs + xm + 2; i++) { in FVRHSFunction_3WaySplit()
252 for (j = 0; j < dof; j++) x[i * dof + j] = x[(xs + xm - 1) * dof + j]; in FVRHSFunction_3WaySplit()
255 for (i = xs - 1; i < xs + xm + 1; i++) { in FVRHSFunction_3WaySplit()
286 for (i = xs; i < xs + xm + 1; i++) { in FVRHSFunction_3WaySplit()
[all …]
H A Dex6.c133 PetscInt i, j, k, dof, xs, xm, Mx; in FVSample_2WaySplit() local
140 PetscCall(DMDAGetCorners(da, &xs, 0, 0, &xm, 0, 0)); in FVSample_2WaySplit()
145 for (i = xs; i < xs + xm; i++) { in FVSample_2WaySplit()
208 PetscInt i, j, k, Mx, dof, xs, xm, sf = ctx->sf, fs = ctx->fs; in FVRHSFunction_2WaySplit() local
229 PetscCall(DMDAGetCorners(da, &xs, 0, 0, &xm, 0, 0)); in FVRHSFunction_2WaySplit()
232 for (i = xs - 2; i < 0; i++) { in FVRHSFunction_2WaySplit()
235 for (i = Mx; i < xs + xm + 2; i++) { in FVRHSFunction_2WaySplit()
236 for (j = 0; j < dof; j++) x[i * dof + j] = x[(xs + xm - 1) * dof + j]; in FVRHSFunction_2WaySplit()
239 for (i = xs - 1; i < xs + xm + 1; i++) { in FVRHSFunction_2WaySplit()
269 for (i = xs; i < xs + xm + 1; i++) { in FVRHSFunction_2WaySplit()
[all …]

1234567