Changeset a6574a

Timestamp:

Aug 16, 2017, 9:26:38 AM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_with_BondGraph_continued

Parents:

08df4a

git-author:

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

git-committer:

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

Message:

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

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

File:

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

src/Dynamics/ForceAnnealing.hpp ¶

@@ -591 +591 @@
 
     // 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) {
@@ -600 +601 @@
       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;