Changeset d34446 for src/Dynamics/ForceAnnealing.hpp

Timestamp:

Nov 12, 2017, 8:40:54 AM (7 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_with_BondGraph_continued_betteresults

Children:

e96f35

Parents:

59634f

git-author:

Frederik Heber <frederik.heber@…> (06/20/17 20:26:12)

git-committer:

Frederik Heber <frederik.heber@…> (11/12/17 08:40:54)

Message:

ForceAnnealing can now be used either atom- or bond-centered.

• atom is old annealing method, bond graph is new annealing method where the bond neighborhood is shifted as well to anneal to force projected onto the BondVector.
• In the first step we always use atom-centered annealing where we use the given deltat step width. Afterwards Barzilai-Borwein may then proceed.

File:

• src/Dynamics/ForceAnnealing.hpp (modified) (3 diffs)

src/Dynamics/ForceAnnealing.hpp

@@ -125 +125 @@
     }
 
+
     LOG(1, "STATUS: Largest remaining force components at step #"
         << currentStep << " are " << maxComponents);
@@ -219 +220 @@
    *
    *
-   * \param CurrentTimeStep current time step (i.e. t where \f$ t + \Delta t \f$ is in the sense of the velocity verlet)
+   * \param CurrentTimeStep current time step (i.e. \f$ t + \Delta t \f$ in the sense of the velocity verlet)
    * \param offset offset in matrix file to the first force component
    * \param maxComponents to be filled with maximum force component over all atoms
@@ -538 +539 @@
           walker->getPositionAtStep(CurrentTimeStep-1>=0 ? CurrentTimeStep - 1 : 0)
           + update);
+      walker->setAtomicVelocity(update);
 //      walker->setAtomicForce( RemnantGradient_per_atom[walker->getId()] );
     }