Changeset 26f238 for src/Fragmentation/Summation/SetValues/SamplingGrid.cpp

Timestamp:

Jun 9, 2016, 5:55:43 AM (9 years ago)

Author:

Frederik Heber <heber@…>

Children:

52ecd0

Parents:

6c0b3f

git-author:

Frederik Heber <heber@…> (05/26/16 15:36:26)

git-committer:

Frederik Heber <heber@…> (06/09/16 05:55:43)

Message:

Added padWithZerosForEvenNumberedSamples(), required by downsample().

• added unit test function.

File:

• src/Fragmentation/Summation/SetValues/SamplingGrid.cpp (modified) (1 diff)

src/Fragmentation/Summation/SetValues/SamplingGrid.cpp

@@ -780 +780 @@
 }
 
+void SamplingGrid::padWithZerosForEvenNumberedSamples()
+{
+  size_t wbegin_index[NDIM];
+  size_t wend_index[NDIM];
+  size_t wlength_index[NDIM];
+  getDiscreteWindowIndices(wbegin_index, wlength_index, wend_index);
+
+  // calculate new window (always extend it such that both indices are even)
+  bool changed = false;
+  size_t wnewbegin_index[NDIM];
+  size_t wnewend_index[NDIM];
+  for(size_t i=0;i<NDIM;++i) {
+    // We require begin and end of window on even indices (and inside grid).
+    wnewbegin_index[i] = wbegin_index[i];
+    if ((wnewbegin_index[i] % (size_t)2) != 0) {
+      if (wnewbegin_index[i] > 0)
+        --wnewbegin_index[i];
+      else
+        wnewbegin_index[i] = 0;
+      changed = true;
+    }
+    wnewend_index[i] = wend_index[i];
+    if ((wnewend_index[i] % (size_t)2) != 0) {
+      if (wnewend_index[i] < getGridPointsPerAxis())
+        ++wnewend_index[i];
+      else
+        wnewend_index[i] = getGridPointsPerAxis();
+      changed = true;
+    }
+    ASSERT( (wbegin_index[i] >= 0) && (wend_index[i] <= getGridPointsPerAxis()),
+        "SamplingGrid::padWithZerosForEvenNumberedSamples() - indices "
+        +toString(wbegin_index[i])+" and "+toString(wend_index[i])+" larger than grid "
+        +toString(getGridPointsPerAxis())+".");
+  }
+  if (changed) {
+    double begin_newwindow[NDIM];
+    double end_newwindow[NDIM];
+    for(size_t i=0;i<NDIM;++i) {
+      const double delta = getDeltaPerAxis(i);
+      begin_newwindow[i] = begin[i]+delta*wnewbegin_index[i];
+      end_newwindow[i] = begin[i]+delta*wnewend_index[i];
+    }
+    // extend window
+    extendWindow(begin_newwindow, end_newwindow);
+  }
+  ASSERT( getWindowGridPoints() % (size_t)8 == 0,
+      "SamplingGrid::padWithZerosForEvenNumberedSamples() - new size "
+      +toString(getWindowGridPoints())+" is still not divisible by 8.");
+}
+
 void SamplingGrid::downsample(
     SamplingGrid& instance,