olb::particles::sorting Namespace Reference

Functions
template<typename T , typename PARTICLETYPE , typename ... NESTED_FIELDS>
std::size_t	partitionParticleSystem (ParticleSystem< T, PARTICLETYPE > &particleSystem)

Function Documentation

partitionParticleSystem()

template<typename T , typename PARTICLETYPE , typename ... NESTED_FIELDS>

std::size_t olb::particles::sorting::partitionParticleSystem

(

ParticleSystem< T, PARTICLETYPE > &

particleSystem

)

Definition at line 394 of file particleUtilities.h.

                                                                                   {
  using namespace descriptors;
  //Find first occurence of element belonging to partition B
  std::size_t iPfirstB=particleSystem.size(); //Fallback (simple to check)
  bool valB = true;
  for (std::size_t iP=0; iP<particleSystem.size(); ++iP){
    auto particle = particleSystem.get(iP);
    bool active = particle.template getField<NESTED_FIELDS...>();
    if (active==valB){
      iPfirstB=iP;
      break;
    }
  }
  //Find succeeding elements belonging to partition A and move them to othe beginning
  for (std::size_t iP=iPfirstB; iP<particleSystem.size(); ++iP){
    auto particle = particleSystem.get(iP);
    bool active = particle.template getField<NESTED_FIELDS...>();
    if (active!=valB){
      particleSystem.swapParticles(iP,iPfirstB);
      ++iPfirstB;
    }
  }
  //Return splitpoint between partitions
  return iPfirstB;
}

References olb::particles::ParticleSystem< T, PARTICLETYPE >::get(), partitionParticleSystem(), olb::particles::ParticleSystem< T, PARTICLETYPE >::size(), and olb::particles::ParticleSystem< T, PARTICLETYPE >::swapParticles().

Here is the call graph for this function:

Here is the caller graph for this function: