582   PetscInt           Np;  in CheckNonNegativeWeights()  local
587   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in CheckNonNegativeWeights()
588 …for (PetscInt p = 0; p < Np; ++p) PetscCheck(w[p] >= 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "P…  in CheckNonNegativeWeights()
754   PetscInt   *species, dim, Np, gNp;  in MonitorEField()  local
764   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in MonitorEField()
772   for (PetscInt p = 0; p < Np; ++p) {  in MonitorEField()
1113   PetscInt   dim, Np;  in MonitorInitialConditions()  local
1120     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in MonitorInitialConditions()
1135     for (PetscInt p = 0; p < Np; ++p) {  in MonitorInitialConditions()
1775   PetscInt      dim, xcStart, xcEnd, vcStart, vcEnd, Ns, Np, Npc, debug;  in InitializeParticles_Centroid()  local
1802   Np  = (xcEnd - xcStart) * Npc;  in InitializeParticles_Centroid()
1803   PetscCall(DMSwarmSetLocalSizes(sw, Np, 0));  in InitializeParticles_Centroid()
1806     PetscCallMPI(MPIU_Allreduce(&Np, &gNp, 1, MPIU_INT, MPIU_SUM, comm));  in InitializeParticles_Centroid()
2008   PetscInt   Np, dim;  in InitializeConstants()  local
2013   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in InitializeConstants()
2017   for (PetscInt p = 0; p < Np; ++p) {  in InitializeConstants()
2147   PetscInt   *species, dim, Np, Ns;  in ComputeFieldAtParticles_Coulomb()  local
2154   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Coulomb()
2160   for (PetscInt p = 0; p < Np; ++p) {  in ComputeFieldAtParticles_Coulomb()
2165     for (PetscInt q = 0; q < Np; ++q) {  in ComputeFieldAtParticles_Coulomb()
2196   PetscInt   dim, cStart, cEnd, Np;  in ComputeFieldAtParticles_Primal()  local
2202   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Primal()
2317   PetscInt   dim, cStart, cEnd, Np;  in ComputeFieldAtParticles_Mixed()  local
2323   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Mixed()
2437   PetscInt   dim, Np;  in ComputeFieldAtParticles()  local
2443   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles()
2453   PetscCall(PetscArrayzero(E, Np * dim));  in ComputeFieldAtParticles()
2483   PetscInt           dim, d, Np, p;  in RHSFunction()  local
2490   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in RHSFunction()
2494   Np /= 2 * dim;  in RHSFunction()
2495   for (p = 0; p < Np; ++p) {  in RHSFunction()
2517   PetscInt         dim, d, Np, p, rStart;  in RHSJacobian()  local
2522   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in RHSJacobian()
2526   Np /= 2 * dim;  in RHSJacobian()
2527   for (p = 0; p < Np; ++p) {  in RHSJacobian()
2549   PetscInt           Np, p, d, dim;  in RHSFunctionX()  local
2555   PetscCall(VecGetLocalSize(Xres, &Np));  in RHSFunctionX()
2558   Np /= dim;  in RHSFunctionX()
2559   for (p = 0; p < Np; ++p) {  in RHSFunctionX()
2576   PetscInt           Np, p, dim, d;  in RHSFunctionV()  local
2589   PetscCall(VecGetLocalSize(Vres, &Np));  in RHSFunctionV()
2592   Np /= dim;  in RHSFunctionV()
2593   for (p = 0; p < Np; ++p) {  in RHSFunctionV()
2607     for (PetscInt p = 0; p < Np; ++p) {  in RHSFunctionV()
2624   PetscInt    dim, d, Np, p, rStart;  in RHSJacobianS()  local
2629   PetscCall(VecGetLocalSize(U, &Np));  in RHSJacobianS()
2631   Np /= 2 * dim;  in RHSJacobianS()
2632   for (p = 0; p < Np; ++p) {  in RHSJacobianS()
2649   PetscInt           dim, Np, cStart, cEnd;  in RHSObjectiveF()  local
2666   PetscCall(VecGetLocalSize(U, &Np));  in RHSObjectiveF()
2667   Np /= 2 * dim;  in RHSObjectiveF()
2697   PetscInt           dim, d, Np, p;  in RHSFunctionG()  local
2702   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in RHSFunctionG()
2714   for (p = 0; p < Np; ++p) {  in RHSFunctionG()
2732   PetscInt dim, Np;  in CreateSolution()  local
2737   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in CreateSolution()
2740   PetscCall(VecSetSizes(u, 2 * Np * dim, PETSC_DECIDE));  in CreateSolution()
2758     PetscInt dim, Np;  in SetProblem()  local
2761     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
2763     PetscCall(MatSetSizes(J, 2 * Np * dim, 2 * Np * dim, PETSC_DECIDE, PETSC_DECIDE));  in SetProblem()
2776     PetscInt        dim, Np, rstart;  in SetProblem()  local
2780     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
2782     PetscCall(ISCreateStride(PETSC_COMM_WORLD, Np, (rstart / dim) + 0, 2, &istmp));  in SetProblem()
2784     PetscCall(ISCreateBlock(PETSC_COMM_WORLD, dim, Np, idx, PETSC_COPY_VALUES, &isx));  in SetProblem()
2787     PetscCall(ISCreateStride(PETSC_COMM_WORLD, Np, (rstart / dim) + 1, 2, &istmp));  in SetProblem()
2789     PetscCall(ISCreateBlock(PETSC_COMM_WORLD, dim, Np, idx, PETSC_COPY_VALUES, &isv));  in SetProblem()
2839   PetscInt  Np;  in line()  local
2847   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in line()
2849   x[0] = ((high[0] - low[0]) / Np) * (p + 0.5);  in line()
2945       PetscInt   Np, dim;  in MigrateParticles()  local
2947       PetscCall(DMSwarmGetLocalSize(sw, &Np));  in MigrateParticles()
2952       for (PetscInt p = 0; p < Np; ++p) {  in MigrateParticles()