Home
last modified time | relevance | path

Searched refs:hx (Results 1 – 25 of 108) sorted by relevance

12345

/petsc/src/dm/impls/stag/tests/
H A Dex30.c69 PetscReal hx, hy, hz; in CreateMat() local
81 hx = 1.0 / N[0]; in CreateMat()
164 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 1.0 / (hz * hz); in CreateMat()
177 valA[2] = 1.0 / (hx * hx); in CreateMat()
183 valA[3] = 1.0 / (hx * hx); in CreateMat()
196 valA[5] = 1.0 / hx; in CreateMat()
202 valA[6] = -1.0 / hx; in CreateMat()
210 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 1.0 / (hz * hz); in CreateMat()
223 valA[2] = 1.0 / (hx * hx); in CreateMat()
229 valA[3] = 1.0 / (hx * hx); in CreateMat()
[all …]
H A Dex15.c245 PetscReal hx, hy, dv; in CreateSystem2d() local
254 hx = 1.0 / N[0]; in CreateSystem2d()
256 dv = hx * hy; in CreateSystem2d()
315 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem2d()
331 valA[3] = dv * 1.0 / (hx * hx); in CreateSystem2d()
349 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem2d()
364 valA[3] = dv * 1.0 / (hx * hx); in CreateSystem2d()
382 valA[0] = -dv * 2.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem2d()
397 valA[3] = dv * 1.0 / (hx * hx); in CreateSystem2d()
402 valA[4] = dv * 1.0 / (hx * hx); in CreateSystem2d()
[all …]
H A Dex18.c256 PetscReal hx, hy, dv, bogusScale; in CreateSystem() local
281 hx = 1.0 / N[0]; in CreateSystem()
283 dv = hx * hy; in CreateSystem()
355 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem()
371 valA[3] = dv * 1.0 / (hx * hx); in CreateSystem()
392 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem()
407 valA[3] = dv * 1.0 / (hx * hx); in CreateSystem()
428 valA[0] = -dv * 2.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem()
443 valA[3] = dv * 1.0 / (hx * hx); in CreateSystem()
448 valA[4] = dv * 1.0 / (hx * hx); in CreateSystem()
[all …]
/petsc/src/ts/tutorials/
H A Dex34.c52 const PetscScalar hx = i + 1 == info.xs + info.xm ? x[i] - x[i - 1] : x[i + 1] - x[i]; in FormRHSFunction() local
54 f[i].u = hx * (u[i].v); in FormRHSFunction()
55 …f[i].v = -hx * (PetscSqr(user->gammaTilde) * u[i].u + (PetscSqr(user->gamma) / user->xi) * (u[i].… in FormRHSFunction()
56 f[i].th = -hx * (u[i].v + 1) * (u[i].th + (1 + user->epsilon) * PetscLogScalar(u[i].v + 1)); in FormRHSFunction()
88 const PetscScalar hx = x[i + 1] - x[i]; in FormIFunction() local
89 f[i].u = hx * udot[i].u; in FormIFunction()
90 f[i].v = hx * udot[i].v - PetscSqr(user->c) * (u[i + 1].u - u[i].u) / hx; in FormIFunction()
91 f[i].th = hx * udot[i].th; in FormIFunction()
93 const PetscScalar hx = x[i] - x[i - 1]; in FormIFunction() local
94 f[i].u = hx * udot[i].u; in FormIFunction()
[all …]
H A Dex25.c45 PetscReal ftime, hx, dt; in main() local
106 hx = 1.0 / (PetscReal)(mx / 2 - 1); in main()
107 dt = 0.4 * PetscSqr(hx) / user.alpha; /* Diffusive stability limit */ in main()
142 PetscReal hx; in FormIFunction() local
148 hx = 1.0 / (PetscReal)(info.mx - 1); in FormIFunction()
168 f[i].u = hx * (x[i].u - user->uleft); in FormIFunction()
169 f[i].v = hx * (x[i].v - user->vleft); in FormIFunction()
171 f[i].u = hx * (x[i].u - user->uright); in FormIFunction()
172 f[i].v = hx * (x[i].v - user->vright); in FormIFunction()
174 f[i].u = hx * xdot[i].u - user->alpha * (x[i - 1].u - 2. * x[i].u + x[i + 1].u) / hx; in FormIFunction()
[all …]
H A Dex8.py15 mx = da.getSizes(); mx = mx[0]; hx = 1.0/mx
24 if xs == 0: f[0] = xx[0]/hx; lxs+=1;
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])
30 mx = da.getSizes(); mx = mx[0]; hx = 1.0/mx
38 if xs == 0: P.setValues([0],[0],1.0/hx); lxs+=1
39 if lxe == mx: P.setValues([mx-1],[mx-1],1.0/hx); lxe-=1
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 Dex35.cxx97 PetscReal hx, dt, ftime; in main() local
143 hx = 1.0 / (PetscReal)(mx / 2 - 1); in main()
144 dt = 0.4 * PetscSqr(hx) / user->alpha; /* Diffusive stability limit */ in main()
199 PetscReal hx; in FormIJacobian() local
211 hx = 1.0 / user->n; in FormIJacobian()
226 {hx, 0 }, in FormIJacobian()
227 {0, hx} in FormIJacobian()
232 …const PetscScalar dxxL = -user->alpha / hx, dxx0 = 2. * user->alpha / hx, dxxR = -user->alpha / hx; in FormIJacobian()
234 {{dxxL, 0}, {a * hx + dxx0, 0}, {dxxR, 0}}, in FormIJacobian()
235 {{0, dxxL}, {0, a * hx + dxx0}, {0, dxxR}} in FormIJacobian()
[all …]
/petsc/src/ts/tests/
H A Dex25.c63 PetscReal ftime, hx, dt, xmax, xmin; in Brusselator() local
124 hx = 1.0 / (PetscReal)(mx / 2 - 1); in Brusselator()
125 dt = 0.4 * PetscSqr(hx) / user.alpha; /* Diffusive stability limit */ in Brusselator()
164 PetscReal hx; in FormIFunction() local
170 hx = 1.0 / (PetscReal)(info.mx - 1); in FormIFunction()
190 f[i].u = hx * (x[i].u - user->uleft); in FormIFunction()
191 f[i].v = hx * (x[i].v - user->vleft); in FormIFunction()
193 f[i].u = hx * (x[i].u - user->uright); in FormIFunction()
194 f[i].v = hx * (x[i].v - user->vright); in FormIFunction()
196 f[i].u = hx * xdot[i].u - user->alpha * (x[i - 1].u - 2. * x[i].u + x[i + 1].u) / hx; in FormIFunction()
[all …]
H A Dex21.c37 PetscReal lambda, hx, hy; in FormIFunctionLocal() local
42 hx = 1.0 / (PetscReal)(info->mx - 1); in FormIFunctionLocal()
62 uxx = (uw - 2.0 * u + ue) / (hx * hx); in FormIFunctionLocal()
78 PetscScalar v[5], lambda, hx, hy; in FormIJacobianLocal() local
82 hx = 1.0 / (PetscReal)(info->mx - 1); in FormIJacobianLocal()
103 v[k] = -1.0 / (hx * hx); in FormIJacobianLocal()
108 v[k] = shift + 2.0 / (hx * hx) + 2.0 / (hy * hy) - lambda * PetscExpScalar(x[j][i]); in FormIJacobianLocal()
113 v[k] = -1.0 / (hx * hx); in FormIJacobianLocal()
/petsc/src/dm/impls/stag/tutorials/
H A Dex3.c172 PetscReal hx, hy, hz; in CreateSystem() local
185 hx = 1.0 / N[0]; in CreateSystem()
282 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 1.0 / (hz * hz); in CreateSystem()
295 valA[2] = 1.0 / (hx * hx); in CreateSystem()
301 valA[3] = 1.0 / (hx * hx); in CreateSystem()
314 valA[5] = 1.0 / hx; in CreateSystem()
320 valA[6] = -1.0 / hx; in CreateSystem()
328 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 1.0 / (hz * hz); in CreateSystem()
341 valA[2] = 1.0 / (hx * hx); in CreateSystem()
347 valA[3] = 1.0 / (hx * hx); in CreateSystem()
[all …]
H A Dex2.c152 PetscReal hx, hy; in CreateSystem() local
171 hx = 1.0 / N[0]; in CreateSystem()
237 valA[0] = -1.0 / (hx * hx) - 2.0 / (hy * hy); in CreateSystem()
253 valA[3] = 1.0 / (hx * hx); in CreateSystem()
271 valA[0] = -1.0 / (hx * hx) - 2.0 / (hy * hy); in CreateSystem()
286 valA[3] = 1.0 / (hx * hx); in CreateSystem()
304 valA[0] = -2.0 / (hx * hx) - 2.0 / (hy * hy); in CreateSystem()
319 valA[3] = 1.0 / (hx * hx); in CreateSystem()
324 valA[4] = 1.0 / (hx * hx); in CreateSystem()
369 valA[0] = -2.0 / (hx * hx) - 1.0 / (hy * hy); in CreateSystem()
[all …]
/petsc/src/ts/tutorials/multirate/
H A Dex5.c191 PetscReal hx, cfl_idt = 0; in FVRHSFunctionslow() local
200 hx = (ctx->xmax - ctx->xmin) / Mx; in FVRHSFunctionslow()
223 …eristic)(ctx->physics.user, dof, &x[i * dof], ctx->R, ctx->Rinv, ctx->speeds, ctx->xmin + hx * i)); in FVRHSFunctionslow()
239 info.hx = hx; in FVRHSFunctionslow()
241 for (j = 0; j < dof; j++) ctx->cslope[j] /= hx; /* rescale to a slope */ in FVRHSFunctionslow()
257 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hx / 2; in FVRHSFunctionslow()
258 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hx / 2; in FVRHSFunctionslow()
260 …riemann)(ctx->physics.user, dof, uL, uR, ctx->flux, &maxspeed, ctx->xmin + hx * i, ctx->xmin, ctx-… in FVRHSFunctionslow()
261 … cfl_idt = PetscMax(cfl_idt, PetscAbsScalar(maxspeed / hx)); /* Max allowable value of 1/Delta t */ in FVRHSFunctionslow()
263 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hx; in FVRHSFunctionslow()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/bratu2d/
H A Dbratu2df90.f9012 real(kind=8) :: hx, hy local
22 hx = 1.0/(m-1) ! x grid spacing
26 (2*u - uE - uW) * (hy/hx) &
27 + (2*u - uN - uS) * (hx/hy) &
28 - alpha * exp(u) * (hx*hy)
H A Dbratu2dnpy.py14 hx = 1.0/(nx-1) # x grid spacing
18 (2*u - uE - uW) * (hy/hx) \
19 + (2*u - uN - uS) * (hx/hy) \
20 - alpha * exp(u) * (hx*hy)
/petsc/src/tao/bound/tutorials/
H A Djbearing2.c177 PetscReal hx, hy, ehxhy; in ComputeB() local
184 hx = two * pi / (nx + 1.0); in ComputeB()
186 ehxhy = ecc * hx * hy; in ComputeB()
197 temp = PetscSinScalar((i + 1) * hx); in ComputeB()
214 PetscReal hx, hy, hxhy, hxhx, hyhy; in FormFunctionGradient() local
225 hx = two * pi / (nx + 1.0); in FormFunctionGradient()
227 hxhy = hx * hy; in FormFunctionGradient()
228 hxhx = one / (hx * hx); in FormFunctionGradient()
247 xi = (i + 1) * hx; in FormFunctionGradient()
248 trule1 = hxhy * (p(xi, ecc) + p(xi + hx, ecc) + p(xi, ecc)) / six; /* L(i,j) */ in FormFunctionGradient()
[all …]
/petsc/src/ts/tutorials/phasefield/
H A Dheat.c154 PetscReal hx, sx; in FormFunction() local
164 hx = 1.0 / (PetscReal)Mx; in FormFunction()
165 sx = 1.0 / (hx * hx); in FormFunction()
210 PetscReal hx, x, r; in FormInitialSolution() local
218 hx = 1.0 / (PetscReal)Mx; in FormInitialSolution()
248 x = i * hx; in FormInitialSolution()
280 PetscReal x, y, hx, pause, sx, len, max, xx[2], yy[2]; in MyMonitor() local
297 hx = 1.0 / (PetscReal)Mx; in MyMonitor()
298 sx = 1.0 / (hx * hx); in MyMonitor()
315 xs = xx[0] / hx; in MyMonitor()
[all …]
/petsc/src/snes/tutorials/
H A Dex46.c125 PetscReal D, K, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
131 hx = 1.0 / (PetscReal)(info->mx - 1); in FormFunctionLocal()
133 hxdhy = hx / hy; in FormFunctionLocal()
134 hydhx = hy / hx; in FormFunctionLocal()
146 ux = (x[j][i + 1] - x[j][i]) / hx; in FormFunctionLocal()
150 …- (K * funcA(x[j][i], user) * PetscSqrtScalar(ux * ux + uy * uy) + funcU(&coords[j][i])) * hx * hy; in FormFunctionLocal()
166 PetscScalar D, K, A, v[5], hx, hy, hxdhy, hydhx, ux, uy; in FormJacobianLocal() local
173 hx = 1.0 / (PetscReal)(info->mx - 1); in FormJacobianLocal()
175 hxdhy = hx / hy; in FormJacobianLocal()
176 hydhx = hy / hx; in FormJacobianLocal()
[all …]
H A Dex47cu.cu79 PetscScalar *xx, *ff, hx; in ComputeFunction() local
89 hx = 1.0 / (PetscReal)(Mx - 1); in ComputeFunction()
117 thrust::constant_iterator<PetscScalar>(hx))), in ComputeFunction()
126 thrust::constant_iterator<PetscScalar>(hx))), in ComputeFunction()
141 if (i == 0 || i == Mx - 1) ff[i] = xx[i] / hx; in ComputeFunction()
142 else ff[i] = (2.0 * xx[i] - xx[i - 1] - xx[i + 1]) / hx - hx * PetscExpScalar(xx[i]); in ComputeFunction()
154 PetscScalar hx, *xx; in ComputeJacobian() local
159 hx = 1.0 / (PetscReal)(Mx - 1); in ComputeJacobian()
168 PetscCall(MatSetValue(J, i, i, 1.0 / hx, INSERT_VALUES)); in ComputeJacobian()
170 PetscCall(MatSetValue(J, i, i - 1, -1.0 / hx, INSERT_VALUES)); in ComputeJacobian()
[all …]
H A Dex55k.kokkos.cxx40 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocalVec() local
50 hx = 1.0 / (PetscReal)(info->mx - 1); in FormFunctionLocalVec()
52 hxdhy = hx / hy; in FormFunctionLocalVec()
53 hydhx = hy / hx; in FormFunctionLocalVec()
66 c.x = i * hx; in FormFunctionLocalVec()
79 c.x = (i - 1) * hx; in FormFunctionLocalVec()
84 c.x = (i + 1) * hx; in FormFunctionLocalVec()
89 c.x = i * hx; in FormFunctionLocalVec()
94 c.x = i * hx; in FormFunctionLocalVec()
102 c.x = i * hx; in FormFunctionLocalVec()
[all …]
H A Dex25.c66 PetscScalar hx, hy; in FormFunctionLocal() local
71 hx = 1.0 / (PetscReal)(info->mx - 1); in FormFunctionLocal()
78 … f[j][i] = t[j][i] - (1.0 - (2.0 * hx * (PetscReal)i - 1.0) * (2.0 * hx * (PetscReal)i - 1.0)); in FormFunctionLocal()
82 gradright = (t[j][i + 1] - t[j][i]) / hx; in FormFunctionLocal()
83 gradleft = (t[j][i] - t[j][i - 1]) / hx; in FormFunctionLocal()
85 gradx = .5 * (t[j][i + 1] - t[j][i - 1]) / hx; in FormFunctionLocal()
94 …f[j][i] = (coeffup * gradup - coeffdown * graddown) * hx + (coeffright * gradright - coeffleft * g… in FormFunctionLocal()
H A Dex4.c118 PetscScalar hx, hy; in FormObjectiveLocal() local
124 hx = 1.0 / (mx + 1); in FormObjectiveLocal()
126 area = 0.5 * hx * hy; in FormObjectiveLocal()
153 d1 /= hx; in FormObjectiveLocal()
154 d2 /= hx; in FormObjectiveLocal()
169 d2 = (user->left[j + 1] - x[j][0]) / hx; in FormObjectiveLocal()
175 d2 = (user->bottom[i + 1] - user->bottom[i + 2]) / hx; in FormObjectiveLocal()
182 d1 = (x[j][mx - 1] - user->right[j + 1]) / hx; in FormObjectiveLocal()
190 d4 = (user->top[i + 1] - user->top[i]) / hx; in FormObjectiveLocal()
196 d2 = (user->bottom[0] - user->bottom[1]) / hx; in FormObjectiveLocal()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/wrap-swig/
H A DBratu3D.c25 PetscReal lambda,temp1,hx,hy,hz,tempk,tempj; in FormInitGuess() local
31 hx = 1.0/(PetscReal)(Mx-1); in FormInitGuess()
69 x[k][j][i] = temp1*sqrt(PetscMin((PetscReal)(PetscMin(i,Mx-i-1))*hx,tempj)); in FormInitGuess()
85 PetscReal two = 2.0,lambda,hx,hy,hz,hxhzdhy,hyhzdhx,hxhydhz,sc; in FormFunction() local
92 hx = 1.0/(PetscReal)(Mx-1); in FormFunction()
95 sc = hx*hy*hz*lambda; in FormFunction()
96 hxhzdhy = hx*hz/hy; in FormFunction()
97 hyhzdhx = hy*hz/hx; in FormFunction()
98 hxhydhz = hx*hy/hz; in FormFunction()
162 PetscReal lambda,hx,hy,hz,hxhzdhy,hyhzdhx,hxhydhz,sc; in FormJacobian() local
[all …]
/petsc/src/binding/petsc4py/demo/legacy/wrap-cython/
H A DBratu3Dimpl.c25 PetscReal lambda,temp1,hx,hy,hz,tempk,tempj; in FormInitGuess() local
32 hx = 1.0/(PetscReal)(Mx-1); in FormInitGuess()
70 x[k][j][i] = temp1*sqrt(PetscMin((PetscReal)(PetscMin(i,Mx-i-1))*hx,tempj)); in FormInitGuess()
86 PetscReal two = 2.0,lambda,hx,hy,hz,hxhzdhy,hyhzdhx,hxhydhz,sc; in FormFunction() local
93 hx = 1.0/(PetscReal)(Mx-1); in FormFunction()
96 sc = hx*hy*hz*lambda; in FormFunction()
97 hxhzdhy = hx*hz/hy; in FormFunction()
98 hyhzdhx = hy*hz/hx; in FormFunction()
99 hxhydhz = hx*hy/hz; in FormFunction()
166 PetscReal lambda,hx,hy,hz,hxhzdhy,hyhzdhx,hxhydhz,sc; in FormJacobian() local
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex44f.F9050 PetscScalar hx
54 hx = 1.0_PETSC_REAL_KIND/(mx - 1)
57 xx(i) = i*hx
71 PetscScalar hx, one
76 hx = 1.0_PETSC_REAL_KIND/(mx - 1)
81 PetscCall(MatSetValue(J, i, i - 1, -hx, INSERT_VALUES, ierr))
82 PetscCall(MatSetValue(J, i, i + 1, -hx, INSERT_VALUES, ierr))
83 PetscCall(MatSetValue(J, i, i, 2*hx, INSERT_VALUES, ierr))
/petsc/src/binding/petsc4py/demo/python_types/
H A Dmat.py43 hx, hy = (1.0 / m for m in [mx, my])
68 u_xx = (-u_e + 2 * u - u_w) * hy / hx
69 u_yy = (-u_n + 2 * u - u_s) * hx / hy
75 hx, hy = (1.0 / m for m in [mx, my])
83 d[i, j] = 2 * hy / hx + 2 * hx / hy
89 hx, hy = (1.0 / m for m in [mx, my])
93 b[i, j] = 1 * hx * hy

12345