Home
last modified time | relevance | path

Searched refs:hy (Results 1 – 25 of 80) sorted by relevance

1234

/petsc/src/dm/impls/stag/tests/
H A Dex30.c69 PetscReal hx, hy, hz; in CreateMat() local
82 hy = 1.0 / N[1]; in CreateMat()
164 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 1.0 / (hz * hz); in CreateMat()
171 valA[1] = 1.0 / (hy * hy); in CreateMat()
210 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 1.0 / (hz * hz); in CreateMat()
217 valA[1] = 1.0 / (hy * hy); in CreateMat()
256 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 2.0 / (hz * hz); in CreateMat()
263 valA[1] = 1.0 / (hy * hy); in CreateMat()
309 valA[0] = -2.0 / (hx * hx) + -2.0 / (hy * hy) - 1.0 / (hz * hz); in CreateMat()
315 valA[1] = 1.0 / (hy * hy); in CreateMat()
[all …]
H A Dex15.c245 PetscReal hx, hy, dv; in CreateSystem2d() local
255 hy = 1.0 / N[1]; in CreateSystem2d()
256 dv = hx * hy; in CreateSystem2d()
315 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem2d()
320 valA[1] = dv * 1.0 / (hy * hy); in CreateSystem2d()
325 valA[2] = dv * 1.0 / (hy * hy); in CreateSystem2d()
336 valA[4] = dv * 1.0 / hy; in CreateSystem2d()
341 valA[5] = -dv * 1.0 / hy; in CreateSystem2d()
349 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem2d()
354 valA[1] = dv * 1.0 / (hy * hy); in CreateSystem2d()
[all …]
H A Dex18.c256 PetscReal hx, hy, dv, bogusScale; in CreateSystem() local
282 hy = 1.0 / N[1]; in CreateSystem()
283 dv = hx * hy; in CreateSystem()
355 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem()
360 valA[1] = dv * 1.0 / (hy * hy); in CreateSystem()
365 valA[2] = dv * 1.0 / (hy * hy); in CreateSystem()
378 valA[4] = dv * 1.0 / hy; in CreateSystem()
383 valA[5] = -dv * 1.0 / hy; in CreateSystem()
392 valA[0] = -dv * 1.0 / (hx * hx) - dv * 2.0 / (hy * hy); in CreateSystem()
397 valA[1] = dv * 1.0 / (hy * hy); in CreateSystem()
[all …]
/petsc/src/dm/impls/stag/tutorials/
H A Dex3.c172 PetscReal hx, hy, hz; in CreateSystem() local
186 hy = 1.0 / N[1]; in CreateSystem()
282 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 1.0 / (hz * hz); in CreateSystem()
289 valA[1] = 1.0 / (hy * hy); in CreateSystem()
328 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 1.0 / (hz * hz); in CreateSystem()
335 valA[1] = 1.0 / (hy * hy); in CreateSystem()
374 valA[0] = -2.0 / (hx * hx) + -1.0 / (hy * hy) - 2.0 / (hz * hz); in CreateSystem()
381 valA[1] = 1.0 / (hy * hy); in CreateSystem()
427 valA[0] = -2.0 / (hx * hx) + -2.0 / (hy * hy) - 1.0 / (hz * hz); in CreateSystem()
433 valA[1] = 1.0 / (hy * hy); in CreateSystem()
[all …]
H A Dex2.c152 PetscReal hx, hy; in CreateSystem() local
172 hy = 1.0 / N[1]; in CreateSystem()
237 valA[0] = -1.0 / (hx * hx) - 2.0 / (hy * hy); in CreateSystem()
242 valA[1] = 1.0 / (hy * hy); in CreateSystem()
247 valA[2] = 1.0 / (hy * hy); in CreateSystem()
258 valA[4] = 1.0 / hy; in CreateSystem()
263 valA[5] = -1.0 / hy; in CreateSystem()
271 valA[0] = -1.0 / (hx * hx) - 2.0 / (hy * hy); in CreateSystem()
276 valA[1] = 1.0 / (hy * hy); in CreateSystem()
281 valA[2] = 1.0 / (hy * hy); in CreateSystem()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/bratu2d/
H A Dbratu2df90.f9012 real(kind=8) :: hx, hy local
23 hy = 1.0/(n-1) ! y 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.py15 hy = 1.0/(ny-1) # y grid spacing
18 (2*u - uE - uW) * (hy/hx) \
19 + (2*u - uN - uS) * (hx/hy) \
20 - alpha * exp(u) * (hx*hy)
/petsc/src/snes/tutorials/
H A Dex46.c125 PetscReal D, K, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
132 hy = 1.0 / (PetscReal)(info->my - 1); in FormFunctionLocal()
133 hxdhy = hx / hy; in FormFunctionLocal()
134 hydhx = hy / hx; in FormFunctionLocal()
147 uy = (x[j + 1][i] - x[j][i]) / hy; 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
174 hy = 1.0 / (PetscReal)(info->my - 1); in FormJacobianLocal()
175 hxdhy = hx / hy; in FormJacobianLocal()
176 hydhx = hy / hx; in FormJacobianLocal()
[all …]
H A Dex55k.kokkos.cxx40 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocalVec() local
51 hy = 1.0 / (PetscReal)(info->my - 1); in FormFunctionLocalVec()
52 hxdhy = hx / hy; in FormFunctionLocalVec()
53 hydhx = hy / hx; in FormFunctionLocalVec()
67 c.y = j * hy; in FormFunctionLocalVec()
80 c.y = j * hy; in FormFunctionLocalVec()
85 c.y = j * hy; in FormFunctionLocalVec()
90 c.y = (j - 1) * hy; in FormFunctionLocalVec()
95 c.y = (j + 1) * hy; in FormFunctionLocalVec()
103 c.y = j * hy; in FormFunctionLocalVec()
[all …]
H A Dex4.c118 PetscScalar hx, hy; in FormObjectiveLocal() local
125 hy = 1.0 / (my + 1); in FormObjectiveLocal()
126 area = 0.5 * hx * hy; in FormObjectiveLocal()
155 d3 /= hy; in FormObjectiveLocal()
156 d4 /= hy; in FormObjectiveLocal()
168 d3 = (user->left[j + 1] - user->left[j + 2]) / hy; in FormObjectiveLocal()
176 d3 = (user->bottom[i + 1] - x[0][i]) / hy; in FormObjectiveLocal()
183 d4 = (user->right[j] - user->right[j + 1]) / hy; in FormObjectiveLocal()
189 d1 = (x[my - 1][i] - user->top[i + 1]) / hy; in FormObjectiveLocal()
195 d1 = (user->left[0] - user->left[1]) / hy; in FormObjectiveLocal()
[all …]
H A Dex5f.F90119 PetscReal temp1, temp, one, hx, hy
126 hy = one/((real(my) - 1))
130 temp = (real(min(j - 1, my - j)))*hy
169 PetscScalar two, one, hx, hy
184 hy = one/(real(my) - 1)
185 sc = hx*hy*lambda
186 hxdhy = hx/hy
187 hydhx = hy/hx
257 PetscScalar two, one, hx, hy, v(5)
267 hy = one/(real(my) - 1)
[all …]
H A Dex5.c76 PetscReal lambda, temp1, temp, hx, hy; in FormInitialGuess() local
84 hy = 1.0 / (PetscReal)(My - 1); in FormInitialGuess()
108 temp = (PetscReal)(PetscMin(j, My - j - 1)) * hy; in FormInitialGuess()
253 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
260 hy = 1.0 / (PetscReal)(info->my - 1); in FormFunctionLocal()
261 hxdhy = hx / hy; in FormFunctionLocal()
262 hydhx = hy / hx; in FormFunctionLocal()
270 c.y = j * hy; in FormFunctionLocal()
283 c.y = j * hy; in FormFunctionLocal()
288 c.y = j * hy; in FormFunctionLocal()
[all …]
H A Dex55.c35 PetscReal lambda, temp1, temp, hx, hy; in FormInitialGuess() local
43 hy = 1.0 / (PetscReal)(My - 1); in FormInitialGuess()
67 temp = (PetscReal)(PetscMin(j, My - j - 1)) * hy; in FormInitialGuess()
196 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
203 hy = 1.0 / (PetscReal)(info->my - 1); in FormFunctionLocal()
204 hxdhy = hx / hy; in FormFunctionLocal()
205 hydhx = hy / hx; in FormFunctionLocal()
213 c.y = j * hy; in FormFunctionLocal()
226 c.y = j * hy; in FormFunctionLocal()
231 c.y = j * hy; in FormFunctionLocal()
[all …]
H A Dex5m.m1 function localF = ex5m(localX,hx,hy,lambda)
8 sc = hx*hy*lambda; hxdhy = hx/hy; hydhx = hy/hx;
/petsc/src/ts/tests/
H A Dex21.c37 PetscReal lambda, hx, hy; in FormIFunctionLocal() local
43 hy = 1.0 / (PetscReal)(info->my - 1); in FormIFunctionLocal()
63 uyy = (un - 2.0 * u + us) / (hy * hy); in FormIFunctionLocal()
78 PetscScalar v[5], lambda, hx, hy; in FormIJacobianLocal() local
83 hy = 1.0 / (PetscReal)(info->my - 1); in FormIJacobianLocal()
99 v[k] = -1.0 / (hy * hy); in FormIJacobianLocal()
108 v[k] = shift + 2.0 / (hx * hx) + 2.0 / (hy * hy) - lambda * PetscExpScalar(x[j][i]); in FormIJacobianLocal()
117 v[k] = -1.0 / (hy * hy); in FormIJacobianLocal()
/petsc/src/binding/petsc4py/demo/legacy/wrap-swig/
H A DBratu3D.c25 PetscReal lambda,temp1,hx,hy,hz,tempk,tempj; in FormInitGuess() local
32 hy = 1.0/(PetscReal)(My-1); in FormInitGuess()
63 tempj = PetscMin((PetscReal)(PetscMin(j,My-j-1))*hy,tempk); in FormInitGuess()
85 PetscReal two = 2.0,lambda,hx,hy,hz,hxhzdhy,hyhzdhx,hxhydhz,sc; in FormFunction() local
93 hy = 1.0/(PetscReal)(My-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
33 hy = 1.0/(PetscReal)(My-1); in FormInitGuess()
64 tempj = PetscMin((PetscReal)(PetscMin(j,My-j-1))*hy,tempk); in FormInitGuess()
86 PetscReal two = 2.0,lambda,hx,hy,hz,hxhzdhy,hyhzdhx,hxhydhz,sc; in FormFunction() local
94 hy = 1.0/(PetscReal)(My-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/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
/petsc/src/tao/unconstrained/tutorials/
H A Deptorsion3.c51 PetscReal hx, hy; /* mesh spacing in x- and y-directions */ member
105 user.hy = one / (my + 1); in main()
205 PetscReal hx = user->hx, hy = user->hy, temp; in FormInitialGuess() local
212 temp = PetscMin(j + 1, ny - j) * hy; in FormInitialGuess()
260 PetscReal hx = user->hx, hy = user->hy, area, three = 3.0, p5 = 0.5; in FormFunction() local
281 dvdy = (vt - v) / hy; in FormFunction()
298 dvdy = (v - vb) / hy; in FormFunction()
308 area = p5 * hx * hy; in FormFunction()
332 PetscReal hx = user->hx, hy = user->hy; in FormGradient() local
355 dvdy = (vt - v) / hy; in FormGradient()
[all …]
H A Deptorsion1.c51 PetscReal hx, hy; /* mesh spacing in x- and y-directions */ member
103 user.hy = one / (my + 1); in main()
189 PetscReal hx = user->hx, hy = user->hy, temp; in FormInitialGuess() local
196 temp = PetscMin(j + 1, ny - j) * hy; in FormInitialGuess()
244 PetscReal hx = user->hx, hy = user->hy, area, three = 3.0, p5 = 0.5; in FormFunction() local
265 dvdy = (vt - v) / hy; in FormFunction()
282 dvdy = (v - vb) / hy; in FormFunction()
292 area = p5 * hx * hy; in FormFunction()
316 PetscReal hx = user->hx, hy = user->hy; in FormGradient() local
339 dvdy = (vt - v) / hy; in FormGradient()
[all …]
/petsc/src/ts/tutorials/
H A Dex1f.F9036 PetscReal hx, hy, hxdhy, hydhx
49 hy = 1.0/real(my - 1)
50 sc = hx*hy
51 hxdhy = hx/hy
52 hydhx = hy/hx
93 PetscReal hx, hy, hxdhy, hydhx
108 hy = 1.0/real(my - 1)
109 sc = hx*hy
110 hxdhy = hx/hy
111 hydhx = hy/hx
[all …]
H A Dex1.c183 PetscReal temp1, temp, hx, hy; in FormInitialGuess() local
192 hy = one / (PetscReal)(my - 1); in FormInitialGuess()
197 temp = (PetscReal)(PetscMin(j, my - j - 1)) * hy; in FormInitialGuess()
217 PetscReal hx, hy, hxdhy, hydhx; in FormFunction() local
227 hy = one / (PetscReal)(my - 1); in FormFunction()
228 sc = hx * hy; in FormFunction()
229 hxdhy = hx / hy; in FormFunction()
230 hydhx = hy / hx; in FormFunction()
271 PetscReal hx, hy, hxdhy, hydhx; in FormJacobian() local
279 hy = 1.0 / (PetscReal)(my - 1); in FormJacobian()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/
H A DBratu2D.F9064 PetscReal hx, hy, temp, temp1, one
68 hy = one/(dble(grd%my - 1))
72 temp = dble(min(j - 1, grd%my - j))*hy
96 PetscReal hx, hy, hxdhy, hydhx, sc, one, two
102 hy = one/dble(grd%my - 1)
103 sc = hx*hy
104 hxdhy = hx/hy
105 hydhx = hy/hx
135 PetscReal hx, hy, hxdhy, hydhx, sc, v(5), one, two
142 hy = one/dble(grd%my - 1)
[all …]
/petsc/src/binding/petsc4py/demo/legacy/bratu3d/
H A Dbratu3d.py42 hx, hy, hz = [1.0/(m-1) for m in [mx, my, mz]]
50 min_j = min(j,my-j-1)*hy
69 hx, hy, hz = [1.0/m for m in [mx, my, mz]]
70 hxhyhz = hx*hy*hz
71 hxhzdhy = hx*hz/hy;
72 hyhzdhx = hy*hz/hx;
73 hxhydhz = hx*hy/hz;
103 hx, hy, hz = [1.0/m for m in [mx, my, mz]]
104 hxhyhz = hx*hy*hz
105 hxhzdhy = hx*hz/hy;
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex46.c64 PetscReal hx = 1. / info.mx, hy = 1. / info.my; in main() local
72 v[ncols++] = 2 * (hx / hy + hy / hx); in main()
77 v[ncols++] = -hy / hx; in main()
82 v[ncols++] = -hy / hx; in main()
87 v[ncols++] = -hx / hy; in main()
92 v[ncols++] = -hx / hy; in main()

1234