364   PetscReal   hx, hy, sx, sy;  in IFunctionLocalPassive()  local
369   sx = 1.0 / (hx * hx);  in IFunctionLocalPassive()
383       uxx       = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx;  in IFunctionLocalPassive()
386       vxx       = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx;  in IFunctionLocalPassive()
404   PetscReal     hx, hy, sx, sy;  in IFunctionActive()  local
414   sx  = 1.0 / (hx * hx);  in IFunctionActive()
479       uxx         = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx;  in IFunctionActive()
482       vxx         = (-2.0 * vc + u_a[j][i - 1].v + u_a[j][i + 1].v) * sx;  in IFunctionActive()
531   PetscReal   hx, hy, sx, sy;  in RHSFunctionPassive()  local
540   sx = 1.0 / (hx * hx);  in RHSFunctionPassive()
575       uxx       = (-2.0 * uc + u[j][i - 1].u + u[j][i + 1].u) * sx;  in RHSFunctionPassive()
578       vxx       = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx;  in RHSFunctionPassive()
611   PetscReal hx, hy, sx, sy;  in RHSFunctionActive()  local
621   sx = 1.0 / (hx * hx);  in RHSFunctionActive()
693       uxx         = (-2.0 * uc + u_a[j][i - 1].u + u_a[j][i + 1].u) * sx;  in RHSFunctionActive()
696       vxx         = (-2.0 * vc + u_a[j][i - 1].v + u_a[j][i + 1].v) * sx;  in RHSFunctionActive()
791   PetscReal   hx, hy, sx, sy;  in IJacobianByHand()  local
804   sx = 1.0 / (hx * hx);  in IJacobianByHand()
865       entries[2]   = -appctx->D1 * sx;  in IJacobianByHand()
868       entries[3]   = -appctx->D1 * sx;  in IJacobianByHand()
871       entries[4]   = 2.0 * appctx->D1 * (sx + sy) + vc * vc + appctx->gamma + a;  in IJacobianByHand()
885       entries[2]   = -appctx->D2 * sx;  in IJacobianByHand()
887       entries[3]   = -appctx->D2 * sx;  in IJacobianByHand()
889 …   entries[4]   = 2.0 * appctx->D2 * (sx + sy) - 2.0 * uc * vc + appctx->gamma + appctx->kappa + a;  in IJacobianByHand()
921   PetscReal   hx, hy, sx, sy;  in RHSJacobianByHand()  local
934   sx = 1.0 / (hx * hx);  in RHSJacobianByHand()
996       entries[2]   = appctx->D1 * sx;  in RHSJacobianByHand()
999       entries[3]   = appctx->D1 * sx;  in RHSJacobianByHand()
1002       entries[4]   = -2.0 * appctx->D1 * (sx + sy) - vc * vc - appctx->gamma;  in RHSJacobianByHand()
1016       entries[2]   = appctx->D2 * sx;  in RHSJacobianByHand()
1018       entries[3]   = appctx->D2 * sx;  in RHSJacobianByHand()
1020       entries[4]   = -2.0 * appctx->D2 * (sx + sy) + 2.0 * uc * vc - appctx->gamma - appctx->kappa;  in RHSJacobianByHand()