Changeset a1eeb3

Timestamp:

Jun 9, 2016, 6:32:25 PM (9 years ago)

Author:

Frederik Heber <heber@…>

Parents:

e2fd2c

git-author:

Frederik Heber <heber@…> (06/07/16 08:00:51)

git-committer:

Frederik Heber <heber@…> (06/09/16 18:32:25)

Message:

FragmentationLongRangeResults::operator() does not need to calculate potential difference on full grid.

File:

• src/Fragmentation/Summation/Containers/FragmentationLongRangeResults.cpp (modified) (2 diffs)

src/Fragmentation/Summation/Containers/FragmentationLongRangeResults.cpp

@@ -172 +172 @@
       // first we obtain the full e-e energy as potential times charge on the
       // respective level.
-      const SamplingGrid &charge_weight =
-          boost::fusion::at_key<MPQCDataFused::sampled_grid>(Result_Grid_fused[level-1]);
+      const SamplingGrid &charge_weight = full_sample[level-1];
       SamplingGrid full_sample_solution = fullsolutionData[level-1].sampled_potential;
       full_sample_solution *= charge_weight;
@@ -183 +182 @@
           boost::fusion::at_key<VMGDataFused::sampled_potential>(Result_GridLongRange_fused[level-1]);
       double electron_short_range_energy = short_range_correction.integral();
-      full_sample_solution -= short_range_correction;
       electron_solution_energy -= electron_short_range_energy;
-      ASSERT( fabs(electron_solution_energy - full_sample_solution.integral()) < 1e-7,
-          "FragmentationLongRangeResults::operator() - integral and energy are not exchangeable.");
+#ifndef NDEBUG
+      {
+        static const double round_offset(
+            (std::numeric_limits<size_t>::round_style == std::round_toward_zero) ?
+                0.5 : 0.); // need offset to get to round_toward_nearest behavior
+        // we can only check equivalence if both have same level, otherwise
+        // short_range_correction.integral() has higher precision because of finger grid
+        const int surplus_level = full_sample_solution.getSurplusLevel(short_range_correction)+round_offset;
+        if (full_sample_solution.level == short_range_correction.level+surplus_level) {
+          SamplingGrid check_difference_full_solution = full_sample_solution;
+          check_difference_full_solution -= short_range_correction;
+          ASSERT( fabs(electron_solution_energy - check_difference_full_solution.integral()) < 1e-7,
+              "FragmentationLongRangeResults::operator() - integral and energy are not exchangeable.");
+        }
+      }
+#endif
 
       // then, we obtain the e-n+n-n full solution in the same way