459   PetscInt      dim, xcStart, xcEnd, vcStart, vcEnd, Ns, Np, Npc, debug;  in InitializeParticles_Centroid()  local
485   Npc = (vcEnd - vcStart) * Ns;  in InitializeParticles_Centroid()
486   Np  = (xcEnd - xcStart) * Npc;  in InitializeParticles_Centroid()
504       for (PetscInt q = 0; q < Npc / Ns; ++q, ++p) {  in InitializeParticles_Centroid()
522     PetscInt      *pidx, Npc;  in InitializeParticles_Centroid()  local
525     PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx));  in InitializeParticles_Centroid()
528       for (PetscInt q = 0; q < Npc / Ns; ++q) {  in InitializeParticles_Centroid()
555     PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx));  in InitializeParticles_Centroid()
614     PetscInt          *pidx, Npc;  in InitializeWeights()  local
621     PetscCall(DMSwarmSortGetPointsPerCell(sw, c, &Npc, &pidx));  in InitializeWeights()
622 …Npc == (vcEnd - vcStart) * Ns, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Number of particles %" Petsc…  in InitializeWeights()
676     PetscCall(DMSwarmSortRestorePointsPerCell(sw, c, &Npc, &pidx));  in InitializeWeights()