| /petsc/src/snes/tests/ |
| H A D | ex12f.F90 | 32 PetscScalar, pointer :: vf2(:), vxx(:), vff(:)
47 PetscCall(VecGetArrayRead(ctx%xl, vxx, ierr))
71 vff(1) = vxx(1)
77 … vff(i - s + 1) = d*(vxx(i) - 2.0*vxx(i + 1) + vxx(i + 2)) + vxx(i + 1)*vxx(i + 1) - vF2(i - s + 1)
81 vff(n - s) = vxx(n) - 1.0
85 PetscCall(VecRestoreArrayRead(ctx%xl, vxx, ierr))
148 PetscScalar, pointer :: vxx(:)
156 PetscCall(VecGetArrayRead(x, vxx, ierr))
181 A = -2.0*d + 2.0*vxx(i + 1)
184 PetscCall(VecRestoreArrayRead(x, vxx, ierr))
|
| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | ex5adj_mf.c | 65 PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb; in MyRHSMatMultTranspose() local
104 vxx = (-2.0 * vcb + x[j][i - 1].v + x[j][i + 1].v) * sx; in MyRHSMatMultTranspose()
107 …y[j][i].v = appctx->D2 * (vxx + vyy) - 2.0 * uc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appc… in MyRHSMatMultTranspose()
124 PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb; in MyIMatMultTranspose() local
163 vxx = (-2.0 * vcb + x[j][i - 1].v + x[j][i + 1].v) * sx; in MyIMatMultTranspose()
166 …y[j][i].v = appctx->D2 * (vxx + vyy) - 2.0 * uc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appc… in MyIMatMultTranspose()
185 PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb; in MyIMatMult() local
224 vxx = (-2.0 * vcb + x[j][i - 1].v + x[j][i + 1].v) * sx; in MyIMatMult()
227 …y[j][i].v = appctx->D2 * (vxx + vyy) + vc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appctx->ka… in MyIMatMult()
|
| H A D | reaction_diffusion.c | 33 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in RHSFunction() local
75 vxx = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx; in RHSFunction()
78 f[j][i].v = appctx->D2 * (vxx + vyy) + uc * vc * vc - (appctx->gamma + appctx->kappa) * vc; in RHSFunction()
241 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in IFunction() local
284 vxx = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx; in IFunction()
287 …f[j][i].v = udot[j][i].v - (appctx->D2 * (vxx + vyy) + uc * vc * vc - (appctx->gamma + appctx->kap… in IFunction()
|
| /petsc/src/ts/tutorials/autodiff/ |
| H A D | adr_ex5adj_mf.cxx | 253 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in IFunctionLocalPassive() local
274 vxx = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx; in IFunctionLocalPassive()
277 …f[j][i].v = udot[j][i].v - appctx->D2 * (vxx + vyy) - uc * vc * vc + (appctx->gamma + appctx->kapp… in IFunctionLocalPassive()
293 adouble uc, uxx, uyy, vc, vxx, vyy; in IFunctionActive() local
370 vxx = (-2.0 * vc + u_a[j][i - 1].v + u_a[j][i + 1].v) * sx; in IFunctionActive()
373 …f_a[j][i].v = udot[j][i].v - appctx->D2 * (vxx + vyy) - uc * vc * vc + (appctx->gamma + appctx->ka… in IFunctionActive()
|
| H A D | adr_ex5adj.cxx | 365 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in IFunctionLocalPassive() local
386 vxx = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx; in IFunctionLocalPassive()
389 …f[j][i].v = udot[j][i].v - appctx->D2 * (vxx + vyy) - uc * vc * vc + (appctx->gamma + appctx->kapp… in IFunctionLocalPassive()
405 adouble uc, uxx, uyy, vc, vxx, vyy; in IFunctionActive() local
482 vxx = (-2.0 * vc + u_a[j][i - 1].v + u_a[j][i + 1].v) * sx; in IFunctionActive()
485 …f_a[j][i].v = udot[j][i].v - appctx->D2 * (vxx + vyy) - uc * vc * vc + (appctx->gamma + appctx->ka… in IFunctionActive()
532 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in RHSFunctionPassive() local
578 vxx = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx; in RHSFunctionPassive()
581 f[j][i].v = appctx->D2 * (vxx + vyy) + uc * vc * vc - (appctx->gamma + appctx->kappa) * vc; in RHSFunctionPassive()
613 adouble uc, uxx, uyy, vc, vxx, vyy; in RHSFunctionActive() local
[all …]
|