158   PetscInt     i, mstart, mend, xm, M;  in InitialConditions()  local
163   PetscCall(DMDAGetCorners(da, &mstart, 0, 0, &xm, 0, 0));  in InitialConditions()
166   mend = mstart + xm;  in InitialConditions()
183 …for (i = mstart; i < mend; i++) u[i] = PetscSinScalar(PETSC_PI * i * 6. * h) + 3. * PetscSinScalar…  in InitialConditions()
206   PetscInt     i, mstart, mend, xm, M;  in Solution()  local
211   PetscCall(DMDAGetCorners(da, &mstart, 0, 0, &xm, 0, 0));  in Solution()
214   mend = mstart + xm;  in Solution()
228 …for (i = mstart; i < mend; i++) u[i] = PetscSinScalar(sc1 * (PetscReal)i + appctx->a * PETSC_PI * …  in Solution()
260   PetscInt    mstart, mend;  in RHSMatrixHeat()  local
268   PetscCall(DMDAGetCorners(da, &mstart, 0, 0, &xm, 0, 0));  in RHSMatrixHeat()
270   mend = mstart + xm;  in RHSMatrixHeat()
282   if (!mstart) {  in RHSMatrixHeat()
286     PetscCall(MatSetValues(A, 1, &mstart, 3, idx, v, INSERT_VALUES));  in RHSMatrixHeat()
287     mstart++;  in RHSMatrixHeat()
302   for (i = mstart; i < mend; i++) {  in RHSMatrixHeat()
311   PetscCall(DMDAGetCorners(da, &mstart, 0, 0, &xm, 0, 0));  in RHSMatrixHeat()
312   mend = mstart + xm;  in RHSMatrixHeat()
317     if (!mstart) {  in RHSMatrixHeat()
320       PetscCall(MatSetValues(A, 1, &mstart, 2, idx, v, ADD_VALUES));  in RHSMatrixHeat()
321       mstart++;  in RHSMatrixHeat()
331     for (i = mstart; i < mend; i++) {  in RHSMatrixHeat()
340     if (!mstart) {  in RHSMatrixHeat()
343       PetscCall(MatSetValues(A, 1, &mstart, 2, idx, v, ADD_VALUES));  in RHSMatrixHeat()
344       mstart++;  in RHSMatrixHeat()
354     for (i = mstart; i < mend; i++) {  in RHSMatrixHeat()