13 template<
typename P,
typename I>
21 template<
typename P,
typename InteractionType>
30 template<
typename T,
typename InteractionType>
33 console.
XDebug() <<
"ParallelInteractionStorage_NE::Update\n";
37 if(m_update_timestamp != this->m_ppa->getTimeStamp()){
39 this->m_interactions.clear();
40 m_set.erase(m_set.begin(),m_set.end());
42 typename ParallelParticleArray<T>::PairListHandle plh =
45 for(
typename ParallelParticleArray<T>::PairListIterator iter=plh->begin();
50 tv.push_back(iter->first->getID());
51 tv.push_back(iter->second->getID());
53 if(!m_exIG->isIn(tv)){
54 this->m_interactions.push_back(InteractionType(iter->first,iter->second,m_param));
55 m_set.insert(pair<int,int>(iter->first->getID(),iter->second->getID()));
59 this->m_interactions.push_back(InteractionType(iter->first,iter->second,m_param));
60 m_set.insert(pair<int,int>(iter->first->getID(),iter->second->getID()));
66 typename ParallelParticleArray<T>::PairListHandle plh =
69 for(
typename ParallelParticleArray<T>::PairListIterator iter=plh->begin();
75 tv.push_back(iter->first->getID());
76 tv.push_back(iter->second->getID());
78 if(!m_exIG->isIn(tv)){
79 this->m_interactions.push_back(InteractionType(iter->first,iter->second,m_param));
80 m_set.insert(pair<int,int>(iter->first->getID(),iter->second->getID()));
84 this->m_interactions.push_back(InteractionType(iter->first,iter->second,m_param));
85 m_set.insert(pair<int,int>(iter->first->getID(),iter->second->getID()));
90 m_update_timestamp = this->m_ppa->getTimeStamp();
92 console.
XDebug() <<
"added " << count_l <<
" pairs to EIG\n";
93 console.
XDebug() <<
"end ParallelInteractionStorage_NE::Update\n";
102 template<
typename P,
typename InteractionType>
107 res=m_set.find(make_pair(pids[0],pids[1]))!=m_set.end();
115 template<
typename P,
typename InteractionType>
120 << this->m_interactions.size()
121 <<
" interaction forces\n";
124 typename list<InteractionType>::iterator it = this->m_interactions.begin();
125 it != this->m_interactions.end();
132 template<
typename P,
typename InteractionType>
137 << this->m_interactions.size()
138 <<
" interaction heat transfers\n" ;
141 typename list<InteractionType>::iterator it = this->m_interactions.begin();
142 it != this->m_interactions.end();