99   PetscInt           dim, d, Np, p;  in RHSFunction()  local
106   PetscCall(VecGetLocalSize(U, &Np));  in RHSFunction()
109   Np /= 2 * dim;  in RHSFunction()
110   for (p = 0; p < Np; ++p) {  in RHSFunction()
135   PetscInt         dim, d, Np, p, rStart;  in RHSJacobian()  local
140   PetscCall(VecGetLocalSize(U, &Np));  in RHSJacobian()
144   Np /= 2 * dim;  in RHSJacobian()
145   for (p = 0; p < Np; ++p) {  in RHSJacobian()
168   PetscInt           Np, p, dim, d;  in RHSFunctionX()  local
173   PetscCall(VecGetLocalSize(Xres, &Np));  in RHSFunctionX()
176   Np /= dim;  in RHSFunctionX()
177   for (p = 0; p < Np; ++p)  in RHSFunctionX()
190   PetscInt           Np, p, dim, d;  in RHSFunctionV()  local
197   PetscCall(VecGetLocalSize(Vres, &Np));  in RHSFunctionV()
201   Np /= dim;  in RHSFunctionV()
202   for (p = 0; p < Np; ++p) {  in RHSFunctionV()
220   PetscInt dim, Np;  in CreateSolution()  local
225   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in CreateSolution()
228   PetscCall(VecSetSizes(u, 2 * Np * dim, PETSC_DECIDE));  in CreateSolution()
246     PetscInt dim, Np;  in SetProblem()  local
249     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
251     PetscCall(MatSetSizes(J, 2 * Np * dim, 2 * Np * dim, PETSC_DECIDE, PETSC_DECIDE));  in SetProblem()
264     PetscInt        dim, Np, rstart;  in SetProblem()  local
268     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
270     PetscCall(ISCreateStride(PETSC_COMM_WORLD, Np, (rstart / dim) + 0, 2, &istmp));  in SetProblem()
272     PetscCall(ISCreateBlock(PETSC_COMM_WORLD, dim, Np, idx, PETSC_COPY_VALUES, &isx));  in SetProblem()
275     PetscCall(ISCreateStride(PETSC_COMM_WORLD, Np, (rstart / dim) + 1, 2, &istmp));  in SetProblem()
277     PetscCall(ISCreateBlock(PETSC_COMM_WORLD, dim, Np, idx, PETSC_COPY_VALUES, &isv));  in SetProblem()
431   PetscInt           dim, Np, p;  in ComputeError()  local
441   PetscCall(VecGetLocalSize(U, &Np));  in ComputeError()
444   Np /= 2 * dim;  in ComputeError()
445   for (p = 0; p < Np; ++p) {  in ComputeError()
481   PetscInt           dim, Np, p;  in EnergyMonitor()  local
488     PetscCall(VecGetLocalSize(U, &Np));  in EnergyMonitor()
489     Np /= 2 * dim;  in EnergyMonitor()
491     for (p = 0; p < Np; ++p) {  in EnergyMonitor()