Changeset ab63b6 for src/interface/interface_particles.cpp

Timestamp:

Dec 8, 2011, 12:53:48 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

138f86

Parents:

facba0

Message:

Ensure that interpolating splines are normalized wrt the grid.

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

File:

• src/interface/interface_particles.cpp (modified) (1 diff)

src/interface/interface_particles.cpp

@@ -100 +100 @@
 #endif
 
-  for (iter=particles.begin(); iter!=particles.end(); ++iter) {
-
-    // Compute global and local grid index of the lower left corner of the grid cell containing the particle
-    index_global = static_cast<Index>((iter->Pos() - temp_grid->Extent().Begin()) / temp_grid->Extent().MeshWidth());
-    index_local = index_global - temp_grid->Global().BeginLocal() + temp_grid->Local().Begin();
-
-    // Iterate over all grid points which lie in the support of the interpolating B-Spline
-    for (index.X()=-1*near_field_cells-1; index.X()<=near_field_cells+1; ++index.X())
-      for (index.Y()=-1*near_field_cells-1; index.Y()<=near_field_cells+1; ++index.Y())
-        for (index.Z()=-1*near_field_cells-1; index.Z()<=near_field_cells+1; ++index.Z()) {
-
-          // Compute distance from grid point to particle
-          length = ( iter->Pos() - grid.Extent().Begin() - grid.Extent().MeshWidth() * (index_global+index) ).Length();
-
-          // If grid point lies in the support of the B-Spline, do charge interpolation
-          if (length < r_cut)
-            (*temp_grid)(index_local+index) += spl.EvaluateSpline(length) * iter->Charge();
-
-        }
-
-  }
+  for (iter=particles.begin(); iter!=particles.end(); ++iter)
+    spl.SetSpline(*temp_grid, *iter, near_field_cells);
 
   // Communicate charges over halo