65   PetscInt       dim, d, cStart, cEnd, c, Np, p;  in SetInitialCoordinates()  local
73   Np = user->particlesPerCell;  in SetInitialCoordinates()
85     if (Np == 1) {  in SetInitialCoordinates()
97       for (p = 0; p < Np; ++p) {  in SetInitialCoordinates()
98         const PetscInt n   = c * Np + p;  in SetInitialCoordinates()
126   PetscInt     dim, d, cStart, cEnd, c, Np, p, n;  in SetInitialConditions()  local
131   Np = user->particlesPerCell;  in SetInitialConditions()
135 …rt) * Np, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "TS solution local size %" PetscInt_FMT " != %" Pets…  in SetInitialConditions()
139     for (p = 0; p < Np; ++p) {  in SetInitialConditions()
140       const PetscInt n = c * Np + p;  in SetInitialConditions()
153   PetscInt      dim, cStart, cEnd, c, Np = user->particlesPerCell, p;  in CreateParticles()  local
169   PetscCall(DMSwarmSetLocalSizes(*sw, (cEnd - cStart) * Np, 0));  in CreateParticles()
173     for (p = 0; p < Np; ++p) {  in CreateParticles()
174       const PetscInt n = c * Np + p;  in CreateParticles()
280   PetscInt           dim, d, Np, Ncp, p, cStart, cEnd, c;  in RHSFunctionParticles()  local
284   PetscCall(VecGetLocalSize(U, &Np));  in RHSFunctionParticles()
293   Np /= dim;  in RHSFunctionParticles()
294   Ncp = Np / (cEnd - cStart);  in RHSFunctionParticles()
297   for (p = 0; p < Np; ++p) {  in RHSFunctionParticles()
328   for (p = 0; p < Np; ++p) {  in RHSFunctionParticles()
356   PetscInt           dim, Np, p;  in HGMonitor()  local
371     PetscCall(VecGetLocalSize(U, &Np));  in HGMonitor()
372     Np /= dim;  in HGMonitor()
375     for (p = 0; p < Np; ++p) PetscCall(PetscDrawHGAddValue(user->drawhg, u[p]));  in HGMonitor()
387   PetscInt           Np, p;  in SPMonitor()  local
399     PetscCall(VecGetLocalSize(U, &Np));  in SPMonitor()
402     PetscCall(PetscMalloc1(Np, &v));  in SPMonitor()
403     for (p = 0; p < Np; ++p) v[p] = PetscRealPart(u[p]);  in SPMonitor()
405     for (p = 0; p < Np - 1; ++p) PetscCall(PetscDrawSPAddPoint(user->drawsp, &coords[p], &v[p]));  in SPMonitor()
420   PetscInt           dim, Np, p, cStart, cEnd;  in KSConv()  local
431     PetscCall(VecGetLocalSize(U, &Np));  in KSConv()
434     PetscCall(PetscMalloc1(Np, &v));  in KSConv()
438     for (p = 0; p < Np; ++p) v[p] = PetscRealPart(u[p]);  in KSConv()
441     for (p = 0; p < Np; ++p) {  in KSConv()
449     for (p = 0; p < Np; ++p) {  in KSConv()