Changes in src/molecule_geometry.cpp [52d777:f429d7]

File:

• src/molecule_geometry.cpp (modified) (12 diffs)

src/molecule_geometry.cpp

@@ -19 +19 @@
 #include "World.hpp"
 #include "Plane.hpp"
+#include "Matrix.hpp"
+#include "Box.hpp"
 #include <boost/foreach.hpp>
 
@@ -33 +35 @@
   const Vector *Center = DetermineCenterOfAll();
   const Vector *CenterBox = DetermineCenterOfBox();
-  double * const cell_size = World::getInstance().getDomain();
-  double *M = ReturnFullMatrixforSymmetric(cell_size);
-  double *Minv = InverseMatrix(M);
+  Box &domain = World::getInstance().getDomain();
 
   // go through all atoms
   ActOnAllVectors( &Vector::SubtractVector, *Center);
   ActOnAllVectors( &Vector::SubtractVector, *CenterBox);
-  ActOnAllVectors( &Vector::WrapPeriodically, (const double *)M, (const double *)Minv);
-
-  delete[](M);
-  delete[](Minv);
+  BOOST_FOREACH(atom* iter, atoms){
+    *iter->node = domain.WrapPeriodically(*iter->node);
+  }
+
   delete(Center);
-  delete(CenterBox);
   return status;
 };
@@ -56 +55 @@
 {
   bool status = true;
-  double * const cell_size = World::getInstance().getDomain();
-  double *M = ReturnFullMatrixforSymmetric(cell_size);
-  double *Minv = InverseMatrix(M);
+  Box &domain = World::getInstance().getDomain();
 
   // go through all atoms
-  ActOnAllVectors( &Vector::WrapPeriodically, (const double *)M, (const double *)Minv);
-
-  delete[](M);
-  delete[](Minv);
+  BOOST_FOREACH(atom* iter, atoms){
+    *iter->node = domain.WrapPeriodically(*iter->node);
+  }
+
   return status;
 };
@@ -153 +150 @@
 {
   Vector *a = new Vector(0.5,0.5,0.5);
-
-  const double *cell_size = World::getInstance().getDomain();
-  double *M = ReturnFullMatrixforSymmetric(cell_size);
-  a->MatrixMultiplication(M);
-  delete[](M);
-
+  const Matrix &M = World::getInstance().getDomain().getM();
+  (*a) *= M;
   return a;
 };
@@ -244 +237 @@
 void molecule::TranslatePeriodically(const Vector *trans)
 {
-  double * const cell_size = World::getInstance().getDomain();
-  double *M = ReturnFullMatrixforSymmetric(cell_size);
-  double *Minv = InverseMatrix(M);
+  Box &domain = World::getInstance().getDomain();
 
   // go through all atoms
   ActOnAllVectors( &Vector::AddVector, *trans);
-  ActOnAllVectors( &Vector::WrapPeriodically, (const double *)M, (const double *)Minv);
-
-  delete[](M);
-  delete[](Minv);
+  BOOST_FOREACH(atom* iter, atoms){
+    *iter->node = domain.WrapPeriodically(*iter->node);
+  }
+
 };
 
@@ -264 +255 @@
   OBSERVE;
   Plane p(*n,0);
-  BOOST_FOREACH( atom* iter, atoms ){
+  BOOST_FOREACH(atom* iter, atoms ){
     (*iter->node) = p.mirrorVector(*iter->node);
   }
@@ -274 +265 @@
 void molecule::DeterminePeriodicCenter(Vector &center)
 {
-  double * const cell_size = World::getInstance().getDomain();
-  double *matrix = ReturnFullMatrixforSymmetric(cell_size);
-  double *inversematrix = InverseMatrix(matrix);
+  const Matrix &matrix = World::getInstance().getDomain().getM();
+  const Matrix &inversematrix = World::getInstance().getDomain().getM();
   double tmp;
   bool flag;
@@ -288 +278 @@
       if ((*iter)->type->Z != 1) {
 #endif
-        Testvector = (*iter)->x;
-        Testvector.MatrixMultiplication(inversematrix);
+        Testvector = inversematrix * (*iter)->x;
         Translationvector.Zero();
         for (BondList::const_iterator Runner = (*iter)->ListOfBonds.begin(); Runner != (*iter)->ListOfBonds.end(); (++Runner)) {
@@ -306 +295 @@
         }
         Testvector += Translationvector;
-        Testvector.MatrixMultiplication(matrix);
+        Testvector *= matrix;
         Center += Testvector;
         Log() << Verbose(1) << "vector is: " << Testvector << endl;
@@ -313 +302 @@
         for (BondList::const_iterator Runner = (*iter)->ListOfBonds.begin(); Runner != (*iter)->ListOfBonds.end(); (++Runner)) {
           if ((*Runner)->GetOtherAtom((*iter))->type->Z == 1) {
-            Testvector = (*Runner)->GetOtherAtom((*iter))->x;
-            Testvector.MatrixMultiplication(inversematrix);
+            Testvector = inversematrix * (*Runner)->GetOtherAtom((*iter))->x;
             Testvector += Translationvector;
-            Testvector.MatrixMultiplication(matrix);
+            Testvector *= matrix;
             Center += Testvector;
             Log() << Verbose(1) << "Hydrogen vector is: " << Testvector << endl;
@@ -325 +313 @@
     }
   } while (!flag);
-  delete[](matrix);
-  delete[](inversematrix);
 
   Center.Scale(1./static_cast<double>(getAtomCount()));
@@ -388 +374 @@
     DoLog(1) && (Log() << Verbose(1) << "Transforming molecule into PAS ... ");
     // the eigenvectors specify the transformation matrix
-    ActOnAllVectors( &Vector::MatrixMultiplication, (const double *) evec->data );
+    Matrix M = Matrix(evec->data);
+    BOOST_FOREACH(atom* iter, atoms){
+      (*iter->node) *= M;
+    }
     DoLog(0) && (Log() << Verbose(0) << "done." << endl);