Changeset 66e1c6

Timestamp:

Sep 24, 2017, 12:24:55 PM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Parents:

ab3827

git-author:

Frederik Heber <frederik.heber@…> (08/16/17 09:25:45)

git-committer:

Frederik Heber <frederik.heber@…> (09/24/17 12:24:55)

Message:

annealWithBondGraph_..() prints the largest update on default verbosity.

• this helps in debugging when the fragmentation during StructOpt fails.

File:

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

src/Dynamics/ForceAnnealing.hpp

@@ -590 +590 @@
 
     // apply the gathered updates and set remnant gradients for atomic annealing
+    Vector LargestUpdate;
     for (std::map<atomId_t, Vector>::const_iterator iter = GatheredUpdates.begin();
         iter != GatheredUpdates.end(); ++iter) {
@@ -599 +600 @@
       LOG(3, "DEBUG: Applying update " << update << " to atom #" << atomid
           << ", namely " << *walker);
+      for (size_t i=0;i<NDIM;++i)
+        LargestUpdate[i] = std::max(LargestUpdate[i], fabs(update[i]));
       walker->setPositionAtStep(_TimeStep,
           walker->getPositionAtStep(CurrentTimeStep) + update); // - CommonTranslation);
     }
+    LOG(1, "STATUS: Largest absolute update components are " << LargestUpdate);
 
     return maxComponents;