Changeset f58e56 for src/Actions

Timestamp:

Oct 17, 2013, 7:24:06 AM (11 years ago)

Author:

Frederik Heber <heber@…>

Branches:

Action_Thermostats, Add_AtomRandomPerturbation, Add_FitFragmentPartialChargesAction, Add_RotateAroundBondAction, Add_SelectAtomByNameAction, Added_ParseSaveFragmentResults, AddingActions_SaveParseParticleParameters, Adding_Graph_to_ChangeBondActions, Adding_MD_integration_tests, Adding_ParticleName_to_Atom, Adding_StructOpt_integration_tests, AtomFragments, Automaking_mpqc_open, AutomationFragmentation_failures, Candidate_v1.5.4, Candidate_v1.6.0, Candidate_v1.6.1, ChangeBugEmailaddress, ChangingTestPorts, ChemicalSpaceEvaluator, CombiningParticlePotentialParsing, Combining_Subpackages, Debian_Package_split, Debian_package_split_molecuildergui_only, Disabling_MemDebug, Docu_Python_wait, EmpiricalPotential_contain_HomologyGraph, EmpiricalPotential_contain_HomologyGraph_documentation, Enable_parallel_make_install, Enhance_userguide, Enhanced_StructuralOptimization, Enhanced_StructuralOptimization_continued, Example_ManyWaysToTranslateAtom, Exclude_Hydrogens_annealWithBondGraph, FitPartialCharges_GlobalError, Fix_BoundInBox_CenterInBox_MoleculeActions, Fix_ChargeSampling_PBC, Fix_ChronosMutex, Fix_FitPartialCharges, Fix_FitPotential_needs_atomicnumbers, Fix_ForceAnnealing, Fix_IndependentFragmentGrids, Fix_ParseParticles, Fix_ParseParticles_split_forward_backward_Actions, Fix_PopActions, Fix_QtFragmentList_sorted_selection, Fix_Restrictedkeyset_FragmentMolecule, Fix_StatusMsg, Fix_StepWorldTime_single_argument, Fix_Verbose_Codepatterns, Fix_fitting_potentials, Fixes, ForceAnnealing_goodresults, ForceAnnealing_oldresults, ForceAnnealing_tocheck, ForceAnnealing_with_BondGraph, ForceAnnealing_with_BondGraph_continued, ForceAnnealing_with_BondGraph_continued_betteresults, ForceAnnealing_with_BondGraph_contraction-expansion, FragmentAction_writes_AtomFragments, FragmentMolecule_checks_bonddegrees, GeometryObjects, Gui_Fixes, Gui_displays_atomic_force_velocity, ImplicitCharges, IndependentFragmentGrids, IndependentFragmentGrids_IndividualZeroInstances, IndependentFragmentGrids_IntegrationTest, IndependentFragmentGrids_Sole_NN_Calculation, JobMarket_RobustOnKillsSegFaults, JobMarket_StableWorkerPool, JobMarket_unresolvable_hostname_fix, MoreRobust_FragmentAutomation, ODR_violation_mpqc_open, PartialCharges_OrthogonalSummation, PdbParser_setsAtomName, PythonUI_with_named_parameters, QtGui_reactivate_TimeChanged_changes, Recreated_GuiChecks, Rewrite_FitPartialCharges, RotateToPrincipalAxisSystem_UndoRedo, SaturateAtoms_findBestMatching, SaturateAtoms_singleDegree, StoppableMakroAction, Subpackage_CodePatterns, Subpackage_JobMarket, Subpackage_LinearAlgebra, Subpackage_levmar, Subpackage_mpqc_open, Subpackage_vmg, Switchable_LogView, ThirdParty_MPQC_rebuilt_buildsystem, TrajectoryDependenant_MaxOrder, TremoloParser_IncreasedPrecision, TremoloParser_MultipleTimesteps, TremoloParser_setsAtomName, Ubuntu_1604_changes, stable

Children:

47cee7

Parents:

0588e9

git-author:

Frederik Heber <heber@…> (09/29/13 09:45:35)

git-committer:

Frederik Heber <heber@…> (10/17/13 07:24:06)

Message:

Extracted cycle contribution calculation into distinct function.

File:

• src/Actions/FragmentationAction/AnalyseFragmentationResultsAction.cpp (modified) (2 diffs)

src/Actions/FragmentationAction/AnalyseFragmentationResultsAction.cpp

@@ -374 +374 @@
 }
 
+
+/** Determines the largest cycle in the container and returns its size.
+ *
+ * \param cycles set of cycles
+ * \return size if largest cycle
+ */
+size_t getMaxCycleSize(const KeySetsContainer &cycles)
+{
+  // gather cycle sizes
+  std::vector<size_t> cyclesizes(cycles.KeySets.size());
+  std::transform(
+      cycles.KeySets.begin(), cycles.KeySets.end(),
+      cyclesizes.begin(),
+      boost::bind(&KeySetsContainer::IntVector::size, boost::lambda::_1)
+  );
+  // get maximum
+  std::vector<size_t>::const_iterator maximum_size =
+      std::max_element(cyclesizes.begin(), cyclesizes.end());
+  if (maximum_size != cyclesizes.end())
+    return *maximum_size;
+  else
+    return 0;
+}
+
+void calculateCycleFullContribution(
+    const std::map<JobId_t, MPQCData> &shortrangedata,
+    const KeySetsContainer &keysets,
+    const KeySetsContainer &forcekeysets,
+    const KeySetsContainer &cycles,
+    const FragmentationShortRangeResults &shortrangeresults)
+{
+  // copy the shortrangeresults such that private MaxLevel is set in
+  // FragmentationShortRangeResults
+  FragmentationShortRangeResults cycleresults(shortrangeresults);
+  // get largest size
+  const size_t maximum_size = getMaxCycleSize(cycles);
+
+  /// The idea here is that (Orthogonal)Summation will place a result
+  /// consisting of level 1,2, and 3 fragment and a level 6 ring nonetheless
+  /// in level 6. If we want to have this result already at level 3, we
+  /// have to specifically inhibit all fragments from later levels but the
+  /// cycles and then pick the result from the last level and placing it at
+  /// the desired one
+
+  // loop from level 1 to max ring size and gather contributions
+  for (size_t level = 1; level <= maximum_size; ++level) {
+    // create ValueMask for this level by stepping through each keyset and checking size
+    std::vector<bool> localValueMask(shortrangedata.size(), false);
+    size_t index=0;
+    // TODO: if only KeySetsContainer was usable as a compliant STL container, might be able to use set_difference or alike.
+    KeySetsContainer::ArrayOfIntVectors::const_iterator keysetsiter = keysets.KeySets.begin();
+    KeySetsContainer::ArrayOfIntVectors::const_iterator cyclesiter = cycles.KeySets.begin();
+    for (; (keysetsiter != keysets.KeySets.end()) && (cyclesiter != cycles.KeySets.end());) {
+      if (cyclesiter->size() > keysetsiter->size()) {
+        // add if not greater than level in size
+        if ((*keysetsiter).size() <= level)
+          localValueMask[index] = true;
+        ++keysetsiter;
+        ++index;
+      } else if (cyclesiter->size() < keysetsiter->size()) {
+        ++cyclesiter;
+      } else { // both sets have same size
+        if (*cyclesiter > *keysetsiter) {
+          // add if not greater than level in size
+          if ((*keysetsiter).size() <= level)
+            localValueMask[index] = true;
+          ++keysetsiter;
+          ++index;
+        } else if (*cyclesiter < *keysetsiter) {
+          ++cyclesiter;
+        } else {
+          // also always add all cycles
+          localValueMask[index] = true;
+          ++cyclesiter;
+          ++keysetsiter;
+          ++index;
+        }
+      }
+    }
+    // activate rest if desired by level
+    for (; keysetsiter != keysets.KeySets.end(); ++keysetsiter) {
+      if ((*keysetsiter).size() <= level)
+        localValueMask[index] = true;
+      ++index;
+    }
+    LOG(2, "DEBUG: ValueMask for cycle correction at level " << level << " is "
+        << localValueMask << ".");
+    // create FragmentationShortRangeResults
+    FragmentationShortRangeResults localresults(shortrangedata, keysets, forcekeysets, localValueMask);
+    // and perform summation
+    localresults(shortrangedata);
+    // finally, extract the corrections from last level
+    cycleresults.Result_Energy_fused[level-1] =
+        localresults.Result_Energy_fused.back();
+    cycleresults.Result_Time_fused[level-1] =
+        localresults.Result_Time_fused.back();
+    cycleresults.Result_Force_fused[level-1] =
+        localresults.Result_Force_fused.back();
+  }
+  printReceivedCycleResults(cycleresults);
+}
 
 ActionState::ptr FragmentationAnalyseFragmentationResultsAction::performCall() {
@@ -444 +545 @@
   }
 
-  // copy the shortrangeresults such that private MaxLevel is set in
-  // FragmentationShortRangeResults
-  FragmentationShortRangeResults cycleresults(shortrangeresults);
   // calculate energy if cycles are calculated fully at each level already
   if (!cycles.KeySets.empty()) {
-    // gather cycle sizes
-    std::vector<size_t> cyclesizes(cycles.KeySets.size());
-    std::transform(
-        cycles.KeySets.begin(), cycles.KeySets.end(),
-        cyclesizes.begin(),
-        boost::bind(&KeySetsContainer::IntVector::size, boost::lambda::_1)
-    );
-    // get maximum
-    std::vector<size_t>::const_iterator maximum_size =
-        std::max_element(cyclesizes.begin(), cyclesizes.end());
-
-    // loop from level 1 to max ring size and gather corrections
-    for (size_t level = 1; level <= *maximum_size; ++level) {
-      // create ValueMask for this level by stepping through each keyset and checking size
-      std::vector<bool> localValueMask(shortrangedata.size(), false);
-      size_t index=0;
-      KeySetsContainer::ArrayOfIntVectors::const_iterator keysetsiter = keysets.KeySets.begin();
-      KeySetsContainer::ArrayOfIntVectors::const_iterator cyclesiter = cycles.KeySets.begin();
-      for (; (keysetsiter != keysets.KeySets.end()) && (cyclesiter != cycles.KeySets.end());) {
-        if (cyclesiter->size() > keysetsiter->size()) {
-          // add if not greater than level in size
-          if ((*keysetsiter).size() <= level)
-            localValueMask[index] = true;
-          ++keysetsiter;
-          ++index;
-        } else if (cyclesiter->size() < keysetsiter->size()) {
-          ++cyclesiter;
-        } else { // both sets have same size
-          if (*cyclesiter > *keysetsiter) {
-            // add if not greater than level in size
-            if ((*keysetsiter).size() <= level)
-              localValueMask[index] = true;
-            ++keysetsiter;
-            ++index;
-          } else if (*cyclesiter < *keysetsiter) {
-            ++cyclesiter;
-          } else {
-            // also always add all cycles
-            localValueMask[index] = true;
-            ++cyclesiter;
-            ++keysetsiter;
-            ++index;
-          }
-        }
-      }
-      // activate rest if desired by level
-      for (; keysetsiter != keysets.KeySets.end(); ++keysetsiter) {
-        if ((*keysetsiter).size() <= level)
-          localValueMask[index] = true;
-        ++index;
-      }
-      LOG(2, "DEBUG: ValueMask for cycle correction at level " << level << " is "
-          << localValueMask << ".");
-      // create FragmentationShortRangeResults
-      FragmentationShortRangeResults localresults(shortrangedata, keysets, forcekeysets, localValueMask);
-      // and perform summation
-      localresults(shortrangedata);
-      // finally, extract the corrections
-      cycleresults.Result_Energy_fused[level-1] =
-          localresults.Result_Energy_fused.back();
-      cycleresults.Result_Time_fused[level-1] =
-          localresults.Result_Time_fused.back();
-      cycleresults.Result_Force_fused[level-1] =
-          localresults.Result_Force_fused.back();
-    }
-    printReceivedCycleResults(cycleresults);
+    calculateCycleFullContribution(
+        shortrangedata,
+        keysets,
+        forcekeysets,
+        cycles,
+        shortrangeresults);
   }