Home
last modified time | relevance | path

Searched refs:hydhx (Results 1 – 25 of 34) sorted by relevance

12

/petsc/src/snes/tutorials/
H A Dex35.c170 PetscScalar v[5], hx, hy, hxdhy, hydhx; in FormMatrix() local
178 hydhx = hy / hx; in FormMatrix()
198 v[0] = 2.0 * (hydhx + hxdhy); in FormMatrix()
207 v[1] = -hydhx; in FormMatrix()
210 v[2] = 2.0 * (hydhx + hxdhy); in FormMatrix()
213 v[3] = -hydhx; in FormMatrix()
230 PetscCall(MatZeroRowsColumnsStencil(jac, nrows, rows, 2.0 * (hydhx + hxdhy), NULL, NULL)); in FormMatrix()
248 PetscReal hx, hy, hxdhy, hydhx; in NonlinearGS() local
262 hydhx = hy / hx; in NonlinearGS()
293 uxx = (2.0 * u - x[j][i - 1] - x[j][i + 1]) * hydhx; in NonlinearGS()
[all …]
H A Dex18.c137 PetscScalar hx, hy, hxdhy, hydhx; in FormFunction() local
151 hydhx = hy / hx; in FormFunction()
282 f[j][i] = -hydhx * (fe - fw) - hxdhy * (fn - fs); in FormFunction()
296 PetscScalar one = 1.0, hx, hy, hxdhy, hydhx, t0, tn, ts, te, tw; in FormJacobian() local
311 hydhx = hy / hx; in FormJacobian()
363 v[1] = -hydhx * (dw - gw); in FormJacobian()
366 v[2] = hxdhy * (ds + dn + gs - gn) + hydhx * (dw + de + gw - ge); in FormJacobian()
369 v[3] = -hydhx * (de + ge); in FormJacobian()
402 v[0] = hxdhy * (dn - gn) + hydhx * (dw + de + gw - ge); in FormJacobian()
405 v[1] = -hydhx * (de + ge); in FormJacobian()
[all …]
H A Dex55k.kokkos.cxx40 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocalVec() local
53 hydhx = hy / hx; in FormFunctionLocalVec()
69 fv(j, i) = 2.0 * (hydhx + hxdhy) * (xv(j, i) - mms_solution); in FormFunctionLocalVec()
99 uxx = (2.0 * u - uw - ue) * hydhx; in FormFunctionLocalVec()
119 PetscReal lambda, hx, hy, hxdhy, hydhx, sc, lobj = 0; in FormObjectiveLocalVec() local
133 hydhx = hy / hx; in FormObjectiveLocalVec()
144 update += PetscRealPart((hydhx + hxdhy) * xv(j, i) * xv(j, i)); in FormObjectiveLocalVec()
159 uxux = u * (2. * u - ue - uw) * hydhx; in FormObjectiveLocalVec()
178 PetscScalar lambda, hx, hy, hxdhy, hydhx, sc; in FormJacobianLocalVec() local
195 hydhx = hy / hx; in FormJacobianLocalVec()
[all …]
H A Dex46.c125 PetscReal D, K, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
134 hydhx = hy / hx; in FormFunctionLocal()
148 uxx = (2.0 * u - x[j][i - 1] - x[j][i + 1]) * hydhx; in FormFunctionLocal()
166 PetscScalar D, K, A, v[5], hx, hy, hxdhy, hydhx, ux, uy; in FormJacobianLocal() local
176 hydhx = hy / hx; in FormJacobianLocal()
208 v[1] = -D * hydhx; in FormJacobianLocal()
211 …v[2] = D * 2.0 * (hydhx + hxdhy) + K * (funcADer(x[j][i], user) * normGradZ - A / normGradZ) *… in FormJacobianLocal()
214 v[3] = -D * hydhx + K * A * hx * hy / (2.0 * normGradZ); in FormJacobianLocal()
H A Dex58.c148 PetscScalar hx, hy, hydhx, hxdhy; in FormGradient() local
162 hydhx = hy / hx; in FormGradient()
215 df1dxc = d1 * hydhx; in FormGradient()
216 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormGradient()
218 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormGradient()
219 df5dxc = d2 * hydhx; in FormGradient()
276 PetscScalar hx, hy, hydhx, hxdhy; in FormJacobian() local
291 hydhx = hy / hx; in FormJacobian()
360 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … in FormJacobian()
361 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in FormJacobian()
[all …]
H A Dex5f.F90170 PetscScalar hxdhy, hydhx, sc
187 hydhx = hy/hx
197 uxx = hydhx*(two*u - x(i - 1, j) - x(i + 1, j))
258 PetscScalar hxdhy, hydhx, sc
270 hydhx = hy/hx
301 v(2) = -hydhx
302 v(3) = two*(hydhx + hxdhy) - sc*lambda*exp(x(i, j))
303 v(4) = -hydhx
H A Dex4.c216 PetscScalar hx, hy, hydhx, hxdhy; in FormFunctionLocal() local
224 hydhx = hy / hx; in FormFunctionLocal()
264 df1dxc = d1 * hydhx; in FormFunctionLocal()
265 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormFunctionLocal()
267 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormFunctionLocal()
268 df5dxc = d2 * hydhx; in FormFunctionLocal()
308 PetscScalar hx, hy, hydhx, hxdhy; in FormJacobianLocal() local
317 hydhx = hy / hx; in FormJacobianLocal()
371 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … in FormJacobianLocal()
372 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in FormJacobianLocal()
[all …]
H A Dex5.c253 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
262 hydhx = hy / hx; in FormFunctionLocal()
272 f[j][i] = 2.0 * (hydhx + hxdhy) * (x[j][i] - mms_solution); in FormFunctionLocal()
302 uxx = (2.0 * u - uw - ue) * hydhx; in FormFunctionLocal()
320 PetscReal lambda, hx, hy, hxdhy, hydhx, sc, lobj = 0; in FormObjectiveLocal() local
332 hydhx = hy / hx; in FormObjectiveLocal()
339 lobj += PetscRealPart((hydhx + hxdhy) * x[j][i] * x[j][i]); in FormObjectiveLocal()
354 uxux = u * (2. * u - ue - uw) * hydhx; in FormObjectiveLocal()
373 PetscScalar lambda, v[5], hx, hy, hxdhy, hydhx, sc; in FormJacobianLocal() local
388 hydhx = hy / hx; in FormJacobianLocal()
[all …]
H A Dex55.c196 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
205 hydhx = hy / hx; in FormFunctionLocal()
215 f[j][i] = 2.0 * (hydhx + hxdhy) * (x[j][i] - mms_solution); in FormFunctionLocal()
245 uxx = (2.0 * u - uw - ue) * hydhx; in FormFunctionLocal()
263 PetscReal lambda, hx, hy, hxdhy, hydhx, sc, lobj = 0; in FormObjectiveLocal() local
275 hydhx = hy / hx; in FormObjectiveLocal()
282 lobj += PetscRealPart((hydhx + hxdhy) * x[j][i] * x[j][i]); in FormObjectiveLocal()
297 uxux = u * (2. * u - ue - uw) * hydhx; in FormObjectiveLocal()
316 PetscScalar lambda, v[5], hx, hy, hxdhy, hydhx, sc; in FormJacobianLocal() local
331 hydhx = hy / hx; in FormJacobianLocal()
[all …]
H A Dex5f90t.F90229 PetscScalar two, one, hx, hy, hxdhy, hydhx, sc
239 hydhx = hy/hx
249 uxx = hydhx*(two*u - x(i - 1, j) - x(i + 1, j))
397 PetscScalar hydhx, sc, v(5)
408 hydhx = hy/hx
434 v(2) = -hydhx
435 v(3) = two*(hydhx + hxdhy) - sc*ctx%lambda*exp(x(i, j))
436 v(4) = -hydhx
H A Dex5f90.F90239 PetscScalar two, one, hx, hy, hxdhy, hydhx, sc
249 hydhx = hy/hx
259 uxx = hydhx*(two*u - x(i - 1, j) - x(i + 1, j))
408 PetscScalar hydhx, sc, v(5)
419 hydhx = hy/hx
445 v(2) = -hydhx
446 v(3) = two*(hydhx + hxdhy) - sc*ctx%lambda*exp(x(i, j))
447 v(4) = -hydhx
H A Dex19.c237 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in FormFunctionLocal() local
258 hydhx = hy * dhx; in FormFunctionLocal()
334 uxx = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx; in FormFunctionLocal()
340 uxx = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx; in FormFunctionLocal()
346 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormFunctionLocal()
352 uxx = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx; in FormFunctionLocal()
378 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in NonlinearGS() local
421 hydhx = hy * dhx; in NonlinearGS()
523 uxx = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx; in NonlinearGS()
526 dfudu = 2.0 * (hydhx + hxdhy); in NonlinearGS()
[all …]
H A Dex5m.m8 sc = hx*hy*lambda; hxdhy = hx/hy; hydhx = hy/hx;
17 uxx = (2.0*u - localX(1:m-2,2:n-1) - localX(3:m,2:n-1))*hydhx;
/petsc/src/ksp/ksp/tests/
H A Dex29.c124 PetscScalar two = 2.0, one = 1.0, v[5], hx, hy, hxdhy, hydhx, value; in FormJacobian_Grid() local
132 hydhx = hy / hx; in FormJacobian_Grid()
149 v[1] = -hydhx; in FormJacobian_Grid()
151 v[2] = two * (hydhx + hxdhy); in FormJacobian_Grid()
153 v[3] = -hydhx; in FormJacobian_Grid()
159 value = .5 * two * (hydhx + hxdhy); in FormJacobian_Grid()
162 value = .25 * two * (hydhx + hxdhy); in FormJacobian_Grid()
H A Dex26.c153 PetscScalar two = 2.0, one = 1.0, v[5], hx, hy, hxdhy, hydhx, value; in FormJacobian_Grid() local
162 hydhx = hy / hx; in FormJacobian_Grid()
179 v[1] = -hydhx; in FormJacobian_Grid()
181 v[2] = two * (hydhx + hxdhy); in FormJacobian_Grid()
183 v[3] = -hydhx; in FormJacobian_Grid()
189 value = .5 * two * (hydhx + hxdhy); in FormJacobian_Grid()
192 value = .25 * two * (hydhx + hxdhy); in FormJacobian_Grid()
H A Dex14.c344 PetscReal two = 2.0, one = 1.0, lambda, hx, hy, hxdhy, hydhx, sc; in ComputeFunction() local
355 hydhx = hy / hx; in ComputeFunction()
390 uxx = (two * u - x[row - 1] - x[row + 1]) * hydhx; in ComputeFunction()
428 PetscScalar two = 2.0, one = 1.0, lambda, v[5], hx, hy, hxdhy, hydhx, sc, *x; in ComputeJacobian() local
438 hydhx = hy / hx; in ComputeJacobian()
491 v[1] = -hydhx; in ComputeJacobian()
493 v[2] = two * (hydhx + hxdhy) - sc * lambda * PetscExpScalar(x[row]); in ComputeJacobian()
495 v[3] = -hydhx; in ComputeJacobian()
H A Dex19.c178 PetscScalar two = 2.0, one = 1.0, v[5], hx, hy, hxdhy, hydhx, value; in FormJacobian_Grid() local
186 hydhx = hy / hx; in FormJacobian_Grid()
203 v[1] = -hydhx; in FormJacobian_Grid()
205 v[2] = two * (hydhx + hxdhy); in FormJacobian_Grid()
207 v[3] = -hydhx; in FormJacobian_Grid()
213 value = .5 * two * (hydhx + hxdhy); in FormJacobian_Grid()
216 value = .25 * two * (hydhx + hxdhy); in FormJacobian_Grid()
/petsc/src/tao/complementarity/tutorials/
H A Dminsurf1.c141 PetscReal hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy; in FormConstraints() local
203 df1dxc = d1 * hydhx; in FormConstraints()
204 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormConstraints()
206 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormConstraints()
207 df5dxc = d2 * hydhx; in FormConstraints()
263 PetscReal hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy; in FormJacobian() local
335 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … in FormJacobian()
336 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in FormJacobian()
343 …hydhx * (1.0 + d7 * d7) / (f1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0… in FormJacobian()
/petsc/src/ts/tutorials/
H A Dex1f.F9036 PetscReal hx, hy, hxdhy, hydhx
52 hydhx = hy/hx
67 uxx = (-ur + two*u - ul)*hydhx
93 PetscReal hx, hy, hxdhy, hydhx
111 hydhx = hy/hx
125 v(2) = hydhx
127 v(3) = -two*(hydhx + hxdhy) + sc*lambda*exp(xx(row(1)))
129 v(4) = hydhx
H A Dex1.c217 PetscReal hx, hy, hxdhy, hydhx; in FormFunction() local
230 hydhx = hy / hx; in FormFunction()
246 uxx = (-ur + two * u - ul) * hydhx; in FormFunction()
271 PetscReal hx, hy, hxdhy, hydhx; in FormJacobian() local
282 hydhx = hy / hx; in FormJacobian()
294 v[1] = hydhx; in FormJacobian()
296 v[2] = -two * (hydhx + hxdhy) + sc * lambda * PetscExpScalar(x[row]); in FormJacobian()
298 v[3] = hydhx; in FormJacobian()
/petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/
H A DBratu2D.F9096 PetscReal hx, hy, hxdhy, hydhx, sc, one, two
105 hydhx = hy/hx
115 uxx = (two*u - x(i - 1, j) - x(i + 1, j))*hydhx
135 PetscReal hx, hy, hxdhy, hydhx, sc, v(5), one, two
145 hydhx = hy/hx
159 v(2) = -hydhx
160 v(3) = two*(hydhx + hxdhy) - lambda*exp(x(i, j))*sc
161 v(4) = -hydhx
/petsc/src/tao/unconstrained/tutorials/
H A Dminsurf1.c130 …PetscReal hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy, are… in FormFunctionGradient() local
189 df1dxc = d1 * hydhx; in FormFunctionGradient()
190 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormFunctionGradient()
192 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormFunctionGradient()
193 df5dxc = d2 * hydhx; in FormFunctionGradient()
308 PetscReal hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy; in QuadraticH() local
381 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … in QuadraticH()
382 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in QuadraticH()
389 …hydhx * (1.0 + d7 * d7) / (f1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0… in QuadraticH()
/petsc/src/ksp/ksp/tutorials/
H A Dex14f.F90152 PetscReal hy, hxdhy, hydhx, sc
164 hydhx = hy/hx
198 uxx = (two*u - xx(row - 1) - xx(row + 1))*hydhx
239 PetscScalar hydhx, sc
252 hydhx = hy/hx
299 v(2) = -hydhx
301 v(3) = two*(hydhx + hxdhy) - sc*lambda*exp(xx(row))
303 v(4) = -hydhx
/petsc/src/snes/tests/
H A Dex1f.F90254 PetscScalar hxdhy, hydhx, sc
265 hydhx = hy/hx
275 uxx = hydhx*(two*u - x(i - 1, j) - x(i + 1, j))
374 PetscScalar hxdhy, hydhx, sc, v(5)
386 hydhx = hy/hx
403 v(2) = -hydhx
404 v(3) = two*(hydhx + hxdhy) - sc*lambda*exp(x(i, j))
405 v(4) = -hydhx
H A Dex1.c371 PetscReal two = 2.0, one = 1.0, lambda, hx, hy, hxdhy, hydhx; in FormFunction() local
383 hydhx = hy / hx; in FormFunction()
406 uxx = (-ur + two * u - ul) * hydhx; in FormFunction()
439 PetscReal hx, hy, hxdhy, hydhx; in FormJacobian() local
449 hydhx = hy / hx; in FormJacobian()
468 v[1] = -hydhx; in FormJacobian()
470 v[2] = two * (hydhx + hxdhy) - sc * lambda * PetscExpScalar(x[row]); in FormJacobian()
472 v[3] = -hydhx; in FormJacobian()

12