Changeset ccbdf4 for src/Dynamics/BondVectors.cpp

Timestamp:

Aug 2, 2017, 8:25:57 PM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_with_BondGraph_continued

Children:

c4a725

Parents:

90e540

Message:

tempcommit: BondVectors::getRemnant...() now requires atom's gradient.

File:

• src/Dynamics/BondVectors.cpp (modified) (4 diffs)

src/Dynamics/BondVectors.cpp

@@ -287 +287 @@
 Vector BondVectors::getRemnantGradientForAtomAtStep(
     const atom &_walker,
+    const Vector &_walkerGradient,
     const std::vector<Vector> _BondVectors,
     const BondVectors::weights_t &_weights,
@@ -292 +293 @@
     forcestore_t _forcestore) const
 {
-  const Vector &walkerGradient = _walker.getAtomicForceAtStep(_step);
   BondVectors::weights_t::const_iterator weightiter = _weights.begin();
   std::vector<Vector>::const_iterator vectoriter = _BondVectors.begin();
@@ -303 +303 @@
     const Vector &BondVector = *vectoriter;
 
-    const double temp = (*weightiter)*walkerGradient.ScalarProduct(BondVector);
+    const double temp = (*weightiter)*_walkerGradient.ScalarProduct(BondVector);
     _forcestore(_walker, current_bond, _step, temp);
     LOG(4, "DEBUG: BondVector " << BondVector << " receives projected force of "
-        << temp);
+        << (*weightiter) << "*" << _walkerGradient << "*" << BondVector << " = " << temp);
     forcesum += temp * BondVector;
   }
@@ -314 +314 @@
       "BondVectors::getRemnantGradientForAtomAtStep() - vectoriter is not at end when it should be.");
 
-  return walkerGradient-forcesum;
+  return _walkerGradient-forcesum;
 }