68   PetscInt       dim, d, cStart, cEnd, c, Np, p;  in SetInitialCoordinates()  local
79   Np = user->N;  in SetInitialCoordinates()
92     if (Np == 1) {  in SetInitialCoordinates()
98       for (p = 0; p < Np; ++p) {  in SetInitialCoordinates()
99         const PetscInt n   = c * Np + p;  in SetInitialCoordinates()
115     for (p = 0; p < Np; ++p) {  in SetInitialCoordinates()
140   PetscInt     dim, d, cStart, cEnd, c, Np, p, n;  in SetInitialConditions()  local
145   Np = user->N;  in SetInitialConditions()
152     for (p = 0; p < Np; ++p) {  in SetInitialConditions()
153       const PetscInt n = c * Np + p;  in SetInitialConditions()
166   PetscInt      dim, cStart, cEnd, c, Np = user->N, p;  in CreateParticles()  local
189   PetscCall(DMSwarmSetLocalSizes(*sw, (cEnd - cStart) * Np, 0));  in CreateParticles()
193     for (p = 0; p < Np; ++p) {  in CreateParticles()
194       const PetscInt n = c * Np + p;  in CreateParticles()
287 static PetscErrorCode ComputeGradS(PetscInt dim, PetscInt Np, const PetscReal vp[], const PetscReal…  in ComputeGradS()  argument
304         for (q = 0; q < Np; ++q) sum += Gaussian(dim, &velocity[q * dim], epsilon, vc_l);  in ComputeGradS()
340   PetscInt           dim, Np, p, q;  in RHSFunctionParticles()  local
346   PetscCall(VecGetLocalSize(U, &Np));  in RHSFunctionParticles()
351   Np /= dim;  in RHSFunctionParticles()
353   for (p = 0; p < Np; ++p) {  in RHSFunctionParticles()
356     PetscCall(ComputeGradS(dim, Np, &velocity[p * dim], velocity, gradS_p, user));  in RHSFunctionParticles()
357     for (q = 0; q < Np; ++q) {  in RHSFunctionParticles()
361       PetscCall(ComputeGradS(dim, Np, &velocity[q * dim], velocity, gradS_q, user));  in RHSFunctionParticles()