Changeset 2d3854 for src/units/particle/interface_particles.cpp

Timestamp:

Aug 30, 2012, 12:30:33 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

b57b9b

Parents:

185795

Message:

Remove average force term in order to avoid particle drift.

git-svn-id: ​https://svn.version.fz-juelich.de/scafacos/trunk@2000 5161e1c8-67bf-11de-9fd5-51895aff932f

File:

• src/units/particle/interface_particles.cpp (modified) (1 diff)

src/units/particle/interface_particles.cpp

@@ -216 +216 @@
       }
 
+  /* Remove average force term */
+  Vector average_force = 0.0;
+  for (std::list<Particle::Particle>::const_iterator iter=particles.begin(); iter!=particles.end(); ++iter)
+    average_force += iter->Charge() * iter->Field();
+  const vmg_int& npl = MG::GetFactory().GetObjectStorageVal<vmg_int>("PARTICLE_NUM_LOCAL");
+  const vmg_int num_particles_global = comm.GlobalSum(npl);
+  average_force /= num_particles_global;
+  comm.GlobalSumArray(average_force.vec(), 3);
+  for (std::list<Particle::Particle>::iterator iter=particles.begin(); iter!=particles.end(); ++iter)
+    iter->Field() -= average_force / iter->Charge();
+
   comm.CommParticlesBack(particles);