365   PetscScalar uc, uxx, uyy, vc, vxx, vyy;  in IFunctionLocalPassive()  local
385       vc        = u[j][i].v;  in IFunctionLocalPassive()
386       vxx       = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx;  in IFunctionLocalPassive()
387       vyy       = (-2.0 * vc + u[j - 1][i].v + u[j + 1][i].v) * 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 …].v = udot[j][i].v - appctx->D2 * (vxx + vyy) - uc * vc * vc + (appctx->gamma + appctx->kappa) * v…  in IFunctionLocalPassive()
405   adouble       uc, uxx, uyy, vc, vxx, vyy;  in IFunctionActive()  local
481       vc          = u_a[j][i].v;  in IFunctionActive()
482       vxx         = (-2.0 * vc + u_a[j][i - 1].v + u_a[j][i + 1].v) * sx;  in IFunctionActive()
483       vyy         = (-2.0 * vc + u_a[j - 1][i].v + u_a[j + 1][i].v) * sy;  in IFunctionActive()
484 … f_a[j][i].u = udot[j][i].u - appctx->D1 * (uxx + uyy) + uc * vc * vc - appctx->gamma * (1.0 - uc);  in IFunctionActive()
485 …].v = udot[j][i].v - appctx->D2 * (vxx + vyy) - uc * vc * vc + (appctx->gamma + appctx->kappa) * v…  in IFunctionActive()
532   PetscScalar uc, uxx, uyy, vc, vxx, vyy;  in RHSFunctionPassive()  local
577       vc        = u[j][i].v;  in RHSFunctionPassive()
578       vxx       = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx;  in RHSFunctionPassive()
579       vyy       = (-2.0 * vc + u[j - 1][i].v + u[j + 1][i].v) * sy;  in RHSFunctionPassive()
580       f[j][i].u = appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc);  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
695       vc          = u_a[j][i].v;  in RHSFunctionActive()
696       vxx         = (-2.0 * vc + u_a[j][i - 1].v + u_a[j][i + 1].v) * sx;  in RHSFunctionActive()
697       vyy         = (-2.0 * vc + u_a[j - 1][i].v + u_a[j + 1][i].v) * sy;  in RHSFunctionActive()
698       f_a[j][i].u = appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc);  in RHSFunctionActive()
699       f_a[j][i].v = appctx->D2 * (vxx + vyy) + uc * vc * vc - (appctx->gamma + appctx->kappa) * vc;  in RHSFunctionActive()
792   PetscScalar uc, vc;  in IJacobianByHand()  local
848       vc = u[j][i].v;  in IJacobianByHand()
871       entries[4]   = 2.0 * appctx->D1 * (sx + sy) + vc * vc + appctx->gamma + a;  in IJacobianByHand()
874       entries[5]   = 2.0 * uc * vc;  in IJacobianByHand()
889 …   entries[4]   = 2.0 * appctx->D2 * (sx + sy) - 2.0 * uc * vc + appctx->gamma + appctx->kappa + a;  in IJacobianByHand()
891       entries[5]   = -vc * vc;  in IJacobianByHand()
922   PetscScalar uc, vc;  in RHSJacobianByHand()  local
979       vc = u[j][i].v;  in RHSJacobianByHand()
1002       entries[4]   = -2.0 * appctx->D1 * (sx + sy) - vc * vc - appctx->gamma;  in RHSJacobianByHand()
1005       entries[5]   = -2.0 * uc * vc;  in RHSJacobianByHand()
1020       entries[4]   = -2.0 * appctx->D2 * (sx + sy) + 2.0 * uc * vc - appctx->gamma - appctx->kappa;  in RHSJacobianByHand()
1022       entries[5]   = vc * vc;  in RHSJacobianByHand()