53   PetscInt xs, xe, gxs, gxe, dof, gdof
105   PetscCallA(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
106   dof = i2*(xe - xs + 1)
147     PetscCallA(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
148     dof = i2*(xe - xs + 1)
150     call SaveSolutionToDisk(da, X, gdof, xs, xe)
165   subroutine GetLayout(da, mx, xs, xe, gxs, gxe, ierr)  argument
170     PetscInt mx, xs, xe, gxs, gxe
174 …PetscCall(DMDAGetCorners(da, xs, PETSC_NULL_INTEGER, PETSC_NULL_INTEGER, xm, PETSC_NULL_INTEGER, P…
176     xs = xs + 1
178     xe = xs + xm - 1
182   subroutine FormIFunctionLocal(mx, xs, xe, gxs, gxe, x, xdot, f, a, k, s, ierr)  argument
184     PetscInt mx, xs, xe, gxs, gxe
185     PetscScalar x(2, xs:xe)
186     PetscScalar xdot(2, xs:xe)
187     PetscScalar f(2, xs:xe)
191     do i = xs, xe
209     PetscInt mx, xs, xe, gxs, gxe
213     PetscCall(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
220     PetscCall(FormIFunctionLocal(mx, xs, xe, gxs, gxe, xx, xxdot, ff, ctx%a, ctx%k, ctx%s, ierr))
227   subroutine FormRHSFunctionLocal(mx, xs, xe, gxs, gxe, t, x, f, a, k, s, ierr)  argument
229     PetscInt mx, xs, xe, gxs, gxe
231     PetscScalar x(2, gxs:gxe), f(2, xs:xe)
254     do i = xs, xe
271     do i = xs, xe
290     PetscInt mx, xs, xe, gxs, gxe
294     PetscCall(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
308     PetscCall(FormRHSFunctionLocal(mx, xs, xe, gxs, gxe, t, xx, ff, ctx%a, ctx%k, ctx%s, ierr))
332     PetscInt mx, xs, xe, gxs, gxe
338     PetscCall(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
343     do i = xs, xe
360   subroutine FormInitialSolutionLocal(mx, xs, xe, gxs, gxe, x, a, k, s, ierr)  argument
362     PetscInt mx, xs, xe, gxs, gxe
363     PetscScalar x(2, xs:xe)
371     do i = xs, xe
394     PetscInt mx, xs, xe, gxs, gxe
398     PetscCall(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
403     PetscCall(FormInitialSolutionLocal(mx, xs, xe, gxs, gxe, xx, ctx%a, ctx%k, ctx%s, ierr))
449     PetscInt mx, xs, xe, gxs, gxe
458     PetscCall(GetLayout(da, mx, xs, xe, gxs, gxe, ierr))
464     do i = xs, xe
486   subroutine SaveSolutionToDisk(da, X, gdof, xs, xe)  argument
492     PetscInt xs, xe, two
495     PetscScalar data2(2, xs:xe), data(gdof)
501     data2 = reshape(xx(gdof:gdof), (/two, xe - xs + 1/))