Changeset 8a91f5 for src/Fragmentation/Summation/OrthogonalSummation_impl.hpp

Timestamp:

Oct 20, 2016, 8:52:32 PM (9 years ago)

Author:

Frederik Heber <heber@…>

Branches:

Fix_IndependentFragmentGrids

Children:

8d6a44

Parents:

1e6cf0

git-author:

Frederik Heber <heber@…> (09/22/16 05:46:50)

git-committer:

Frederik Heber <heber@…> (10/20/16 20:52:32)

Message:

FIX: (Orthogonal)Summation::Sum() summed small fragments then large ones.

• added rbegin() and rend() to SetValueMap's interface.
• this fails for the new individual fragment SamplingGrids because there we assume that we always sum a grid of smaller extension onto a larger one. now, we go in reverse order through the list of setvalues and sum from largest down to smallest.

File:

• src/Fragmentation/Summation/OrthogonalSummation_impl.hpp (modified) (1 diff)

src/Fragmentation/Summation/OrthogonalSummation_impl.hpp

@@ -72 +72 @@
 {
   T sum(zeroinstance);
-  for(typename SetValueMap<T>::const_iterator iter = setvalues.begin();
-      (iter != setvalues.end()) && (iter->first->size() <= level); ++iter) {
+  // go from largest fragments to smallest
+  for(typename SetValueMap<T>::const_reverse_iterator iter = setvalues.rbegin();
+      iter != setvalues.rend(); ++iter) {
+    if (iter->first->size() > level)
+      continue;
     const T tempvalue = (iter->second)->getContribution();
     sum += tempvalue;