Changes in src/moleculelist.cpp [35b698:84c494]

File:

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

src/moleculelist.cpp

@@ -24 +24 @@
 #include "periodentafel.hpp"
 #include "Helpers/Assert.hpp"
+#include "Matrix.hpp"
+#include "Box.hpp"
 
 #include "Helpers/Assert.hpp"
@@ -631 +633 @@
   int FragmentCounter = 0;
   ofstream output;
-  double cell_size_backup[6];
-  double * const cell_size = World::getInstance().getDomain();
-
-  // backup cell_size
-  for (int i=0;i<6;i++)
-    cell_size_backup[i] = cell_size[i];
+  Matrix cell_size = World::getInstance().getDomain().getM();
+  Matrix cell_size_backup = cell_size;
+
   // store the fragments as config and as xyz
   for (MoleculeList::iterator ListRunner = ListOfMolecules.begin(); ListRunner != ListOfMolecules.end(); ListRunner++) {
@@ -674 +673 @@
     (*ListRunner)->CenterEdge(&BoxDimension);
     (*ListRunner)->SetBoxDimension(&BoxDimension); // update Box of atoms by boundary
-    int j = -1;
     for (int k = 0; k < NDIM; k++) {
-      j += k + 1;
       BoxDimension[k] = 2.5 * (World::getInstance().getConfig()->GetIsAngstroem() ? 1. : 1. / AtomicLengthToAngstroem);
-      cell_size[j] = BoxDimension[k] * 2.;
-    }
+      cell_size.at(k,k) = BoxDimension[k] * 2.;
+    }
+    World::getInstance().setDomain(cell_size);
     (*ListRunner)->Translate(&BoxDimension);
 
@@ -725 +723 @@
 
   // restore cell_size
-  for (int i=0;i<6;i++)
-    cell_size[i] = cell_size_backup[i];
+  World::getInstance().setDomain(cell_size_backup);
 
   return result;
@@ -887 +884 @@
   // center at set box dimensions
   mol->CenterEdge(&center);
-  World::getInstance().getDomain()[0] = center[0];
-  World::getInstance().getDomain()[1] = 0;
-  World::getInstance().getDomain()[2] = center[1];
-  World::getInstance().getDomain()[3] = 0;
-  World::getInstance().getDomain()[4] = 0;
-  World::getInstance().getDomain()[5] = center[2];
+  Matrix domain;
+  for(int i =0;i<NDIM;++i)
+    domain.at(i,i) = center[i];
+  World::getInstance().setDomain(domain);
   insert(mol);
 }