Changeset c9f4a0 for src/Dynamics/unittests/BondVectorsUnitTest.cpp

Timestamp:

Jul 20, 2017, 9:38:38 AM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_with_BondGraph_continued

Children:

6d7c73

Parents:

d3d964

git-author:

Frederik Heber <frederik.heber@…> (06/29/17 15:46:02)

git-committer:

Frederik Heber <frederik.heber@…> (07/20/17 09:38:38)

Message:

tempcommit: Added failing weights test case to BondVectors.

File:

• src/Dynamics/unittests/BondVectorsUnitTest.cpp (modified) (3 diffs)

src/Dynamics/unittests/BondVectorsUnitTest.cpp

@@ -67 +67 @@
   // failing asserts should be thrown
   ASSERT_DO(Assert::Throw);
+
+  setVerbosity(4);
 
   // create an atom
@@ -258 +260 @@
   // calculate weights
   BondVectors::weights_t weights = bv->getWeightsForAtomAtStep(*atoms[center], WorldTime::getTime());
-  LOG(2, "DEBUG: triangle weights are " << weights);
+  LOG(2, "DEBUG: Triangle weights are " << weights);
   // check number of weights
   CPPUNIT_ASSERT_EQUAL( weights.size(), (size_t)3 );
@@ -265 +267 @@
   CPPUNIT_ASSERT( fabs(weight_sum - 2.) < 1e-10 );
 }
+
+/** Test whether calculating weights works on complex config
+ *
+ */
+void BondVectorsTest::weights_complexTest()
+{
+  // gather atoms
+  atomvector += atoms[center], atoms[left], atoms[right], atoms[top], atoms[topright], atoms [bottomright];
+  // create bonds
+  bondvector +=
+      atoms[center]->addBond(atoms[left]),
+      atoms[center]->addBond(atoms[right]),
+      atoms[center]->addBond(atoms[top]),
+      atoms[center]->addBond(atoms[topright]),
+      atoms[center]->addBond(atoms[bottomright]);
+  // prepare bondvectors
+  bv->setFromAtomRange< std::vector<atom *> >(atomvector.begin(), atomvector.end(), WorldTime::getTime());
+  // calculate weights
+  BondVectors::weights_t weights = bv->getWeightsForAtomAtStep(*atoms[center], WorldTime::getTime());
+  LOG(2, "DEBUG: Complex weights are " << weights);
+  // check number of weights
+  CPPUNIT_ASSERT_EQUAL( weights.size(), (size_t)5 );
+  // check sum of weights: one linear independent, two dependent vectors = 1 + 2*0.5
+  const double weight_sum = std::accumulate(weights.begin(), weights.end(), 0.);
+  CPPUNIT_ASSERT( fabs(weight_sum - 2.) < 1e-10 );
+}