Changeset a150d0 for src/interface/interface_particles.cpp

Timestamp:

Apr 25, 2012, 4:21:08 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

7dd744

Parents:

b51c3b

Message:

vmg: Added parameter to control the interpolation degree (Tensor Lagrange interpolation) for interpolating the potentials from the grid to the particle positions.

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

File:

• src/interface/interface_particles.cpp (modified) (2 diffs)

src/interface/interface_particles.cpp

@@ -102 +102 @@
   vmg_float length;
   Vector dist_vec;
-  InterpolatePolynomial ip(8);
 
 #ifdef DEBUG_OUTPUT
@@ -119 +118 @@
    */
   const vmg_int& num_particles_local = factory.GetObjectStorageVal<vmg_int>("PARTICLE_NUM_LOCAL");
-  const int& near_field_cells = factory.GetObjectStorageVal<int>("PARTICLE_NEAR_FIELD_CELLS");
+  const vmg_int& near_field_cells = factory.GetObjectStorageVal<int>("PARTICLE_NEAR_FIELD_CELLS");
+  const vmg_int& interpolation_degree = factory.GetObjectStorageVal<int>("PARTICLE_INTERPOLATION_DEGREE");
   vmg_float* p = factory.GetObjectStorageArray<vmg_float>("PARTICLE_POTENTIAL_ARRAY");
   vmg_float* f = factory.GetObjectStorageArray<vmg_float>("PARTICLE_FIELD_ARRAY");
+
+  InterpolatePolynomial ip(interpolation_degree);
 
   const vmg_float r_cut = (near_field_cells+0.5) * grid.Extent().MeshWidth().Max();