Changeset 6e10c12

Timestamp:

Oct 10, 2015, 11:41:00 PM (10 years ago)

Author:

Frederik Heber <heber@…>

Children:

784b62

Parents:

199fe3

Message:

tempcommit: Removed setting net force to zero.

File:

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

src/units/particle/interface_particles.cpp

@@ -219 +219 @@
 
   /* 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();
+//  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);