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