364   PetscReal   hx, hy, sx, sy;  in IFunctionLocalPassive()  local
371   sy = 1.0 / (hy * hy);  in IFunctionLocalPassive()
384       uyy       = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy;  in IFunctionLocalPassive()
387       vyy       = (-2.0 * vc + u[j - 1][i].v + u[j + 1][i].v) * sy;  in IFunctionLocalPassive()
404   PetscReal     hx, hy, sx, sy;  in IFunctionActive()  local
416   sy  = 1.0 / (hy * hy);  in IFunctionActive()
480       uyy         = (-2.0 * uc + u_a[j - 1][i].u + u_a[j + 1][i].u) * sy;  in IFunctionActive()
483       vyy         = (-2.0 * vc + u_a[j - 1][i].v + u_a[j + 1][i].v) * sy;  in IFunctionActive()
531   PetscReal   hx, hy, sx, sy;  in RHSFunctionPassive()  local
542   sy = 1.0 / (hy * hy);  in RHSFunctionPassive()
576       uyy       = (-2.0 * uc + u[j - 1][i].u + u[j + 1][i].u) * sy;  in RHSFunctionPassive()
579       vyy       = (-2.0 * vc + u[j - 1][i].v + u[j + 1][i].v) * sy;  in RHSFunctionPassive()
611   PetscReal hx, hy, sx, sy;  in RHSFunctionActive()  local
623   sy = 1.0 / (hy * hy);  in RHSFunctionActive()
694       uyy         = (-2.0 * uc + u_a[j - 1][i].u + u_a[j + 1][i].u) * sy;  in RHSFunctionActive()
697       vyy         = (-2.0 * vc + u_a[j - 1][i].v + u_a[j + 1][i].v) * sy;  in RHSFunctionActive()
791   PetscReal   hx, hy, sx, sy;  in IJacobianByHand()  local
806   sy = 1.0 / (hy * hy);  in IJacobianByHand()
859       entries[0]   = -appctx->D1 * sy;  in IJacobianByHand()
862       entries[1]   = -appctx->D1 * sy;  in IJacobianByHand()
871       entries[4]   = 2.0 * appctx->D1 * (sx + sy) + vc * vc + appctx->gamma + a;  in IJacobianByHand()
881       entries[0]   = -appctx->D2 * sy;  in IJacobianByHand()
883       entries[1]   = -appctx->D2 * sy;  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
936   sy = 1.0 / (hy * hy);  in RHSJacobianByHand()
990       entries[0]   = appctx->D1 * sy;  in RHSJacobianByHand()
993       entries[1]   = appctx->D1 * sy;  in RHSJacobianByHand()
1002       entries[4]   = -2.0 * appctx->D1 * (sx + sy) - vc * vc - appctx->gamma;  in RHSJacobianByHand()
1012       entries[0]   = appctx->D2 * sy;  in RHSJacobianByHand()
1014       entries[1]   = appctx->D2 * sy;  in RHSJacobianByHand()
1020       entries[4]   = -2.0 * appctx->D2 * (sx + sy) + 2.0 * uc * vc - appctx->gamma - appctx->kappa;  in RHSJacobianByHand()