Changeset 2a0a651 for src/Dynamics/ForceAnnealing.hpp

Timestamp:

Nov 29, 2017, 11:06:25 PM (7 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_with_BondGraph_continued_betteresults

Parents:

0dd99b

git-author:

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

git-committer:

Frederik Heber <frederik.heber@…> (11/29/17 23:06:25)

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

@@ -608 +608 @@
 
     // 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) {
@@ -617 +618 @@
       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;