60   MCtx       *mctx;  in MyRHSMatMultTranspose()  local
70   PetscCall(MatShellGetContext(A_shell, &mctx));  in MyRHSMatMultTranspose()
71   appctx = mctx->appctx;  in MyRHSMatMultTranspose()
72   PetscCall(TSGetDM(mctx->ts, &da));  in MyRHSMatMultTranspose()
89   PetscCall(DMDAVecGetArrayRead(da, mctx->U, &u));  in MyRHSMatMultTranspose()
111   PetscCall(DMDAVecRestoreArrayRead(da, mctx->U, &u));  in MyRHSMatMultTranspose()
119   MCtx       *mctx;  in MyIMatMultTranspose()  local
129   PetscCall(MatShellGetContext(A_shell, &mctx));  in MyIMatMultTranspose()
130   appctx = mctx->appctx;  in MyIMatMultTranspose()
131   PetscCall(TSGetDM(mctx->ts, &da));  in MyIMatMultTranspose()
148   PetscCall(DMDAVecGetArrayRead(da, mctx->U, &u));  in MyIMatMultTranspose()
167       y[j][i].u = mctx->shift * ucb - y[j][i].u;  in MyIMatMultTranspose()
168       y[j][i].v = mctx->shift * vcb - y[j][i].v;  in MyIMatMultTranspose()
172   PetscCall(DMDAVecRestoreArrayRead(da, mctx->U, &u));  in MyIMatMultTranspose()
180   MCtx       *mctx;  in MyIMatMult()  local
190   PetscCall(MatShellGetContext(A_shell, &mctx));  in MyIMatMult()
191   appctx = mctx->appctx;  in MyIMatMult()
192   PetscCall(TSGetDM(mctx->ts, &da));  in MyIMatMult()
209   PetscCall(DMDAVecGetArrayRead(da, mctx->U, &u));  in MyIMatMult()
228       y[j][i].u = mctx->shift * ucb - y[j][i].u;  in MyIMatMult()
229       y[j][i].v = mctx->shift * vcb - y[j][i].v;  in MyIMatMult()
233   PetscCall(DMDAVecRestoreArrayRead(da, mctx->U, &u));  in MyIMatMult()
245   MCtx           mctx;  in main()  local
311     mctx.ts     = ts;  in main()
312     mctx.appctx = &appctx;  in main()
313     PetscCall(VecDuplicate(x, &mctx.U));  in main()
319 …PETSC_COMM_WORLD, dof * xm * ym, PETSC_DETERMINE, dof * Mx * My, dof * Mx * My, &mctx, &appctx.A));  in main()
379     PetscCall(VecDestroy(&mctx.U));  in main()
390   MCtx *mctx;  in RHSJacobianShell()  local
393   PetscCall(MatShellGetContext(A, &mctx));  in RHSJacobianShell()
394   PetscCall(VecCopy(U, mctx->U));  in RHSJacobianShell()
400   MCtx *mctx;  in IJacobianShell()  local
403   PetscCall(MatShellGetContext(A, &mctx));  in IJacobianShell()
404   PetscCall(VecCopy(U, mctx->U));  in IJacobianShell()
406   mctx->shift = a;  in IJacobianShell()