65   PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb;  in MyRHSMatMultTranspose()  local
99       ucb       = x[j][i].u;  in MyRHSMatMultTranspose()
100       uxx       = (-2.0 * ucb + x[j][i - 1].u + x[j][i + 1].u) * sx;  in MyRHSMatMultTranspose()
101       uyy       = (-2.0 * ucb + x[j - 1][i].u + x[j + 1][i].u) * sy;  in MyRHSMatMultTranspose()
106       y[j][i].u = appctx->D1 * (uxx + uyy) - ucb * (vc * vc + appctx->gamma) + vc * vc * vcb;  in MyRHSMatMultTranspose()
107 …y[j][i].v = appctx->D2 * (vxx + vyy) - 2.0 * uc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appc…  in MyRHSMatMultTranspose()
124   PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb;  in MyIMatMultTranspose()  local
158       ucb       = x[j][i].u;  in MyIMatMultTranspose()
159       uxx       = (-2.0 * ucb + x[j][i - 1].u + x[j][i + 1].u) * sx;  in MyIMatMultTranspose()
160       uyy       = (-2.0 * ucb + x[j - 1][i].u + x[j + 1][i].u) * sy;  in MyIMatMultTranspose()
165       y[j][i].u = appctx->D1 * (uxx + uyy) - ucb * (vc * vc + appctx->gamma) + vc * vc * vcb;  in MyIMatMultTranspose()
166 …y[j][i].v = appctx->D2 * (vxx + vyy) - 2.0 * uc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appc…  in MyIMatMultTranspose()
167       y[j][i].u = mctx->shift * ucb - y[j][i].u;  in MyIMatMultTranspose()
185   PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb;  in MyIMatMult()  local
219       ucb       = x[j][i].u;  in MyIMatMult()
220       uxx       = (-2.0 * ucb + x[j][i - 1].u + x[j][i + 1].u) * sx;  in MyIMatMult()
221       uyy       = (-2.0 * ucb + x[j - 1][i].u + x[j + 1][i].u) * sy;  in MyIMatMult()
226       y[j][i].u = appctx->D1 * (uxx + uyy) - (vc * vc + appctx->gamma) * ucb - 2.0 * uc * vc * vcb;  in MyIMatMult()
227 …y[j][i].v = appctx->D2 * (vxx + vyy) + vc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appctx->ka…  in MyIMatMult()
228       y[j][i].u = mctx->shift * ucb - y[j][i].u;  in MyIMatMult()