17 #include "Foundation/console.h"
34 console.
XDebug() <<
"making CBWallInteractionGroup \n";
36 m_k=I->getSpringConst();
40 this->m_inner_count=0;
47 console.
XDebug() <<
"calculating " << m_bonded_interactions.size() <<
" bonded wall forces\n" ;
48 console.
XDebug() <<
"calculating " << m_elastic_interactions.size() <<
" elastic wall forces\n" ;
52 it!=m_bonded_interactions.end();
57 it!=m_elastic_interactions.end();
73 console.
XDebug() <<
"CBWallInteractionGroup<T>::applyForce: F = " << F <<
"\n";
79 it!=m_bonded_interactions.end();
82 K+=it->getStiffness();
86 it!=m_elastic_interactions.end();
89 K+=it->getStiffness();
93 double K_global=this->m_comm->sum_all(K);
98 console.
XDebug() <<
"CBWallInteractionGroup<T>::applyForce: unitF = " << O_f <<
"\n";
105 iter!=m_bonded_interactions.end();
109 Vec3 f_i=iter->getForce();
110 F_local+=(f_i*O_f)*O_f;
116 iter!=m_elastic_interactions.end();
120 Vec3 f_i=iter->getForce();
121 F_local+=(f_i*O_f)*O_f;
126 double fgx=this->m_comm->sum_all(F_local.X());
127 double fgy=this->m_comm->sum_all(F_local.Y());
128 double fgz=this->m_comm->sum_all(F_local.Z());
132 d=((F+F_global)*O_f)/K_global;
134 <<
"CBWallInteractionGroup<T>::applyForce: iteration " << it <<
", d = " << fabs(d) <<
"\n";
138 <<
"CBWallInteractionGroup<T>::applyForce: moving wall by " << d*O_f <<
"\n";
139 this->m_wall->moveBy(d*O_f);
141 }
while((it<10)&&(fabs(d)>10e-6));
143 <<
"CBWallInteractionGroup<T>::applyForce: d = " << fabs(d)
144 <<
", num iterations = " << it <<
"\n";
156 console.
XDebug() <<
"CBWallInteractionGroup::Update()\n" ;
160 m_bonded_interactions.erase(m_bonded_interactions.begin(),m_bonded_interactions.end());
161 this->m_inner_count=0;
163 typename ParallelParticleArray<T>::ParticleListHandle plh=
165 for(
typename ParallelParticleArray<T>::ParticleListIterator iter=plh->begin();
168 if(((*iter)->getTag() & m_mask )== (m_tag & m_mask)){
169 bool iflag=PPA->
isInInner((*iter)->getPos());
171 this->m_inner_count+=(iflag ? 1 : 0);
176 m_elastic_interactions.erase(m_elastic_interactions.begin(),m_elastic_interactions.end());
178 for(
typename ParallelParticleArray<T>::ParticleListIterator iter=plh->begin();
181 if(((*iter)->getTag() & m_mask )!=(m_tag & m_mask )){
182 bool iflag=PPA->
isInInner((*iter)->getPos());
184 this->m_inner_count+=(iflag ? 1 : 0);
187 console.
XDebug() <<
"end CBWallInteractionGroup::Update()\n";
191 ostream& operator<<(ostream& ost,const CBWallInteractionGroup<T>& IG)
193 ost <<
"CBWallInteractionGroup" << endl << flush;
194 ost << *(IG.m_wall) << endl << flush;