33   PetscScalar uc, uxx, uyy, vc, vxx, vyy;  in RHSFunction()  local
74       vc        = u[j][i].v;  in RHSFunction()
75       vxx       = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx;  in RHSFunction()
76       vyy       = (-2.0 * vc + u[j - 1][i].v + u[j + 1][i].v) * 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
154       vc = u[j][i].v;  in RHSJacobian()
177       entries[4]   = -2.0 * appctx->D1 * (sx + sy) - vc * vc - appctx->gamma;  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()
197       entries[5]   = vc * vc;  in RHSJacobian()
241   PetscScalar uc, uxx, uyy, vc, vxx, vyy;  in IFunction()  local
283       vc        = u[j][i].v;  in IFunction()
284       vxx       = (-2.0 * vc + u[j][i - 1].v + u[j][i + 1].v) * sx;  in IFunction()
285       vyy       = (-2.0 * vc + u[j - 1][i].v + u[j + 1][i].v) * sy;  in IFunction()
286 … f[j][i].u = udot[j][i].u - (appctx->D1 * (uxx + uyy) - uc * vc * vc + appctx->gamma * (1.0 - uc));  in IFunction()
287 ….v = udot[j][i].v - (appctx->D2 * (vxx + vyy) + uc * vc * vc - (appctx->gamma + appctx->kappa) * v…  in IFunction()
308   PetscScalar uc, vc;  in IJacobian()  local
364       vc = u[j][i].v;  in IJacobian()
387       entries[4]   = 2.0 * appctx->D1 * (sx + sy) + vc * vc + appctx->gamma + a;  in IJacobian()
390       entries[5]   = 2.0 * uc * vc;  in IJacobian()
405 …   entries[4]   = 2.0 * appctx->D2 * (sx + sy) - 2.0 * uc * vc + appctx->gamma + appctx->kappa + a;  in IJacobian()
407       entries[5]   = -vc * vc;  in IJacobian()