Example #1
0
  void
  OPMisc::eventUpdate(const IntEvent& eevent, const PairEventData& PDat)
  {
    stream(eevent.getdt());
    eventUpdate(PDat);

    newEvent(eevent.getParticle1ID(), eevent.getType(), getClassKey(eevent));
    newEvent(eevent.getParticle2ID(), eevent.getType(), getClassKey(eevent));
  }
Example #2
0
  void 
  OPCollMatrix::eventUpdate(const IntEvent& iEvent, const PairEventData&)
  {
    newEvent(iEvent.getParticle1ID(), iEvent.getType(), 
	     getClassKey(iEvent));

    newEvent(iEvent.getParticle2ID(), iEvent.getType(), 
	     getClassKey(iEvent));
  }
Example #3
0
  void 
  OPCollMatrix::eventUpdate(const System& sysEvent, const NEventData& SDat, const double&)
  {
    for (const ParticleEventData& pData : SDat.L1partChanges)
      newEvent(pData.getParticleID(), pData.getType(), getClassKey(sysEvent));
  
    for (const PairEventData& pData : SDat.L2partChanges)
      {
	newEvent(pData.particle1_.getParticleID(), 
		 pData.getType(), getClassKey(sysEvent));

	newEvent(pData.particle2_.getParticleID(), 
		 pData.getType(), getClassKey(sysEvent));  
      } 
  }
Example #4
0
  void
  OPMisc::eventUpdate(const System& eevent, const NEventData& NDat, const double& dt)
  {
    stream(dt);
    eventUpdate(NDat);

    BOOST_FOREACH(const ParticleEventData& pData, NDat.L1partChanges)
      newEvent(pData.getParticleID(), pData.getType(), getClassKey(eevent));

    BOOST_FOREACH(const PairEventData& pData, NDat.L2partChanges)
      {
	newEvent(pData.particle1_.getParticleID(), pData.getType(), getClassKey(eevent));
	newEvent(pData.particle2_.getParticleID(), pData.getType(), getClassKey(eevent));
      }
  }
Example #5
0
  void 
  OPCollMatrix::eventUpdate(const LocalEvent& localEvent, const NEventData& SDat)
  {
    for (const ParticleEventData& pData : SDat.L1partChanges)
      newEvent(pData.getParticleID(), 
	       pData.getType(), getClassKey(localEvent));  
  
    for (const PairEventData& pData : SDat.L2partChanges)
      {
	newEvent(pData.particle1_.getParticleID(), 
		 pData.getType(), getClassKey(localEvent));

	newEvent(pData.particle2_.getParticleID(),
		 pData.getType(), getClassKey(localEvent));
      }
  }
Example #6
0
void
OPRijVij::eventUpdate(const System& sysEvent, const NEventData& SDat, const double&)
{
  BOOST_FOREACH(const PairEventData& pDat, SDat.L2partChanges)
    process2PED(rvdotacc[mapKey(sysEvent.getType(), getClassKey(sysEvent))],
		pDat);
}
Example #7
0
void 
OPRijVij::eventUpdate(const LocalEvent& localEvent, const NEventData& SDat)
{
  BOOST_FOREACH(const PairEventData& pDat, SDat.L2partChanges)
    process2PED(rvdotacc[mapKey(localEvent.getType(), getClassKey(localEvent))],
		pDat);
}
Example #8
0
void 
OPRijVij::eventUpdate(const IntEvent& iEvent, const PairEventData& pDat)
{
  
  process2PED(rvdotacc[mapKey(iEvent.getType(), 
			      getClassKey(iEvent))], pDat);
}
Example #9
0
void 
OPCollMatrix::eventUpdate(const GlobalEvent& globEvent, const NEventData& SDat)
{
  BOOST_FOREACH(const ParticleEventData& pData, SDat.L1partChanges)
    newEvent(pData.getParticle().getID(), pData.getType(), 
	     getClassKey(globEvent));  
  
  BOOST_FOREACH(const PairEventData& pData, SDat.L2partChanges)
    {
      newEvent(pData.particle1_.getParticle().getID(), 
	       pData.getType(), getClassKey(globEvent));

      newEvent(pData.particle2_.getParticle().getID(), 
	       pData.getType(), getClassKey(globEvent));
    }
}
Example #10
0
void 
OPCubeComp::eventUpdate(const IntEvent& iEvent, const PairEventData& pDat)
{
  mapdata& ref = angles[mapKey(iEvent.getType(), getClassKey(iEvent))];

  std::vector<double> vals(NDIM, 0);

  for (size_t i(0); i < NDIM; ++i)
    vals[i] = pDat.rij[i] * pDat.rij[i] / Sim->dynamics.units().unitArea();

  std::sort(vals.begin(), vals.end());
  
  for (size_t i(0); i < NDIM; ++i)
    ref.angles[i].addVal(vals[i]);
}
Example #11
0
void 
OPCubeComp::eventUpdate(const GlobalEvent& globEvent, const NEventData& SDat)
{
  BOOST_FOREACH(const PairEventData& pDat, SDat.L2partChanges)
    {
      mapdata& ref = angles[mapKey(globEvent.getType(), 
				   getClassKey(globEvent))];
      std::vector<double> vals(NDIM, 0);
      
      for (size_t i(0); i < NDIM; ++i)
	vals[i] = pDat.rij[i] * pDat.rij[i] / Sim->dynamics.units().unitArea();
      
      std::sort(vals.begin(), vals.end());
      
      for (size_t i(0); i < NDIM; ++i)
	ref.angles[i].addVal(vals[i]);
    }
}