| /petsc/src/snes/tutorials/ |
| H A D | ex15.c | 411 PetscReal hx, hy, dhx, dhy, sc; in FormFunctionLocal() local
420 dhy = 1 / hy; in FormFunctionLocal()
430 …dhy * (x[j + 1][i] + x[j + 1][i + 1] - x[j - 1][i] - x[j - 1][i + 1]), ux_W = dhx * (x[j][i] - x[j… in FormFunctionLocal()
484 PetscReal hx, hy, hxdhy, hydhx, dhx, dhy, sc; in FormJacobianLocal() local
492 dhy = 1 / hy; in FormJacobianLocal()
516 …dhy * (x[j + 1][i] + x[j + 1][i + 1] - x[j - 1][i] - x[j - 1][i + 1]), ux_W = dhx * (x[j][i] - x[j… in FormJacobianLocal()
729 PetscReal hx, hy, hxdhy, hydhx, dhx, dhy, sc; in NonlinearGS() local
745 dhy = 1 / hy; in NonlinearGS()
781 …const PetscScalar uy_E = 0.25 * dhy * (x[j + 1][i] + x[j + 1][i + 1] - x[j - 1][i] - x[j - 1][i + … in NonlinearGS()
784 …cScalar ux_E = dhx * (u_E - u), ux_W = dhx * (u - u_W), uy_N = dhy * (u_N - u), uy_S = dhy * (u - … in NonlinearGS()
|
| H A D | ex19.c | 237 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in FormFunctionLocal() local
254 dhy = (PetscReal)(info->my - 1); in FormFunctionLocal()
256 hy = 1.0 / dhy; in FormFunctionLocal()
257 hxdhy = hx * dhy; in FormFunctionLocal()
273 f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy; in FormFunctionLocal()
286 f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy; in FormFunctionLocal()
378 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in NonlinearGS() local
417 dhy = (PetscReal)(info.my - 1); in NonlinearGS()
419 hy = 1.0 / dhy; in NonlinearGS()
420 hxdhy = hx * dhy; in NonlinearGS()
[all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex5.c | 610 PetscReal dhx, dhy; in RhsFunc() local
642 dhy = (PetscReal)(My - 1) / (5000 * (Mx - 1)); /* dhy = 1/dy; */ in RhsFunc()
693 …- v_plus * (u - X[j - 1][i].u) * dhy - v_minus * (X[j + 1][i].u - u) * dhy … in RhsFunc()
700 …- v_plus * (v - X[j - 1][i].v) * dhy - v_minus * (X[j + 1][i].v - v) * dhy … in RhsFunc()
701 … (Ts * (X[j + 1][i].p - X[j - 1][i].p) * 0.5 * dhy + p * 0 * (X[j + 1][i].Ts - X[j - 1][i].Ts) * 0… in RhsFunc()
708 …- v_plus * (Ts - X[j - 1][i].Ts) * dhy - v_minus * (X[j + 1][i].Ts - Ts) * dhy /* - v*(dTs/dy) ad… in RhsFunc()
710 + (X[j + 1][i].Ts - 2 * Ts + X[j - 1][i].Ts) * dhy * dhy); in RhsFunc()
714 … - v_plus * (p - X[j - 1][i].p) * dhy - v_minus * (X[j + 1][i].p - p) * dhy; /* - v*(dp/dy) */ in RhsFunc()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex29.c | 225 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx, scale; in FormIFunctionLocal() local
241 dhy = 1. / hy; in FormIFunctionLocal()
257 else if (j == 0) f[j][i].u += (x[j][i].u - x[j + 1][i].u) * dhy; in FormIFunctionLocal()
258 else if (j == info->my - 1) f[j][i].u += (x[j][i].u - x[j - 1][i].u) * dhy; in FormIFunctionLocal()
264 by = 0.5 * (c[j + 1][i].beta - c[j - 1][i].beta) * dhy; in FormIFunctionLocal()
267 uy = 0.5 * (x[j + 1][i].u - x[j - 1][i].u) * dhy; in FormIFunctionLocal()
|
| H A D | ex27.c | 190 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx, scale; in FormIFunctionLocal() local
200 dhy = 1. / hy; in FormIFunctionLocal()
208 vy = 0.0 * dhy; in FormIFunctionLocal()
229 syp = (x[j + 1][i].sp[l] - x[j][i].sp[l]) * dhy; in FormIFunctionLocal()
232 syp = (x[j][i].sp[l] - x[j - 1][i].sp[l]) * dhy; in FormIFunctionLocal()
235 sym = (x[j + 1][i].sp[l] - x[j][i].sp[l]) * dhy; in FormIFunctionLocal()
236 syp = (x[j][i].sp[l] - x[j - 1][i].sp[l]) * dhy; in FormIFunctionLocal()
|
| H A D | ex26.c | 261 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in FormIFunctionLocal() local
278 dhy = (PetscReal)(info->my - 1); in FormIFunctionLocal()
280 hy = 1.0 / dhy; in FormIFunctionLocal()
281 hxdhy = hx * dhy; in FormIFunctionLocal()
297 f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy; in FormIFunctionLocal()
310 f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy; in FormIFunctionLocal()
|
| /petsc/src/snes/tests/ |
| H A D | ex69.c | 212 PetscReal hx, hy, dhx, dhy, hxdhy, hydhx; in FormFunctionLocal() local
235 dhy = (PetscReal)(info->my - 1); in FormFunctionLocal()
237 hy = 1.0 / dhy; in FormFunctionLocal()
238 hxdhy = hx * dhy; in FormFunctionLocal()
254 f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy; in FormFunctionLocal()
267 f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy; in FormFunctionLocal()
|