Changeset 8ddd07 for src/Actions

Timestamp:

Apr 10, 2018, 6:43:11 AM (7 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

AutomationFragmentation_failures, Candidate_v1.6.1, ChemicalSpaceEvaluator, Enhanced_StructuralOptimization_continued, Exclude_Hydrogens_annealWithBondGraph, ForceAnnealing_with_BondGraph, ForceAnnealing_with_BondGraph_contraction-expansion, Gui_displays_atomic_force_velocity, PythonUI_with_named_parameters, StoppableMakroAction, TremoloParser_IncreasedPrecision

Children:

ecfcf6

Parents:

0477b0

git-author:

Frederik Heber <frederik.heber@…> (06/27/17 21:02:47)

git-committer:

Frederik Heber <frederik.heber@…> (04/10/18 06:43:11)

Message:

Sorting given selected atoms by ids in ForceAnnealing.

File:

• src/Actions/MoleculeAction/ForceAnnealingAction.cpp (modified) (2 diffs)

src/Actions/MoleculeAction/ForceAnnealingAction.cpp

@@ -71 +71 @@
     return Action::failure;
   }
+  // first, we need to sort the mixin according to their ids (as selected atoms are sorted
+  // according to their arbitrary address in memory)
+  set.sortByIds();
+
+  // create undo state for all selected atoms (undo info)
+  std::vector<AtomicInfo> UndoInfo;
+  UndoInfo.reserve(set.size());
+  {
+    for (World::AtomSelectionConstIterator iter = World::getInstance().beginAtomSelection();
+        iter != World::getInstance().endAtomSelection();
+        ++iter)
+      UndoInfo.push_back(AtomicInfo(*(iter->second)));
+  }
+
   // we always operate relative to current time step, except on single debug output
   size_t CurrentStep = WorldTime::getInstance().getTime();
@@ -103 +117 @@
     else
       LOG(2, "File " << params.forcesfile.get() << " found and parsed.");
-  }
-
-  // create undo state for all selected atoms (undo info)
-  std::vector<AtomicInfo> UndoInfo;
-  UndoInfo.reserve(set.size());
-  {
-    for (World::AtomSelectionConstIterator iter = World::getInstance().beginAtomSelection();
-        iter != World::getInstance().endAtomSelection();
-        ++iter)
-      UndoInfo.push_back(AtomicInfo(*(iter->second)));
   }