17 #include "Foundation/console.h"
35 console.
XDebug() <<
"making CEWallInteractionGroup \n";
37 m_k=I->getSpringConst();
45 console.
XDebug() <<
"calculating " << m_interactions.size() <<
" elastic wall forces\n" ;
49 it != m_interactions.end();
60 console.
XDebug() <<
"CEWallInteractionGroup::Update()\n" ;
63 <<
"CEWallInteractionGroup::Update: wall origin = " << this->m_wall->getOrigin()
64 <<
", wall normal = " << this->m_wall->getNormal() <<
"\n" ;
68 m_interactions.erase(m_interactions.begin(),m_interactions.end());
69 this->m_inner_count=0;
71 typename ParallelParticleArray<T>::ParticleListHandle plh=
73 for(
typename ParallelParticleArray<T>::ParticleListIterator iter=plh->begin();
76 bool iflag=PPA->
isInInner((*iter)->getPos());
78 this->m_inner_count+=(iflag ? 1 : 0);
81 console.
XDebug() <<
"end CEWallInteractionGroup::Update()\n";
106 iter!=m_interactions.end();
108 k_local+=iter->getStiffness();
112 m_k_global=this->m_comm->sum_all(k_local);
119 iter!=m_interactions.end();
123 Vec3 f_i=iter->getForce();
124 F_local+=(f_i*O_f)*O_f;
130 double fgx=this->m_comm->sum_all(F_local.X());
131 double fgy=this->m_comm->sum_all(F_local.Y());
132 double fgz=this->m_comm->sum_all(F_local.Z());
143 this->m_wall->moveBy(d*O_f);
151 this->m_wall->moveBy(d*O_f);
154 }
while((it<50)&&(ef>1e-3));
159 ostream& operator<<(ostream& ost,const CEWallInteractionGroup<T>& IG)
161 ost <<
"CEWallInteractionGroup" << endl << flush;
162 ost << *(IG.m_wall) << endl << flush;