135   PetscInt           dim, d, Np, p;  in RHSFunction()  local
140   PetscCall(VecGetLocalSize(U, &Np));  in RHSFunction()
143   Np /= 2 * dim;  in RHSFunction()
144   for (p = 0; p < Np; ++p) {  in RHSFunction()
163   PetscInt    dim, d, Np, p, rStart;  in RHSJacobian()  local
168   PetscCall(VecGetLocalSize(U, &Np));  in RHSJacobian()
170   Np /= 2 * dim;  in RHSJacobian()
171   for (p = 0; p < Np; ++p) {  in RHSJacobian()
186   PetscInt           Np, p;  in RHSFunctionX()  local
189   PetscCall(VecGetLocalSize(Xres, &Np));  in RHSFunctionX()
192   for (p = 0; p < Np; ++p) xres[p] = v[p];  in RHSFunctionX()
203   PetscInt           Np, p;  in RHSFunctionV()  local
208   PetscCall(VecGetLocalSize(Vres, &Np));  in RHSFunctionV()
209   for (p = 0; p < Np; ++p) vres[p] = -PetscSqr(omega) * x[p];  in RHSFunctionV()
219   PetscInt    dim, d, Np, p, rStart;  in RHSJacobianS()  local
224   PetscCall(VecGetLocalSize(U, &Np));  in RHSJacobianS()
226   Np /= 2 * dim;  in RHSJacobianS()
227   for (p = 0; p < Np; ++p) {  in RHSJacobianS()
243   PetscInt           dim, Np, p;  in RHSObjectiveF()  local
249   PetscCall(VecGetLocalSize(U, &Np));  in RHSObjectiveF()
250   Np /= 2 * dim;  in RHSObjectiveF()
251   for (p = 0; p < Np; ++p) {  in RHSObjectiveF()
269   PetscInt           dim, d, Np, p;  in RHSFunctionG()  local
276   PetscCall(VecGetLocalSize(U, &Np));  in RHSFunctionG()
277   Np /= 2 * dim;  in RHSFunctionG()
278   for (p = 0; p < Np; ++p) {  in RHSFunctionG()
293   PetscInt dim, Np;  in CreateSolution()  local
298   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in CreateSolution()
301   PetscCall(VecSetSizes(u, 2 * Np * dim, PETSC_DECIDE));  in CreateSolution()
319     PetscInt dim, Np;  in SetProblem()  local
322     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
324     PetscCall(MatSetSizes(J, 2 * Np * dim, 2 * Np * dim, PETSC_DECIDE, PETSC_DECIDE));  in SetProblem()
336     PetscInt        dim, Np;  in SetProblem()  local
339     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
340     PetscCall(ISCreateStride(PETSC_COMM_SELF, Np, 0, 2, &istmp));  in SetProblem()
342     PetscCall(ISCreateBlock(PETSC_COMM_SELF, dim, Np, idx, PETSC_COPY_VALUES, &isx));  in SetProblem()
345     PetscCall(ISCreateStride(PETSC_COMM_SELF, Np, 1, 2, &istmp));  in SetProblem()
347     PetscCall(ISCreateBlock(PETSC_COMM_SELF, dim, Np, idx, PETSC_COPY_VALUES, &isv));  in SetProblem()
402   PetscInt           dim, d, Np, p;  in ComputeError()  local
412   PetscCall(VecGetLocalSize(U, &Np));  in ComputeError()
414   Np /= 2 * dim;  in ComputeError()
415   for (p = 0; p < Np; ++p) {  in ComputeError()
448   PetscInt           dim, Np, p;  in EnergyMonitor()  local
458     PetscCall(VecGetLocalSize(U, &Np));  in EnergyMonitor()
459     Np /= 2 * dim;  in EnergyMonitor()
461     for (p = 0; p < Np; ++p) {  in EnergyMonitor()