Esempio n. 1
0
void 
OPRijVij::eventUpdate(const GlobalEvent& globEvent, const NEventData& SDat)
{
  BOOST_FOREACH(const PairEventData& pDat, SDat.L2partChanges)
    process2PED(rvdotacc[mapKey(globEvent.getType(), getClassKey(globEvent))],
		pDat);
}
Esempio n. 2
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]);
    }
}