olb::particles::statistics Namespace Reference

Functions
template<typename T , typename S , typename PARTICLETYPE , typename F >
void	evaluateParticleSystemStatistics (SuperParticleSystem< T, PARTICLETYPE > &sParticleSystem, std::vector< S > &data, int noSampledQuantities, F f, int root=0)
	Evaluate particle statistics for each cuboid (each ParticleSystem)
template<typename T , typename PARTICLETYPE , typename PCONDITION = conditions::valid_particles>
std::vector< std::size_t >	gatherActivity (SuperParticleSystem< T, PARTICLETYPE > &sParticleSystem)
	Gather number of active particles and total number for each cuboid.
template<typename T , typename PARTICLETYPE >
void	printActivityGathered (SuperParticleSystem< T, PARTICLETYPE > &sParticleSystem)
	Print activity gathered a main processor.

Function Documentation

evaluateParticleSystemStatistics()

template<typename T , typename S , typename PARTICLETYPE , typename F >

void olb::particles::statistics::evaluateParticleSystemStatistics	(	SuperParticleSystem< T, PARTICLETYPE > &	sParticleSystem,
		std::vector< S > &	data,
		int	noSampledQuantities,
		F	f,
		int	root = 0 )

Evaluate particle statistics for each cuboid (each ParticleSystem)

Definition at line 38 of file particleStatistics.h.

{
  //Set number of quantities (including globiC at first position)
  int noQ = noSampledQuantities + 1;
  //Retrieve load balancer and set up samples vector
  auto& loadBalancer = sParticleSystem.getSuperStructure().getLoadBalancer();
  std::vector<S> samples(loadBalancer.size()*noQ);
  //Iterate over particle systems
  communication::forSystemsInSuperParticleSystem( sParticleSystem,
    [&](ParticleSystem<T,PARTICLETYPE>& particleSystem, int iC, int globiC){
    //Set index for samples array
    int idx = iC*noQ;
    samples[idx] = globiC;
    //Perform requested sampling (assuming globiC at first pos)
    f ( particleSystem, globiC, idx+1, samples );
  });
 
  //Execute gathering
  const int numTasks = singleton::mpi().getSize();
  const int recvBufSize = samples.size()*numTasks;
  data.resize(recvBufSize);
  int sendCount = samples.size();
  int recvCounts[numTasks];
  int displs[numTasks];
  std::fill(recvCounts,&recvCounts[numTasks],sendCount);
  std::generate(displs,&displs[numTasks],[&sendCount,n = 0] () mutable { return sendCount*n++; });
#ifdef PARALLEL_MODE_MPI
  singleton::mpi().gatherv<S>(samples.data(), sendCount, data.data(),
                recvCounts, displs, root);
#endif
}

References olb::particles::communication::forSystemsInSuperParticleSystem(), olb::singleton::MpiManager::gatherv(), olb::SuperStructure< T, D >::getLoadBalancer(), olb::singleton::MpiManager::getSize(), olb::particles::SuperParticleSystem< T, PARTICLETYPE >::getSuperStructure(), and olb::singleton::mpi().

Here is the call graph for this function:

Here is the caller graph for this function:

gatherActivity()

template<typename T , typename PARTICLETYPE , typename PCONDITION = conditions::valid_particles>

std::vector< std::size_t > olb::particles::statistics::gatherActivity

(

SuperParticleSystem< T, PARTICLETYPE > &

sParticleSystem

)

Gather number of active particles and total number for each cuboid.

Definition at line 73 of file particleStatistics.h.

{
  //Set sampling properties
  unsigned noSampledQuantities = 2;
  std::vector<std::size_t> dataV;
  //Evaluate particle system statistics
  evaluateParticleSystemStatistics( sParticleSystem, dataV, noSampledQuantities,
    [](ParticleSystem<T,PARTICLETYPE>& particleSystem, int globiC, int idx, auto& samples){
    //Introduce couters for sampled quantities
    std::size_t noP = 0;
    std::size_t noActive = 0;
    //Iterate over particles
    forParticlesInParticleSystem<T,PARTICLETYPE,PCONDITION>( particleSystem,
      [&](Particle<T,PARTICLETYPE>& particle){
      //Evaluate counter increment
      if ( access::isActive(particle) ){ ++noActive; }
      ++noP;
    },globiC); //forParticlesInParticleSystem<T,PARTICLETYPE,PCONDITION>
   //Update samples vector with counters
   samples[idx+0] = noActive;
   samples[idx+1] = noP;
  });
  //Return gathered dataV
  return dataV;
}

References evaluateParticleSystemStatistics(), and olb::particles::access::isActive().

Here is the call graph for this function:

Here is the caller graph for this function:

printActivityGathered()

template<typename T , typename PARTICLETYPE >

void olb::particles::statistics::printActivityGathered

(

SuperParticleSystem< T, PARTICLETYPE > &

sParticleSystem

)

Print activity gathered a main processor.

Definition at line 103 of file particleStatistics.h.

                                                                                  {
  OstreamManager clout ( std::cout, "Activity" );
  auto activityStatistics = gatherActivity( sParticleSystem );
  unsigned noQ = 3;
  for (int idx=0; idx<activityStatistics.size(); idx+=noQ ){
    clout << "ParticleSystem " << activityStatistics[idx] << ":"
          << " active=" << activityStatistics[idx+1]
          << "/" << activityStatistics[idx+2]
          << std::endl;
  }
}

References gatherActivity().

Here is the call graph for this function: