| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | reaction_diffusion.c | 33 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in RHSFunction() local
71 uc = u[j][i].u; in RHSFunction()
72 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in RHSFunction()
73 uyy = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy; in RHSFunction()
77 f[j][i].u = appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc); in RHSFunction()
78 f[j][i].v = appctx->D2 * (vxx + vyy) + uc * vc * vc - (appctx->gamma + appctx->kappa) * vc; in RHSFunction()
98 PetscScalar uc, vc; in RHSJacobian() local
153 uc = u[j][i].u; in RHSJacobian()
180 entries[5] = -2.0 * uc * vc; in RHSJacobian()
195 entries[4] = -2.0 * appctx->D2 * (sx + sy) + 2.0 * uc * vc - appctx->gamma - appctx->kappa; in RHSJacobian()
[all …]
|
| H A D | ex5adj_mf.c | 65 PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb; in MyRHSMatMultTranspose() local
98 uc = u[j][i].u; 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
157 uc = u[j][i].u; 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
218 uc = u[j][i].u; in MyIMatMult()
226 y[j][i].u = appctx->D1 * (uxx + uyy) - (vc * vc + appctx->gamma) * ucb - 2.0 * uc * vc * vcb; in MyIMatMult()
227 …y[j][i].v = appctx->D2 * (vxx + vyy) + vc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appctx->ka… in MyIMatMult()
|
| /petsc/src/ts/tutorials/autodiff/ |
| H A D | adr_ex5adj.cxx | 365 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in IFunctionLocalPassive() local
382 uc = u[j][i].u; in IFunctionLocalPassive()
383 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in IFunctionLocalPassive()
384 uyy = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy; in IFunctionLocalPassive()
388 … f[j][i].u = udot[j][i].u - appctx->D1 * (uxx + uyy) + uc * vc * vc - appctx->gamma * (1.0 - uc); 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
478 uc = u_a[j][i].u; in IFunctionActive()
479 uxx = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx; in IFunctionActive()
480 uyy = (-2.0 * uc + u_a[j - 1][i].u + u_a[j + 1][i].u) * sy; in IFunctionActive()
[all …]
|
| H A D | adr_ex5adj_mf.cxx | 253 PetscScalar uc, uxx, uyy, vc, vxx, vyy; in IFunctionLocalPassive() local
270 uc = u[j][i].u; in IFunctionLocalPassive()
271 uxx = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx; in IFunctionLocalPassive()
272 uyy = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy; in IFunctionLocalPassive()
276 … f[j][i].u = udot[j][i].u - appctx->D1 * (uxx + uyy) + uc * vc * vc - appctx->gamma * (1.0 - uc); 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
366 uc = u_a[j][i].u; in IFunctionActive()
367 uxx = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx; in IFunctionActive()
368 uyy = (-2.0 * uc + u_a[j - 1][i].u + u_a[j + 1][i].u) * sy; in IFunctionActive()
[all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex32.c | 33 PetscReal uc; member
76 for (i = 0; i < lsize; i++) fvalue[i] = PetscRealPart(u[i]) - app->uc; in EventFunction()
205 app.uc = 0.5; in main()
206 PetscCall(PetscOptionsReal("-uc", "", "", app.uc, &app.uc, NULL)); in main()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex36.c | 438 Vec uc; in CoarseTest() local
450 PetscCall(DMCreateGlobalVector(dmc, &uc)); in CoarseTest()
451 PetscCall(PetscObjectSetName((PetscObject)uc, "potential")); in CoarseTest()
452 PetscCall(PetscObjectSetOptionsPrefix((PetscObject)uc, "coarse_")); in CoarseTest()
456 PetscCall(MatMult(P, u, uc)); in CoarseTest()
463 PetscCall(MatMultTranspose(In, u, uc)); in CoarseTest()
464 PetscCall(VecPointwiseMult(uc, sc, uc)); in CoarseTest()
469 PetscCall(CompareView(uc)); in CoarseTest()
472 PetscCall(VecDestroy(&uc)); in CoarseTest()
|
| /petsc/src/ts/tests/ |
| H A D | ex24.c | 22 PetscScalar uc, vc; in RHSJacobian() local
77 uc = u[j][i].u; in RHSJacobian()
104 entries[5] = -2.0 * uc * vc; in RHSJacobian()
118 entries[4] = -2.0 * appctx->D2 * (sx + sy) + 2.0 * uc * vc - appctx->gamma - appctx->kappa; in RHSJacobian()
|
| /petsc/src/ksp/pc/impls/tfs/ |
| H A D | xxt.c | 407 static PetscErrorCode do_xxt_solve(xxt_ADT xxt_handle, PetscScalar *uc) in do_xxt_solve() argument
431 PetscCallBLAS("BLASdot", *uu_ptr++ = BLASdot_(&dlen, uc + off, &i1, x_ptr, &i1)); in do_xxt_solve()
438 PetscCall(PCTFS_rvec_zero(uc, n)); in do_xxt_solve()
445 PetscCallBLAS("BLASaxpy", BLASaxpy_(&dlen, uu_ptr++, x_ptr, &i1, uc + off, &i1)); in do_xxt_solve()
|
| H A D | xyt.c | 516 static PetscErrorCode do_xyt_solve(xyt_ADT xyt_handle, PetscScalar *uc) in do_xyt_solve() argument
542 PetscCallBLAS("BLASdot", *uu_ptr++ = BLASdot_(&dlen, uc + off, &i1, y_ptr, &i1)); in do_xyt_solve()
548 PetscCall(PCTFS_rvec_zero(uc, n)); in do_xyt_solve()
555 PetscCallBLAS("BLASaxpy", BLASaxpy_(&dlen, uu_ptr++, x_ptr, &i1, uc + off, &i1)); in do_xyt_solve()
|