Home
last modified time | relevance | path

Searched refs:hxdhy (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
177 hxdhy = hx / hy; in FormMatrix()
198 v[0] = 2.0 * (hydhx + hxdhy); in FormMatrix()
204 v[0] = -hxdhy; in FormMatrix()
210 v[2] = 2.0 * (hydhx + hxdhy); in FormMatrix()
216 v[4] = -hxdhy; 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
261 hxdhy = hx / hy; in NonlinearGS()
294 uyy = (2.0 * u - x[j - 1][i] - x[j + 1][i]) * hxdhy; in NonlinearGS()
[all …]
H A Dex18.c137 PetscScalar hx, hy, hxdhy, hydhx; in FormFunction() local
150 hxdhy = hx / hy; 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
310 hxdhy = hx / hy; in FormJacobian()
360 v[0] = -hxdhy * (ds - gs); in FormJacobian()
366 v[2] = hxdhy * (ds + dn + gs - gn) + hydhx * (dw + de + gw - ge); in FormJacobian()
372 v[4] = -hxdhy * (dn + gn); in FormJacobian()
402 v[0] = hxdhy * (dn - gn) + hydhx * (dw + de + gw - ge); in FormJacobian()
408 v[2] = -hxdhy * (dn + gn); in FormJacobian()
[all …]
H A Dex55k.kokkos.cxx40 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocalVec() local
52 hxdhy = hx / hy; in FormFunctionLocalVec()
69 fv(j, i) = 2.0 * (hydhx + hxdhy) * (xv(j, i) - mms_solution); in FormFunctionLocalVec()
100 uyy = (2.0 * u - un - us) * hxdhy; in FormFunctionLocalVec()
119 PetscReal lambda, hx, hy, hxdhy, hydhx, sc, lobj = 0; in FormObjectiveLocalVec() local
132 hxdhy = hx / hy; in FormObjectiveLocalVec()
144 update += PetscRealPart((hydhx + hxdhy) * xv(j, i) * xv(j, i)); in FormObjectiveLocalVec()
160 uyuy = u * (2. * u - un - us) * hxdhy; in FormObjectiveLocalVec()
178 PetscScalar lambda, hx, hy, hxdhy, hydhx, sc; in FormJacobianLocalVec() local
194 hxdhy = hx / hy; in FormJacobianLocalVec()
[all …]
H A Dex46.c125 PetscReal D, K, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
133 hxdhy = hx / hy; in FormFunctionLocal()
149 uyy = (2.0 * u - x[j - 1][i] - x[j + 1][i]) * hxdhy; in FormFunctionLocal()
166 PetscScalar D, K, A, v[5], hx, hy, hxdhy, hydhx, ux, uy; in FormJacobianLocal() local
175 hxdhy = hx / hy; in FormJacobianLocal()
205 v[0] = -D * hxdhy; in FormJacobianLocal()
211 …v[2] = D * 2.0 * (hydhx + hxdhy) + K * (funcADer(x[j][i], user) * normGradZ - A / normGradZ) *… in FormJacobianLocal()
217 v[4] = -D * hxdhy + K * A * hx * hy / (2.0 * normGradZ); in FormJacobianLocal()
H A Dex58.c148 PetscScalar hx, hy, hydhx, hxdhy; in FormGradient() local
163 hxdhy = hx / hy; in FormGradient()
216 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormGradient()
217 df3dxc = d3 * hxdhy; in FormGradient()
218 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormGradient()
220 df6dxc = d4 * hxdhy; in FormGradient()
276 PetscScalar hx, hy, hydhx, hxdhy; in FormJacobian() local
292 hxdhy = hx / hy; in FormJacobian()
362 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in FormJacobian()
363 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in FormJacobian()
[all …]
H A Dex5f.F90170 PetscScalar hxdhy, hydhx, sc
186 hxdhy = hx/hy
198 uyy = hxdhy*(two*u - x(i, j - 1) - x(i, j + 1))
258 PetscScalar hxdhy, hydhx, sc
269 hxdhy = hx/hy
300 v(1) = -hxdhy
302 v(3) = two*(hydhx + hxdhy) - sc*lambda*exp(x(i, j))
304 v(5) = -hxdhy
H A Dex4.c216 PetscScalar hx, hy, hydhx, hxdhy; in FormFunctionLocal() local
225 hxdhy = hx / hy; in FormFunctionLocal()
265 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormFunctionLocal()
266 df3dxc = d3 * hxdhy; in FormFunctionLocal()
267 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormFunctionLocal()
269 df6dxc = d4 * hxdhy; in FormFunctionLocal()
308 PetscScalar hx, hy, hydhx, hxdhy; in FormJacobianLocal() local
318 hxdhy = hx / hy; in FormJacobianLocal()
373 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in FormJacobianLocal()
374 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in FormJacobianLocal()
[all …]
H A Dex5.c253 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
261 hxdhy = hx / hy; in FormFunctionLocal()
272 f[j][i] = 2.0 * (hydhx + hxdhy) * (x[j][i] - mms_solution); in FormFunctionLocal()
303 uyy = (2.0 * u - un - us) * hxdhy; in FormFunctionLocal()
320 PetscReal lambda, hx, hy, hxdhy, hydhx, sc, lobj = 0; in FormObjectiveLocal() local
331 hxdhy = hx / hy; in FormObjectiveLocal()
339 lobj += PetscRealPart((hydhx + hxdhy) * x[j][i] * x[j][i]); in FormObjectiveLocal()
355 uyuy = u * (2. * u - un - us) * hxdhy; in FormObjectiveLocal()
373 PetscScalar lambda, v[5], hx, hy, hxdhy, hydhx, sc; in FormJacobianLocal() local
387 hxdhy = hx / hy; in FormJacobianLocal()
[all …]
H A Dex55.c196 PetscReal lambda, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local
204 hxdhy = hx / hy; in FormFunctionLocal()
215 f[j][i] = 2.0 * (hydhx + hxdhy) * (x[j][i] - mms_solution); in FormFunctionLocal()
246 uyy = (2.0 * u - un - us) * hxdhy; in FormFunctionLocal()
263 PetscReal lambda, hx, hy, hxdhy, hydhx, sc, lobj = 0; in FormObjectiveLocal() local
274 hxdhy = hx / hy; in FormObjectiveLocal()
282 lobj += PetscRealPart((hydhx + hxdhy) * x[j][i] * x[j][i]); in FormObjectiveLocal()
298 uyuy = u * (2. * u - un - us) * hxdhy; in FormObjectiveLocal()
316 PetscScalar lambda, v[5], hx, hy, hxdhy, hydhx, sc; in FormJacobianLocal() local
330 hxdhy = hx / hy; in FormJacobianLocal()
[all …]
H A Dex5f90t.F90229 PetscScalar two, one, hx, hy, hxdhy, hydhx, sc
238 hxdhy = hx/hy
250 uyy = hxdhy*(two*u - x(i, j - 1) - x(i, j + 1))
396 PetscScalar two, one, hx, hy, hxdhy
407 hxdhy = hx/hy
433 v(1) = -hxdhy
435 v(3) = two*(hydhx + hxdhy) - sc*ctx%lambda*exp(x(i, j))
437 v(5) = -hxdhy
H A Dex5f90.F90239 PetscScalar two, one, hx, hy, hxdhy, hydhx, sc
248 hxdhy = hx/hy
260 uyy = hxdhy*(two*u - x(i, j - 1) - x(i, j + 1))
407 PetscScalar two, one, hx, hy, hxdhy
418 hxdhy = hx/hy
444 v(1) = -hxdhy
446 v(3) = two*(hydhx + hxdhy) - sc*ctx%lambda*exp(x(i, j))
448 v(5) = -hxdhy
H A Dex19.c237 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in FormFunctionLocal() local
257 hxdhy = hx * dhy; in FormFunctionLocal()
335 uyy = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy; in FormFunctionLocal()
341 uyy = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy; in FormFunctionLocal()
347 uyy = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy; in FormFunctionLocal()
353 uyy = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy; in FormFunctionLocal()
378 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in NonlinearGS() local
420 hxdhy = hx * dhy; in NonlinearGS()
524 uyy = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy; 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;
18 uyy = (2.0*u - localX(2:m-1,1:n-2) - localX(2:m-1,3:n))*hxdhy;
/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
131 hxdhy = hx / hy; in FormJacobian_Grid()
147 v[0] = -hxdhy; in FormJacobian_Grid()
151 v[2] = two * (hydhx + hxdhy); in FormJacobian_Grid()
155 v[4] = -hxdhy; 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
161 hxdhy = hx / hy; in FormJacobian_Grid()
177 v[0] = -hxdhy; in FormJacobian_Grid()
181 v[2] = two * (hydhx + hxdhy); in FormJacobian_Grid()
185 v[4] = -hxdhy; 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
354 hxdhy = hx / hy; in ComputeFunction()
391 uyy = (two * u - x[row - gxm] - x[row + gxm]) * hxdhy; in ComputeFunction()
428 PetscScalar two = 2.0, one = 1.0, lambda, v[5], hx, hy, hxdhy, hydhx, sc, *x; in ComputeJacobian() local
437 hxdhy = hx / hy; in ComputeJacobian()
489 v[0] = -hxdhy; in ComputeJacobian()
493 v[2] = two * (hydhx + hxdhy) - sc * lambda * PetscExpScalar(x[row]); in ComputeJacobian()
497 v[4] = -hxdhy; in ComputeJacobian()
H A Dex19.c178 PetscScalar two = 2.0, one = 1.0, v[5], hx, hy, hxdhy, hydhx, value; in FormJacobian_Grid() local
185 hxdhy = hx / hy; in FormJacobian_Grid()
201 v[0] = -hxdhy; in FormJacobian_Grid()
205 v[2] = two * (hydhx + hxdhy); in FormJacobian_Grid()
209 v[4] = -hxdhy; 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
204 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormConstraints()
205 df3dxc = d3 * hxdhy; in FormConstraints()
206 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormConstraints()
208 df6dxc = d4 * hxdhy; in FormConstraints()
263 PetscReal hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy; in FormJacobian() local
337 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in FormJacobian()
338 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in FormJacobian()
343 …hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5) + hxdhy * (1.0… in FormJacobian()
/petsc/src/ts/tutorials/
H A Dex1f.F9036 PetscReal hx, hy, hxdhy, hydhx
51 hxdhy = hx/hy
68 uyy = (-ut + two*u - ub)*hxdhy
93 PetscReal hx, hy, hxdhy, hydhx
110 hxdhy = hx/hy
123 v(1) = hxdhy
127 v(3) = -two*(hydhx + hxdhy) + sc*lambda*exp(xx(row(1)))
131 v(5) = hxdhy
H A Dex1.c217 PetscReal hx, hy, hxdhy, hydhx; in FormFunction() local
229 hxdhy = hx / hy; in FormFunction()
247 uyy = (-ut + two * u - ub) * hxdhy; in FormFunction()
271 PetscReal hx, hy, hxdhy, hydhx; in FormJacobian() local
281 hxdhy = hx / hy; in FormJacobian()
292 v[0] = hxdhy; in FormJacobian()
296 v[2] = -two * (hydhx + hxdhy) + sc * lambda * PetscExpScalar(x[row]); in FormJacobian()
300 v[4] = hxdhy; in FormJacobian()
/petsc/src/binding/petsc4py/demo/legacy/wrap-f2py/
H A DBratu2D.F9096 PetscReal hx, hy, hxdhy, hydhx, sc, one, two
104 hxdhy = hx/hy
116 uyy = (two*u - x(i, j - 1) - x(i, j + 1))*hxdhy
135 PetscReal hx, hy, hxdhy, hydhx, sc, v(5), one, two
144 hxdhy = hx/hy
158 v(1) = -hxdhy
160 v(3) = two*(hydhx + hxdhy) - lambda*exp(x(i, j))*sc
162 v(5) = -hxdhy
/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
190 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormFunctionGradient()
191 df3dxc = d3 * hxdhy; in FormFunctionGradient()
192 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormFunctionGradient()
194 df6dxc = d4 * hxdhy; in FormFunctionGradient()
308 PetscReal hx = 1.0 / (mx + 1), hy = 1.0 / (my + 1), hydhx = hy / hx, hxdhy = hx / hy; in QuadraticH() local
383 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in QuadraticH()
384 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in QuadraticH()
389 …hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5) + hxdhy * (1.0… in QuadraticH()
/petsc/src/ksp/ksp/tutorials/
H A Dex14f.F90152 PetscReal hy, hxdhy, hydhx, sc
163 hxdhy = hx/hy
199 uyy = (two*u - xx(row - gxm) - xx(row + gxm))*hxdhy
238 PetscScalar v(5), hx, hy, hxdhy
251 hxdhy = hx/hy
297 v(1) = -hxdhy
301 v(3) = two*(hydhx + hxdhy) - sc*lambda*exp(xx(row))
305 v(5) = -hxdhy
/petsc/src/snes/tests/
H A Dex1f.F90254 PetscScalar hxdhy, hydhx, sc
264 hxdhy = hx/hy
276 uyy = hxdhy*(two*u - x(i, j - 1) - x(i, j + 1))
374 PetscScalar hxdhy, hydhx, sc, v(5)
385 hxdhy = hx/hy
402 v(1) = -hxdhy
404 v(3) = two*(hydhx + hxdhy) - sc*lambda*exp(x(i, j))
406 v(5) = -hxdhy
H A Dex1.c371 PetscReal two = 2.0, one = 1.0, lambda, hx, hy, hxdhy, hydhx; in FormFunction() local
382 hxdhy = hx / hy; in FormFunction()
407 uyy = (-ut + two * u - ub) * hxdhy; in FormFunction()
439 PetscReal hx, hy, hxdhy, hydhx; in FormJacobian() local
448 hxdhy = hx / hy; in FormJacobian()
466 v[0] = -hxdhy; in FormJacobian()
470 v[2] = two * (hydhx + hxdhy) - sc * lambda * PetscExpScalar(x[row]); in FormJacobian()
474 v[4] = -hxdhy; in FormJacobian()

12