235   PetscInt    dim, d, Np, p, q;  in ComputeFieldAtParticles_Coulomb()  local
242   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Coulomb()
245   for (p = 0; p < Np; ++p) {  in ComputeFieldAtParticles_Coulomb()
249     for (q = 0; q < Np; ++q) {  in ComputeFieldAtParticles_Coulomb()
271   PetscInt     dim, d, cStart, cEnd, c, Np;  in ComputeFieldAtParticles_Primal()  local
278   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Primal()
381   PetscInt        dim, d, cStart, cEnd, c, Np, pot_field = 1;  in ComputeFieldAtParticles_Mixed()  local
388   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles_Mixed()
491   PetscInt dim, Np;  in ComputeFieldAtParticles()  local
498   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in ComputeFieldAtParticles()
500   PetscCall(PetscArrayzero(E, Np * dim));  in ComputeFieldAtParticles()
528   PetscInt           dim, d, Np, p;  in RHSFunction()  local
536   PetscCall(VecGetLocalSize(U, &Np));  in RHSFunction()
542   Np /= 2 * dim;  in RHSFunction()
543   for (p = 0; p < Np; ++p) {  in RHSFunction()
571   PetscInt         dim, d, Np, p, rStart;  in RHSJacobian()  local
576   PetscCall(VecGetLocalSize(U, &Np));  in RHSJacobian()
580   Np /= 2 * dim;  in RHSJacobian()
581   for (p = 0; p < Np; ++p) {  in RHSJacobian()
604   PetscInt           Np, p, dim, d;  in RHSFunctionX()  local
609   PetscCall(VecGetLocalSize(Xres, &Np));  in RHSFunctionX()
610   Np /= dim;  in RHSFunctionX()
613   for (p = 0; p < Np; ++p) {  in RHSFunctionX()
629   PetscInt           Np, p, dim, d;  in RHSFunctionV()  local
637   PetscCall(VecGetLocalSize(Vres, &Np));  in RHSFunctionV()
644   Np /= dim;  in RHSFunctionV()
645   for (p = 0; p < Np; ++p) {  in RHSFunctionV()
665   PetscInt    dim, d, Np, p, rStart;  in RHSJacobianS()  local
670   PetscCall(VecGetLocalSize(U, &Np));  in RHSJacobianS()
672   Np /= 2 * dim;  in RHSJacobianS()
673   for (p = 0; p < Np; ++p) {  in RHSJacobianS()
689   PetscInt           dim, Np, cStart, cEnd;  in RHSObjectiveF()  local
698   PetscCall(VecGetLocalSize(U, &Np));  in RHSObjectiveF()
710   Np /= 2 * dim;  in RHSObjectiveF()
749   PetscInt           dim, d, Np, p;  in RHSFunctionG()  local
757   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in RHSFunctionG()
766   for (p = 0; p < Np; ++p) {  in RHSFunctionG()
787   PetscInt dim, Np;  in CreateSolution()  local
792   PetscCall(DMSwarmGetLocalSize(sw, &Np));  in CreateSolution()
795   PetscCall(VecSetSizes(u, 2 * Np * dim, PETSC_DECIDE));  in CreateSolution()
813     PetscInt dim, Np;  in SetProblem()  local
816     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
818     PetscCall(MatSetSizes(J, 2 * Np * dim, 2 * Np * dim, PETSC_DECIDE, PETSC_DECIDE));  in SetProblem()
831     PetscInt        dim, Np, rstart;  in SetProblem()  local
835     PetscCall(DMSwarmGetLocalSize(sw, &Np));  in SetProblem()
837     PetscCall(ISCreateStride(PETSC_COMM_WORLD, Np, (rstart / dim) + 0, 2, &istmp));  in SetProblem()
839     PetscCall(ISCreateBlock(PETSC_COMM_WORLD, dim, Np, idx, PETSC_COPY_VALUES, &isx));  in SetProblem()
842     PetscCall(ISCreateStride(PETSC_COMM_WORLD, Np, (rstart / dim) + 1, 2, &istmp));  in SetProblem()
844     PetscCall(ISCreateBlock(PETSC_COMM_WORLD, dim, Np, idx, PETSC_COPY_VALUES, &isv));  in SetProblem()
969   PetscInt           dim, Np, p;  in ComputeError()  local
979   PetscCall(VecGetLocalSize(U, &Np));  in ComputeError()
982   Np /= 2 * dim;  in ComputeError()
983   for (p = 0; p < Np; ++p) {  in ComputeError()
1022   PetscInt           dim, d, Np, p, q;  in EnergyMonitor()  local
1029     PetscCall(VecGetLocalSize(U, &Np));  in EnergyMonitor()
1030     Np /= 2 * dim;  in EnergyMonitor()
1034     for (p = 0; p < Np; ++p) {  in EnergyMonitor()
1043         for (q = p + 1; q < Np; ++q) {  in EnergyMonitor()