| /petsc/src/snes/tutorials/ |
| H A D | ex15.c | 411 PetscReal hx, hy, dhx, dhy, sc; in FormFunctionLocal() local
419 dhx = 1 / hx; in FormFunctionLocal()
430 …dhx * (x[j][i + 1] - x[j][i]), uy_E = 0.25 * dhy * (x[j + 1][i] + x[j + 1][i + 1] - x[j - 1][i] - … in FormFunctionLocal()
484 PetscReal hx, hy, hxdhy, hydhx, dhx, dhy, sc; in FormJacobianLocal() local
491 dhx = 1 / hx; in FormJacobianLocal()
516 …dhx * (x[j][i + 1] - x[j][i]), uy_E = 0.25 * dhy * (x[j + 1][i] + x[j + 1][i + 1] - x[j - 1][i] - … in FormJacobianLocal()
729 PetscReal hx, hy, hxdhy, hydhx, dhx, dhy, sc; in NonlinearGS() local
744 dhx = 1 / hx; in NonlinearGS()
781 …] - x[j - 1][i]), ux_N = 0.25 * dhx * (x[j][i + 1] + x[j + 1][i + 1] - x[j][i - 1] - x[j + 1][i - … in NonlinearGS()
784 …const PetscScalar ux_E = dhx * (u_E - u), ux_W = dhx * (u - u_W), uy_N = dhy * (u_N - u), uy_S = d… in NonlinearGS()
|
| H A D | ex19.c | 237 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in FormFunctionLocal() local
253 dhx = (PetscReal)(info->mx - 1); in FormFunctionLocal()
255 hx = 1.0 / dhx; in FormFunctionLocal()
258 hydhx = hy * dhx; in FormFunctionLocal()
299 f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx; in FormFunctionLocal()
312 f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx; in FormFunctionLocal()
378 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in NonlinearGS() local
416 dhx = (PetscReal)(info.mx - 1); in NonlinearGS()
418 hx = 1.0 / dhx; in NonlinearGS()
421 hydhx = hy * dhx; in NonlinearGS()
[all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex5.c | 610 PetscReal dhx, dhy; in RhsFunc() local
641 …dhx = (PetscReal)(Mx - 1) / (5000 * (Mx - 1)); /* dhx = 1/dx; assume 2D space domain: [0.0, 1.e5] … in RhsFunc()
692 …Frhs[j][i].u = -u_plus * (u - X[j][i - 1].u) * dhx - u_minus * (X[j][i + 1].u - u) * dhx … in RhsFunc()
694 … (Ts * (X[j][i + 1].p - X[j][i - 1].p) * 0.5 * dhx + p * 0 * (X[j][i + 1].Ts - X[j][i - 1].Ts) * 0… in RhsFunc()
699 …Frhs[j][i].v = -u_plus * (v - X[j][i - 1].v) * dhx - u_minus * (X[j][i + 1].v - v) * dhx … in RhsFunc()
707 …- u_plus * (Ts - X[j][i - 1].Ts) * dhx - u_minus * (X[j][i + 1].Ts - Ts) * dhx /* - u*(dTs/dx) ad… in RhsFunc()
709 …+ diffconst * ((X[j][i + 1].Ts - 2 * Ts + X[j][i - 1].Ts) * dhx * dhx /* + D(Ts_xx + Ts_y… in RhsFunc()
713 …Frhs[j][i].p = -u_plus * (p - X[j][i - 1].p) * dhx - u_minus * (X[j][i + 1].p - p) * dhx /* - u*(d… in RhsFunc()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex29.c | 225 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx, scale; in FormIFunctionLocal() local
240 dhx = 1. / hx; in FormIFunctionLocal()
255 if (i == 0) f[j][i].u += (x[j][i].u - x[j][i + 1].u) * dhx; in FormIFunctionLocal()
256 else if (i == info->mx - 1) f[j][i].u += (x[j][i].u - x[j][i - 1].u) * dhx; in FormIFunctionLocal()
263 bx = 0.5 * (c[j][i + 1].beta - c[j][i - 1].beta) * dhx; in FormIFunctionLocal()
266 ux = 0.5 * (x[j][i + 1].u - x[j][i - 1].u) * dhx; in FormIFunctionLocal()
|
| H A D | ex26.c | 261 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in FormIFunctionLocal() local
277 dhx = (PetscReal)(info->mx - 1); in FormIFunctionLocal()
279 hx = 1.0 / dhx; in FormIFunctionLocal()
282 hydhx = hy * dhx; in FormIFunctionLocal()
323 f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx; in FormIFunctionLocal()
336 f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx; in FormIFunctionLocal()
|
| H A D | ex27.c | 190 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx, scale; in FormIFunctionLocal() local
199 dhx = 1. / hx; in FormIFunctionLocal()
219 sxp = (x[j][i + 1].sp[l] - x[j][i].sp[l]) * dhx; in FormIFunctionLocal()
222 sxp = (x[j][i].sp[l] - x[j][i - 1].sp[l]) * dhx; in FormIFunctionLocal()
225 sxm = (x[j][i + 1].sp[l] - x[j][i].sp[l]) * dhx; in FormIFunctionLocal()
226 sxp = (x[j][i].sp[l] - x[j][i - 1].sp[l]) * dhx; in FormIFunctionLocal()
|
| /petsc/src/snes/tests/ |
| H A D | ex69.c | 212 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in FormFunctionLocal() local
234 dhx = (PetscReal)(info->mx - 1); in FormFunctionLocal()
236 hx = 1.0 / dhx; in FormFunctionLocal()
239 hydhx = hy * dhx; in FormFunctionLocal()
280 f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx; in FormFunctionLocal()
293 f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx; in FormFunctionLocal()
|