Changeset 6c0b3f for src/Fragmentation/Summation/SetValues/SamplingGrid.cpp

Timestamp:

Jun 8, 2016, 10:45:28 PM (9 years ago)

Author:

Frederik Heber <heber@…>

Children:

26f238

Parents:

6369bc

git-author:

Frederik Heber <heber@…> (05/26/16 09:46:06)

git-committer:

Frederik Heber <heber@…> (06/08/16 22:45:28)

Message:

SamplingGridUnitTest checks downsample on window.

• FIX: getLengthsOfGrid() asserted that total equals 2^{level which is only true if the window spans the whole grid. We now assert for equal or less.}
• Renamed getLengthsOfGrid -> getLengthsOfWindow.
• Extended downsample() tests to check for level difference of one to three.

File:

• src/Fragmentation/Summation/SetValues/SamplingGrid.cpp (modified) (4 diffs)

src/Fragmentation/Summation/SetValues/SamplingGrid.cpp

@@ -674 +674 @@
 };
 
-static void getLengthsOfGrid(
+static void getLengthsOfWindow(
     int _total[NDIM],
     const SamplingGrid &_grid)
@@ -688 +688 @@
           } else
             _total[index] = 0;
-          ASSERT (_total[index] == ::pow(2, _grid.level),
+          // we can only assert that its atmost the maximum number of grid points
+          ASSERT (_total[index] <= ::pow(2, _grid.level),
               "SamplingGrid::downsample() - total "+toString(_total[index])
-              +" does not match 2^level: "+toString(_grid.level));
+              +" is not equal or less than 2^level: "+toString(_grid.level));
         }
 }
@@ -788 +789 @@
     static_cast<SamplingGridProperties &>(instance) = other;
     instance.setWindowSize(other.begin_window, other.end_window);
-    if (_level >= other.level) {
+    ASSERT( _level <= other.level,
+        "SamplingGrid::downsample() - desired level "+toString(_level)
+        +" larger than level "+toString(other.level)+" of the given values.");
+    if (_level == other.level) {
       instance.sampled_grid = other.sampled_grid;
     } else {
@@ -801 +805 @@
       int length_d[3];
       int length_s[3];
-      getLengthsOfGrid(length_s, other);
+      getLengthsOfWindow(length_s, other);
       for (instance.level = other.level-1; instance.level >= _level; --instance.level) {
-        getLengthsOfGrid(length_d, instance);
+        getLengthsOfWindow(length_d, instance);
         // we always have an eighth of the number of sample points as we stop
+        ASSERT( sourcevalues->size() % 8 == 0,
+            "SamplingGrid::downsample() - at level "+toString( instance.level)
+            +" given grid points "+toString(sourcevalues->size())+" are not even numbered per axis anymore.");
         sampledvalues_t downsampled(sourcevalues->size()/(size_t)8, 0.);
         restrictFullWeight(downsampled, length_d, *sourcevalues, length_s);