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()
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 …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
574       uc        = u[j][i].u;  in RHSFunctionPassive()
575       uxx       = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx;  in RHSFunctionPassive()
576       uyy       = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * 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
692       uc          = u_a[j][i].u;  in RHSFunctionActive()
693       uxx         = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx;  in RHSFunctionActive()
694       uyy         = (-2.0 * uc + u_a[j - 1][i].u + u_a[j + 1][i].u) * 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
847       uc = u[j][i].u;  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()
922   PetscScalar uc, vc;  in RHSJacobianByHand()  local
978       uc = u[j][i].u;  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()