65   PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb;  in MyRHSMatMultTranspose()  local
103       vcb       = x[j][i].v;  in MyRHSMatMultTranspose()
104       vxx       = (-2.0 * vcb + x[j][i - 1].v + x[j][i + 1].v) * sx;  in MyRHSMatMultTranspose()
105       vyy       = (-2.0 * vcb + x[j - 1][i].v + x[j + 1][i].v) * sy;  in MyRHSMatMultTranspose()
106       y[j][i].u = appctx->D1 * (uxx + uyy) - ucb * (vc * vc + appctx->gamma) + vc * vc * vcb;  in MyRHSMatMultTranspose()
107 …tx->D2 * (vxx + vyy) - 2.0 * uc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appctx->kappa) * vcb;  in MyRHSMatMultTranspose()
124   PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb;  in MyIMatMultTranspose()  local
162       vcb       = x[j][i].v;  in MyIMatMultTranspose()
163       vxx       = (-2.0 * vcb + x[j][i - 1].v + x[j][i + 1].v) * sx;  in MyIMatMultTranspose()
164       vyy       = (-2.0 * vcb + x[j - 1][i].v + x[j + 1][i].v) * sy;  in MyIMatMultTranspose()
165       y[j][i].u = appctx->D1 * (uxx + uyy) - ucb * (vc * vc + appctx->gamma) + vc * vc * vcb;  in MyIMatMultTranspose()
166 …tx->D2 * (vxx + vyy) - 2.0 * uc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appctx->kappa) * vcb;  in MyIMatMultTranspose()
168       y[j][i].v = mctx->shift * vcb - y[j][i].v;  in MyIMatMultTranspose()
185   PetscScalar uc, uxx, uyy, vc, vxx, vyy, ucb, vcb;  in MyIMatMult()  local
223       vcb       = x[j][i].v;  in MyIMatMult()
224       vxx       = (-2.0 * vcb + x[j][i - 1].v + x[j][i + 1].v) * sx;  in MyIMatMult()
225       vyy       = (-2.0 * vcb + x[j - 1][i].v + x[j + 1][i].v) * 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 …= appctx->D2 * (vxx + vyy) + vc * vc * ucb + (2.0 * uc * vc - appctx->gamma - appctx->kappa) * vcb;  in MyIMatMult()
229       y[j][i].v = mctx->shift * vcb - y[j][i].v;  in MyIMatMult()