Changeset 14a64b for src/Actions/FragmentationAction

Timestamp:

Nov 7, 2017, 7:33:51 AM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_with_BondGraph_continued_betteresults

Children:

0e6205

Parents:

d9632e

git-author:

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

git-committer:

Frederik Heber <frederik.heber@…> (11/07/17 07:33:51)

Message:

ForceAnnealing no longer resets forces, AnalyseFragmentResults now sets forces.

• before we added onto force vector in AnalyseFragmentResults.
• if we do not reset forces in ForceAnnealing we can inspect the force components in the output file along the annealing run.
• TESTFIX: data files for regression test AnalyseFragmentResults OBC and PBC changed because of modified way of resetting forces.

File:

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

src/Actions/FragmentationAction/AnalyseFragmentationResultsAction.cpp

@@ -585 +585 @@
 }
 
-static void AddForces(
+static void SetForces(
     const IndexedVectors::indexedvectors_t &_forces,
     const bool _IsAngstroem)
@@ -603 +603 @@
     atom *_atom = World::getInstance().getAtom(AtomById(index));
     if(_atom != NULL)
-      _atom->setAtomicForce(_atom->getAtomicForce() + ForceVector);
+      _atom->setAtomicForce(ForceVector);
     else
       ELOG(2, "Could not find atom to id " << index << ".");
@@ -687 +687 @@
             shortrangeresults.Result_Force_fused.back()
             ).getVectors();
-    AddForces(shortrange_forces,IsAngstroem);
+    SetForces(shortrange_forces,IsAngstroem);
   } else {
     LOG(1, "INFO: Full molecule not loaded, hence will not add forces to atoms.");
@@ -760 +760 @@
               longrangeresults.Result_ForcesLongRangeIntegrated_fused.back()
               ).getVectors();
-      AddForces(longrange_forces,IsAngstroem);
+      SetForces(longrange_forces,IsAngstroem);
     } else {
       LOG(1, "INFO: Full molecule not loaded, hence will not add forces to atoms.");