Changeset 1e4f0a for src/Dynamics/ForceAnnealing.hpp

Timestamp:

Sep 13, 2017, 5:20:45 PM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Children:

1db3520

Parents:

15acdb

Message:

ForceAnnealing returns whether to stop and ForceAnnealingAction stops MakroAction if inside such.

• this is done by setting the MakroAction's loop variable to the current step which is now accessible.
• ForceAnnealing currently has a fixed value of 1e-8 when the force components are below to stop annealing.

File:

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

src/Dynamics/ForceAnnealing.hpp

@@ -78 +78 @@
     maxSteps(_maxSteps),
     max_distance(_max_distance),
-    damping_factor(_damping_factor)
+    damping_factor(_damping_factor),
+    FORCE_THRESHOLD(1e-8)
   {}
 
@@ -94 +95 @@
    * \param _TimeStep 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
+   * \return false - need to continue annealing, true - may stop because forces very small
    * \todo This is not yet checked if it is correctly working with DoConstrainedMD set >0.
    */
-  void operator()(
+  bool operator()(
       const int _TimeStep,
       const size_t _offset,
@@ -134 +136 @@
         << currentStep << " are " << maxComponents);
 
+    // check whether are smaller than threshold
+    bool AnnealingFinished = false;
+    double maxcomp = 0.;
+    for (size_t i=0;i<NDIM;++i)
+      maxcomp = std::max(maxcomp, fabs(maxComponents[i]));
+    if (maxcomp < FORCE_THRESHOLD) {
+      LOG(1, "STATUS: Force components are all less than " << FORCE_THRESHOLD
+          << ", stopping.");
+      currentStep = maxSteps;
+      AnnealingFinished = true;
+    }
+
     // are we in final step? Remember to reset static entities
     if (currentStep == maxSteps) {
@@ -139 +153 @@
       reset();
     }
+
+    return AnnealingFinished;
   }
 
@@ -627 +643 @@
   //!> the shifted is dampened by this factor with the power of the bond graph distance to the shift causing atom
   const double damping_factor;
+  //!> threshold for force components to stop annealing
+  const double FORCE_THRESHOLD;
 };