Home
last modified time | relevance | path

Searched refs:sx (Results 1 – 24 of 24) sorted by relevance

/petsc/src/ts/tutorials/phasefield/
H A Dbiharmonic.c201 PetscReal hx, sx; in FormFunction() local
213 sx = 1.0 / (hx * hx); in FormFunction()
240 c = (1. - x[i] * x[i]) * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
241 r = (1. - x[i + 1] * x[i + 1]) * (x[i] + x[i + 2] - 2.0 * x[i + 1]) * sx; in FormFunction()
242 l = (1. - x[i - 1] * x[i - 1]) * (x[i - 2] + x[i] - 2.0 * x[i - 1]) * sx; in FormFunction()
244 c = (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
245 r = (x[i] + x[i + 2] - 2.0 * x[i + 1]) * sx; in FormFunction()
246 l = (x[i - 2] + x[i] - 2.0 * x[i - 1]) * sx; in FormFunction()
248 f[i] = -ctx->kappa * (l + r - 2.0 * c) * sx; in FormFunction()
252 …1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (3. * x[i] * x[i] - 1.) * (x[i - 1] + x[i + 1] - 2.0… in FormFunction()
[all …]
H A Dbiharmonic3.c186 PetscReal hx, sx; in FormFunction() local
199 sx = 1.0 / (hx * hx); in FormFunction()
228 f[i].w = x[i].w + ctx->kappa * (x[i - 1].u + x[i + 1].u - 2.0 * x[i].u) * sx; in FormFunction()
246 f[i].u = xdot[i].u - (x[i - 1].w + x[i + 1].w - 2.0 * x[i].w) * sx; in FormFunction()
253 …u - x[i - 1].u) * .25 * sx - (ctx->theta - ctx->theta_c * (1 - x[i].u * x[i].u)) * (x[i + 1].u + x… in FormFunction()
273 PetscReal hx, xx, r, sx; in FormInitialSolution() local
280 sx = 1.0 / (hx * hx); in FormInitialSolution()
305 for (i = xs; i < xs + xm; i++) x[i].w = -kappa * (x[i - 1].u + x[i + 1].u - 2.0 * x[i].u) * sx; in FormInitialSolution()
H A Dbiharmonic2.c187 PetscReal hx, sx; in FormFunction() local
199 sx = 1.0 / (hx * hx); in FormFunction()
228 f[i].w = x[i].w + ctx->kappa * (x[i - 1].u + x[i + 1].u - 2.0 * x[i].u) * sx; in FormFunction()
244 f[i].u = xdot[i].u - (x[i - 1].w + x[i + 1].w - 2.0 * x[i].w) * sx; in FormFunction()
263 PetscReal hx, xx, r, sx; in FormInitialSolution() local
270 sx = 1.0 / (hx * hx); in FormInitialSolution()
295 for (i = xs; i < xs + xm; i++) x[i].w = -kappa * (x[i - 1].u + x[i + 1].u - 2.0 * x[i].u) * sx; in FormInitialSolution()
H A Dheat.c154 PetscReal hx, sx; in FormFunction() local
165 sx = 1.0 / (hx * hx); in FormFunction()
191 f[i] = ctx->kappa * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
280 PetscReal x, y, hx, pause, sx, len, max, xx[2], yy[2]; in MyMonitor() local
298 sx = 1.0 / (hx * hx); in MyMonitor()
327 …[0] = PetscRealPart(.25 * ctx->kappa * (u[i - 1] - u[i + 1]) * (u[i - 1] - u[i + 1]) * sx); in MyMonitor()
347 yy[0] = PetscRealPart(ctx->kappa * (u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx); in MyMonitor()
387 len = .5 * PetscRealPart(ctx->kappa * (u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx) / max; in MyMonitor()
/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dreaction_diffusion.c32 PetscReal hx, hy, sx, sy; in RHSFunction() local
42 sx = 1.0 / (hx * hx); in RHSFunction()
72 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in RHSFunction()
75 vxx = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx; in RHSFunction()
97 PetscReal hx, hy, sx, sy; in RHSJacobian() local
110 sx = 1.0 / (hx * hx); in RHSJacobian()
171 entries[2] = appctx->D1 * sx; in RHSJacobian()
174 entries[3] = appctx->D1 * sx; in RHSJacobian()
177 entries[4] = -2.0 * appctx->D1 * (sx + sy) - vc * vc - appctx->gamma; in RHSJacobian()
191 entries[2] = appctx->D2 * sx; in RHSJacobian()
[all …]
H A Dex5adj_mf.c64 PetscReal hx, hy, sx, sy; in MyRHSMatMultTranspose() local
76 sx = 1.0 / (hx * hx); in MyRHSMatMultTranspose()
100 uxx = (-2.0 * ucb + x[j][i - 1].u + x[j][i + 1].u) * sx; in MyRHSMatMultTranspose()
104 vxx = (-2.0 * vcb + x[j][i - 1].v + x[j][i + 1].v) * sx; in MyRHSMatMultTranspose()
123 PetscReal hx, hy, sx, sy; in MyIMatMultTranspose() local
135 sx = 1.0 / (hx * hx); in MyIMatMultTranspose()
159 uxx = (-2.0 * ucb + x[j][i - 1].u + x[j][i + 1].u) * sx; in MyIMatMultTranspose()
163 vxx = (-2.0 * vcb + x[j][i - 1].v + x[j][i + 1].v) * sx; in MyIMatMultTranspose()
184 PetscReal hx, hy, sx, sy; in MyIMatMult() local
196 sx = 1.0 / (hx * hx); in MyIMatMult()
[all …]
H A Dex4.c140 PetscReal hx, sx; in IFunction() local
151 sx = 1.0 / (hx * hx); in IFunction()
191 rhoxx = (-2.0 * rho + u[i - 1].rho + u[i + 1].rho) * sx; in IFunction()
193 cxx = (-2.0 * c + u[i - 1].c + u[i + 1].c) * sx; in IFunction()
200 …x = appctx->kappa * ((u[i + 1].c - u[i].c) * u[i + 1].rho - (u[i].c - u[i - 1].c) * u[i].rho) * sx; in IFunction()
/petsc/src/ts/tutorials/
H A Dex17.c161 PetscReal hx, sx; in FormIFunction() local
171 sx = 1.0 / (hx * hx); in FormIFunction()
194 else f[i] = udot[i] + (2. * u[i] - u[i - 1] - u[i + 1]) * sx; in FormIFunction()
198 else f[i] = udot[i] + (2. * u[i] - u[i - 1] - u[i + 1]) * sx; in FormIFunction()
217 PetscReal hx, sx; in FormIJacobian() local
229 sx = 1.0 / (hx * hx); in FormIJacobian()
248 vals[nc++] = -1.0 * sx; in FormIJacobian()
250 vals[nc++] = 2.0 * sx + a; in FormIJacobian()
252 vals[nc++] = -1.0 * sx; in FormIJacobian()
H A Dex13.c121 PetscReal two = 2.0, hx, hy, sx, sy; in RHSFunction() local
131 sx = 1.0 / (hx * hx); in RHSFunction()
159 uxx = (-two * u + uarray[j][i - 1] + uarray[j][i + 1]) * sx; in RHSFunction()
193 PetscReal hx, hy, sx, sy; in RHSJacobian() local
199 sx = 1.0 / (hx * hx); in RHSJacobian()
216 val[nc++] = sx; in RHSJacobian()
219 val[nc++] = sx; in RHSJacobian()
228 val[nc++] = -2 * sx - 2 * sy; in RHSJacobian()
H A Dex15.c160 PetscReal hx, hy, sx, sy; in FormIFunction() local
169 sx = 1.0 / (hx * hx); in FormIFunction()
227 f[j][i] = udot[j][i] - (uxx * sx + uyy * sy); in FormIFunction()
252 PetscScalar vals[5], hx, hy, sx, sy; in FormIJacobian() local
259 sx = 1.0 / (hx * hx); in FormIJacobian()
307 vals[nc++] = -sx; in FormIJacobian()
310 vals[nc++] = 2.0 * (sx + sy) + a; in FormIJacobian()
313 vals[nc++] = -sx; in FormIJacobian()
H A Dex7.c128 PetscReal two = 2.0, hx, hy, sx, sy; in FormFunction() local
138 sx = 1.0 / (hx * hx); in FormFunction()
172 uxx = (two * u - x[j][i - 1] - x[j][i + 1]) * sx; in FormFunction()
175 … 1.0) * ((x[j][i + 1] - x[j][i - 1]) * (x[j][i + 1] - x[j][i - 1]) * .25 * sx + (x[j + 1][i] - x[j… in FormFunction()
H A Dex12.c135 PetscReal hx, hy, /*hxdhy,hydhx,*/ sx, sy; in FormFunction() local
144 sx = 1.0 / (hx * hx); in FormFunction()
182 uxx = (-2.0 * u + x[j][i - 1][0] + x[j][i + 1][0]) * sx; in FormFunction()
/petsc/src/ts/tutorials/autodiff/
H A Dadr_ex5adj.cxx364 PetscReal hx, hy, sx, sy; in IFunctionLocalPassive() local
369 sx = 1.0 / (hx * hx); in IFunctionLocalPassive()
383 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in IFunctionLocalPassive()
386 vxx = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx; in IFunctionLocalPassive()
404 PetscReal hx, hy, sx, sy; in IFunctionActive() local
414 sx = 1.0 / (hx * hx); in IFunctionActive()
479 uxx = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx; in IFunctionActive()
482 vxx = (-2.0 * vc + u_a[j][i - 1].v + u_a[j][i + 1].v) * sx; in IFunctionActive()
531 PetscReal hx, hy, sx, sy; in RHSFunctionPassive() local
540 sx = 1.0 / (hx * hx); in RHSFunctionPassive()
[all …]
H A Dadr_ex5adj_mf.cxx252 PetscReal hx, hy, sx, sy; in IFunctionLocalPassive() local
257 sx = 1.0 / (hx * hx); in IFunctionLocalPassive()
271 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in IFunctionLocalPassive()
274 vxx = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx; in IFunctionLocalPassive()
292 PetscReal hx, hy, sx, sy; in IFunctionActive() local
302 sx = 1.0 / (hx * hx); in IFunctionActive()
367 uxx = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx; in IFunctionActive()
370 vxx = (-2.0 * vc + u_a[j][i - 1].v + u_a[j][i + 1].v) * sx; in IFunctionActive()
/petsc/src/ts/tests/
H A Dex24.c21 PetscReal hx, hy, sx, sy; in RHSJacobian() local
34 sx = 1.0 / (hx * hx); in RHSJacobian()
95 entries[2] = appctx->D1 * sx; in RHSJacobian()
98 entries[3] = appctx->D1 * sx; in RHSJacobian()
101 entries[4] = -2.0 * appctx->D1 * (sx + sy) - vc * vc - appctx->gamma; in RHSJacobian()
114 entries[2] = appctx->D2 * sx; in RHSJacobian()
116 entries[3] = appctx->D2 * sx; in RHSJacobian()
118 entries[4] = -2.0 * appctx->D2 * (sx + sy) + 2.0 * uc * vc - appctx->gamma - appctx->kappa; in RHSJacobian()
/petsc/src/ksp/ksp/tests/
H A Dex51.c12 PetscScalar value[3], alpha, beta, sx; in main() local
24 sx = PetscSinReal(n * PETSC_PI / 2 / (n + 1)); in main()
25 alpha = 4.0 * sx * sx; /* alpha is the largest eigenvalue of the matrix */ in main()
/petsc/src/dm/tests/
H A Dex36.c58 PetscInt sx, nx, sy, ny, sz, nz, dim; in DAApplyConformalMapping() local
79 PetscCall(DMDAGetCorners(da, &sx, &sy, &sz, &nx, &ny, &nz)); in DAApplyConformalMapping()
198 PetscInt sx, nx, sy, ny, sz, nz; in DAApplyTrilinearMapping() local
210 PetscCall(DMDAGetCorners(da, &sx, &sy, &sz, &nx, &ny, &nz)); in DAApplyTrilinearMapping()
212 for (i = sx; i < sx + nx; i++) { in DAApplyTrilinearMapping()
253 PetscInt sx, nx, sy, ny; in DADefineXLinearField2D() local
266 PetscCall(DMDAGetCorners(da, &sx, &sy, 0, &nx, &ny, 0)); in DADefineXLinearField2D()
268 for (i = sx; i < sx + nx; i++) { in DADefineXLinearField2D()
280 PetscInt sx, nx, sy, ny, sz, nz; in DADefineXLinearField3D() local
293 PetscCall(DMDAGetCorners(da, &sx, &sy, &sz, &nx, &ny, &nz)); in DADefineXLinearField3D()
[all …]
/petsc/src/dm/impls/da/
H A Dgrglvis.c38 PetscInt sx, sy, sz, ien, jen, ken; in DMDAGetNumElementsGhosted() local
45 PetscCall(DMDAGetCorners(da, &sx, &sy, &sz, &ien, &jen, &ken)); in DMDAGetNumElementsGhosted()
51 if (!sx) ien--; in DMDAGetNumElementsGhosted()
58 if (sx + ien == M) ien--; in DMDAGetNumElementsGhosted()
97 PetscInt sx, sy, sz, gsx, gsy, gsz, ist, jst, kst, gm, gn, gp; in DMDASampleGLVisFields_Private() local
108 PetscCall(DMDAGetCorners(da, &sx, &sy, &sz, NULL, NULL, NULL)); in DMDASampleGLVisFields_Private()
114 ist = gsx != sx ? 1 : 0; in DMDASampleGLVisFields_Private()
343 PetscInt sx, sy, sz, ie, je, ke, ien, jen, ken, nel; in DMDAView_GLVis_ASCII() local
479 PetscCall(DMDAGetCorners(da, &sx, &sy, &sz, NULL, NULL, NULL)); in DMDAView_GLVis_ASCII()
486 ist = gsx != sx ? 1 : 0; in DMDAView_GLVis_ASCII()
/petsc/src/vec/vec/impls/nest/
H A Dvecnest.c839 static PetscErrorCode VecNestGetSubVec_Nest(Vec X, PetscInt idxm, Vec *sx) in VecNestGetSubVec_Nest() argument
843 PetscCall(VecNestGetSubVecs_Private(X, 1, &idxm, sx)); in VecNestGetSubVec_Nest()
863 PetscErrorCode VecNestGetSubVec(Vec X, PetscInt idxm, Vec *sx) in VecNestGetSubVec() argument
866 PetscUseMethod(X, "VecNestGetSubVec_C", (Vec, PetscInt, Vec *), (X, idxm, sx)); in VecNestGetSubVec()
870 static PetscErrorCode VecNestGetSubVecs_Nest(Vec X, PetscInt *N, Vec **sx) in VecNestGetSubVecs_Nest() argument
877 if (sx) *sx = b->v; in VecNestGetSubVecs_Nest()
903 PetscErrorCode VecNestGetSubVecs(Vec X, PetscInt *N, Vec *sx[]) in VecNestGetSubVecs() argument
906 PetscUseMethod(X, "VecNestGetSubVecs_C", (Vec, PetscInt *, Vec **), (X, N, sx)); in VecNestGetSubVecs()
978 static PetscErrorCode VecNestSetSubVec_Nest(Vec X, PetscInt idxm, Vec sx) in VecNestSetSubVec_Nest() argument
982 PetscCall(VecNestSetSubVec_Private(X, idxm, sx)); in VecNestSetSubVec_Nest()
[all …]
/petsc/include/
H A Dpetscblaslapack.h166 … BLASdot_(const PetscBLASInt *n, const PetscScalar *x, const PetscBLASInt *sx, const PetscScalar *… in BLASdot_() argument
170 if (*sx == 1 && *sy == 1) { in BLASdot_()
173 for (i = 0, j = 0, k = 0; i < *n; i++, j += *sx, k += *sy) sum += PetscConj(x[j]) * y[k]; in BLASdot_()
177 …BLASdotu_(const PetscBLASInt *n, const PetscScalar *x, const PetscBLASInt *sx, const PetscScalar *… in BLASdotu_() argument
181 if (*sx == 1 && *sy == 1) { in BLASdotu_()
184 for (i = 0, j = 0, k = 0; i < *n; i++, j += *sx, k += *sy) sum += x[j] * y[k]; in BLASdotu_()
/petsc/src/ksp/pc/impls/mpi/
H A Dpcmpi.c390 PetscScalar *b, *sx = NULL; in PCMPISolve() local
424 if (pc) PetscCall(VecGetArray(X, &sx)); in PCMPISolve()
425 const void *inaddr[2] = {(const void **)sb, (const void **)sx}; in PCMPISolve()
426 if (pc) PetscCall(VecRestoreArray(X, &sx)); in PCMPISolve()
451 if (pc) PetscCall(VecGetArray(X, &sx)); in PCMPISolve()
452 …PetscCallMPI(MPIU_Gatherv(x, n, MPIU_SCALAR, sx, pc ? km->sendcount : NULL, pc ? km->displ : NULL,… in PCMPISolve()
453 if (pc) PetscCall(VecRestoreArray(X, &sx)); in PCMPISolve()
/petsc/src/sys/classes/draw/impls/image/
H A Ddrawimage.h56 int dx = PetscAbs(x_2 - x_1), sx = (x_2 - x_1) >= 0 ? +1 : -1; in PetscImageDrawLine() local
65 x_1 += sx; in PetscImageDrawLine()
/petsc/src/vec/vec/utils/
H A Dvinv.c1511 const PetscScalar *sx, *tx; in VecDotNorm2() local
1515 PetscCall(VecGetArrayRead(s, &sx)); in VecDotNorm2()
1520 work[0] += sx[i] * txconj; in VecDotNorm2()
1524 PetscCall(VecRestoreArrayRead(s, &sx)); in VecDotNorm2()
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DVec.pyx3476 cdef PetscVec* sx=NULL
3477 CHKERR(VecNestGetSubVecs(self.vec, &N, &sx))
3481 pyvec.vec = sx[i]
3489 sx: Sequence[Vec],
3497 sx
3500 Indices of the component vectors, defaults to ``range(len(sx))``.
3507 if idxm is None: idxm = range(len(sx))
3508 else: assert len(idxm) == len(sx)
3515 for i from 0 <= i < N: csx[i] = (<Vec?>sx[i]).vec