420   PetscInt           Np;  in CheckNonNegativeWeights()  local
425   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in CheckNonNegativeWeights()
426 …for (PetscInt p = 0; p < Np; ++p) PetscCheck(w[p] >= 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "P…  in CheckNonNegativeWeights()
466   PetscInt    dim, Np;  in computeVelocityFEMMoments()
471   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in computeVelocityFEMMoments()
473   for (PetscInt p = 0; p < Np; ++p) {  in computeVelocityFEMMoments()
520   PetscInt   *species, dim, Np, gNp;  in MonitorEField()  local
528   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in MonitorEField()
536   for (PetscInt p = 0; p < Np; ++p) {  in MonitorEField()
597   PetscInt   dim, Np;  in MonitorInitialConditions()  local
604     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in MonitorInitialConditions()
619     for (PetscInt p = 0; p < Np; ++p) {  in MonitorInitialConditions()
1161   PetscInt      dim, xcStart, xcEnd, vcStart, vcEnd, Ns, Np, Npc, debug;  in InitializeParticles_Centroid()  local
1188   Np  = (xcEnd - xcStart) * Npc;  in InitializeParticles_Centroid()
1189   PetscCall(DMSwarmSetLocalSizes(sw, Np, 0));  in InitializeParticles_Centroid()
1192     PetscCallMPI(MPIU_Allreduce(&Np, &gNp, 1, MPIU_INT, MPIU_SUM, comm));  in InitializeParticles_Centroid()
1394   PetscInt   Np, dim;  in InitializeConstants()  local
1399   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in InitializeConstants()
1403   for (PetscInt p = 0; p < Np; ++p) {  in InitializeConstants()
1523   PetscInt   *species, dim, Np, Ns;  in ComputeFieldAtParticles_Coulomb()  local
1530   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Coulomb()
1536   for (PetscInt p = 0; p < Np; ++p) {  in ComputeFieldAtParticles_Coulomb()
1541     for (PetscInt q = 0; q < Np; ++q) {  in ComputeFieldAtParticles_Coulomb()
1572   PetscInt   dim, cStart, cEnd, Np;  in ComputeFieldAtParticles_Primal()  local
1578   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Primal()
1713   PetscInt   dim, cStart, cEnd, Np;  in ComputeFieldAtParticles_Mixed()  local
1719   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Mixed()
1832   PetscInt   dim, Np;  in ComputeFieldAtParticles()  local
1838   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles()
1848   PetscCall(PetscArrayzero(E, Np * dim));  in ComputeFieldAtParticles()
1878   PetscInt           dim, d, Np, p;  in RHSFunction()  local
1885   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in RHSFunction()
1889   Np /= 2 * dim;  in RHSFunction()
1890   for (p = 0; p < Np; ++p) {  in RHSFunction()
1912   PetscInt         dim, d, Np, p, rStart;  in RHSJacobian()  local
1917   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in RHSJacobian()
1921   Np /= 2 * dim;  in RHSJacobian()
1922   for (p = 0; p < Np; ++p) {  in RHSJacobian()
1944   PetscInt           Np, p, d, dim;  in RHSFunctionX()  local
1950   PetscCall(VecGetLocalSize(Xres, &Np));  in RHSFunctionX()
1953   Np /= dim;  in RHSFunctionX()
1954   for (p = 0; p < Np; ++p) {  in RHSFunctionX()
1971   PetscInt           Np, p, dim, d;  in RHSFunctionV()  local
1984   PetscCall(VecGetLocalSize(Vres, &Np));  in RHSFunctionV()
1987   Np /= dim;  in RHSFunctionV()
1988   for (p = 0; p < Np; ++p) {  in RHSFunctionV()
2002     for (PetscInt p = 0; p < Np; ++p) {  in RHSFunctionV()
2019   PetscInt    dim, d, Np, p, rStart;  in RHSJacobianS()  local
2024   PetscCall(VecGetLocalSize(U, &Np));  in RHSJacobianS()
2026   Np /= 2 * dim;  in RHSJacobianS()
2027   for (p = 0; p < Np; ++p) {  in RHSJacobianS()
2044   PetscInt           dim, Np, cStart, cEnd;  in RHSObjectiveF()  local
2061   PetscCall(VecGetLocalSize(U, &Np));  in RHSObjectiveF()
2062   Np /= 2 * dim;  in RHSObjectiveF()
2092   PetscInt           dim, d, Np, p;  in RHSFunctionG()  local
2097   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in RHSFunctionG()
2109   for (p = 0; p < Np; ++p) {  in RHSFunctionG()
2127   PetscInt dim, Np;  in CreateSolution()  local
2132   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in CreateSolution()
2135   PetscCall(VecSetSizes(u, 2 * Np * dim, PETSC_DECIDE));  in CreateSolution()
2153     PetscInt dim, Np;  in SetProblem()  local
2156     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
2158     PetscCall(MatSetSizes(J, 2 * Np * dim, 2 * Np * dim, PETSC_DECIDE, PETSC_DECIDE));  in SetProblem()
2171     PetscInt        dim, Np, rstart;  in SetProblem()  local
2175     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
2177     PetscCall(ISCreateStride(PETSC_COMM_WORLD, Np, (rstart / dim) + 0, 2, &istmp));  in SetProblem()
2179     PetscCall(ISCreateBlock(PETSC_COMM_WORLD, dim, Np, idx, PETSC_COPY_VALUES, &isx));  in SetProblem()
2182     PetscCall(ISCreateStride(PETSC_COMM_WORLD, Np, (rstart / dim) + 1, 2, &istmp));  in SetProblem()
2184     PetscCall(ISCreateBlock(PETSC_COMM_WORLD, dim, Np, idx, PETSC_COPY_VALUES, &isv));  in SetProblem()
2234   PetscInt  Np;  in line()  local
2242   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in line()
2244   x[0] = ((high[0] - low[0]) / Np) * (p + 0.5);  in line()
2339       PetscInt   Np, dim;  in MigrateParticles()  local
2341       PetscCall(DMSwarmGetLocalSize(sw, &Np));  in MigrateParticles()
2346       for (PetscInt p = 0; p < Np; ++p) {  in MigrateParticles()