Changeset 38ac25 for src/Atom

Timestamp:

Nov 7, 2017, 7:34:56 AM (7 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_with_BondGraph_continued_betteresults

Children:

10b1ef

Parents:

0f9726

git-author:

Frederik Heber <frederik.heber@…> (08/03/17 10:47:26)

git-committer:

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

Message:

FIX: ForceAnnealingAction's Undo is again working.

• AtomicInfo::setAtom() and SetAtomsFromAtomicInfo() now accept a time step. This allows to undo changes to the old state (where forces might have been passed to) and the current state (modified by time integration).
• TESTFIX: ForceAnnealing's undo test.conf was no longer valid since the additional time stepping is extracted (and the tests were changed already). I.e. the extra trajectory step is because of step-world-time and not because of force-annealing. Hence, it should not be undone here.

Location:

src/Atom

Files:

• AtomicInfo.cpp (modified) (2 diffs)
• AtomicInfo.hpp (modified) (2 diffs)

src/Atom/AtomicInfo.cpp

@@ -58 +58 @@
 {}
 
-AtomicInfo::AtomicInfo(const atom &_atom) :
-    Position(_atom.getPosition()),
+AtomicInfo::AtomicInfo(const atom &_atom, unsigned int _step) :
+    Position(_atom.getPositionAtStep(_step)),
     Type(_atom.getType()),
     charge(_atom.getCharge()),
-    Velocity(_atom.getAtomicVelocity()),
-    Force(_atom.getAtomicForce()),
+    Velocity(_atom.getAtomicVelocityAtStep(_step)),
+    Force(_atom.getAtomicForceAtStep(_step)),
     FatherId(_atom.father->getId()),
     MolId(0),
@@ -101 +101 @@
 
 
-bool AtomicInfo::setAtom(atom &_atom) const
+bool AtomicInfo::setAtom(atom &_atom, const unsigned int _step) const
 {
-  _atom.setPosition(Position);
+  _atom.setPositionAtStep(_step, Position);
   _atom.setType(Type);
   _atom.setCharge(charge);
-  _atom.setAtomicVelocity(Velocity);
-  _atom.setAtomicForce(Force);
+  _atom.setAtomicVelocityAtStep(_step, Velocity);
+  _atom.setAtomicForceAtStep(_step, Force);
 
   // set old id

src/Atom/AtomicInfo.hpp

@@ -20 +20 @@
 
 #include "Bond/BondInfo.hpp"
+#include "WorldTime.hpp"
 
 class atom;
@@ -31 +32 @@
 public:
   AtomicInfo();
-  AtomicInfo(const atom &_atom);
+  AtomicInfo(const atom &_atom, const unsigned int _step = WorldTime::getTime());
   AtomicInfo(const AtomicInfo &_atominfo);
   ~AtomicInfo();
 
-  bool setAtom(atom &_atom) const;
+  bool setAtom(atom &_atom, const unsigned int _step = WorldTime::getTime()) const;
   atomId_t getId() const;