Changes in src/molecule_geometry.cpp [c94eeb:112b09]

File:

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

src/molecule_geometry.cpp

@@ -19 +19 @@
 #include "World.hpp"
 #include "Plane.hpp"
-#include "Matrix.hpp"
 #include <boost/foreach.hpp>
 
@@ -155 +154 @@
 
   const double *cell_size = World::getInstance().getDomain();
-  double *M_double = ReturnFullMatrixforSymmetric(cell_size);
-  Matrix M = Matrix(M_double);
-  delete[](M_double);
+  double *M = ReturnFullMatrixforSymmetric(cell_size);
   a->MatrixMultiplication(M);
+  delete[](M);
 
   return a;
@@ -276 +274 @@
 {
   double * const cell_size = World::getInstance().getDomain();
-  double *matrix_double = ReturnFullMatrixforSymmetric(cell_size);
-  Matrix matrix = Matrix(matrix_double);
-  delete[](matrix_double);
-  Matrix inversematrix = matrix.invert();
+  double *matrix = ReturnFullMatrixforSymmetric(cell_size);
+  double *inversematrix = InverseMatrix(matrix);
   double tmp;
   bool flag;
@@ -328 +324 @@
     }
   } while (!flag);
+  delete[](matrix);
+  delete[](inversematrix);
 
   Center.Scale(1./static_cast<double>(getAtomCount()));
@@ -389 +387 @@
     DoLog(1) && (Log() << Verbose(1) << "Transforming molecule into PAS ... ");
     // the eigenvectors specify the transformation matrix
-    Matrix M = Matrix(evec->data);
-    ActOnAllVectors( &Vector::MatrixMultiplication, static_cast<const Matrix>(M));
+    ActOnAllVectors( &Vector::MatrixMultiplication, (const double *) evec->data );
     DoLog(0) && (Log() << Verbose(0) << "done." << endl);