Changeset 1050ac for src/Dynamics

Timestamp:

Jul 20, 2017, 9:38:38 AM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_with_BondGraph_continued

Children:

8270be

Parents:

44441a

git-author:

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

git-committer:

Frederik Heber <frederik.heber@…> (07/20/17 09:38:38)

Message:

Sorting given selected atoms by ids in ForceAnnealing.

File:

• src/Dynamics/AtomicForceManipulator.hpp (modified) (1 diff)

src/Dynamics/AtomicForceManipulator.hpp

@@ -82 +82 @@
     Vector tempVector;
     size_t i=0;
+
     for(typename AtomSetMixin<T>::iterator iter = atoms.begin(); iter != atoms.end(); ++iter,++i) {
       for(size_t d=0;d<NDIM;d++) {
         tempVector[d] = Force.Matrix[0][i][d+offset]*(IsAngstroem ? AtomicLengthToAngstroem : 1.);
       }
+      LOG(3, "DEBUG: Adding force vector " << tempVector << " to atom " << **iter);
+      ASSERT( ((*iter)->getId()+1) == Force.Matrix[0][i][0],
+          "AtomicForceManipulator::addForceMatrixToAtomicForce() - mismatch in ids "
+          +toString(((*iter)->getId()+1))+" and file "+toString(Force.Matrix[0][i][0])+".");
       tempVector += (*iter)->getAtomicForceAtStep(TimeStep);
       (*iter)->setAtomicForceAtStep(TimeStep, tempVector);