Changeset 5c22be for src/units/particle/interface_particles.cpp

Timestamp:

Sep 6, 2015, 10:23:11 PM (10 years ago)

Author:

Frederik Heber <heber@…>

Children:

ebc866

Parents:

b4e1be

git-author:

Frederik Heber <heber@…> (03/09/15 16:50:34)

git-committer:

Frederik Heber <heber@…> (09/06/15 22:23:11)

Message:

tempcommit: Small changes to interface_particles.

File:

• src/units/particle/interface_particles.cpp (modified) (3 diffs)

src/units/particle/interface_particles.cpp

@@ -71 +71 @@
   Grid& particle_grid = comm.GetParticleGrid();
 
+//  grid.Clear();
   particle_grid.Clear();
 
@@ -184 +185 @@
           // Subtract self-induced potential
           (*p1)->Pot() -= (*p1)->Charge() * spl.GetAntiDerivativeAtZero();
+//      spl.SubtractSelfInducedForces(particle_grid, **p1);
 
 #ifdef OUTPUT_DEBUG
@@ -233 +235 @@
   const vmg_int& num_particles_local = factory.GetObjectStorageVal<vmg_int>("PARTICLE_NUM_LOCAL");
   const vmg_float* p = factory.GetObjectStorageArray<vmg_float>("PARTICLE_POTENTIAL_ARRAY");
-
+  const vmg_float* f = factory.GetObjectStorageArray<vmg_float>("PARTICLE_FIELD_ARRAY");
+
+  // extract forces
+  if (num_particles_local != 0) {
+    size_t index = 0;
+    comm.PrintOnce(Debug, "%d force vector:        %e %e %e", (index/3)+1, f[index++], f[index++], f[index++]);
+  }
 
   e_long = comm.GlobalSumRoot(e_long);