Changes in src/boundary.cpp [986ed3:952f38]

File:

• src/boundary.cpp (modified) (9 diffs)

src/boundary.cpp

@@ -12 +12 @@
 #include "config.hpp"
 #include "element.hpp"
-#include "helpers.hpp"
-#include "info.hpp"
+#include "Helpers/helpers.hpp"
+#include "Helpers/Info.hpp"
 #include "linkedcell.hpp"
-#include "verbose.hpp"
-#include "log.hpp"
+#include "Helpers/Verbose.hpp"
+#include "Helpers/Log.hpp"
 #include "molecule.hpp"
 #include "tesselation.hpp"
 #include "tesselationhelpers.hpp"
 #include "World.hpp"
-#include "Plane.hpp"
-#include "Matrix.hpp"
+#include "LinearAlgebra/Plane.hpp"
+#include "LinearAlgebra/Matrix.hpp"
 #include "Box.hpp"
 
@@ -301 +301 @@
  * \param *out output stream for debugging
  * \param *mol molecule structure with Atom's and Bond's.
+ * \param *BoundaryPts set of boundary points to use or NULL
  * \param *TesselStruct Tesselation filled with points, lines and triangles on boundary on return
  * \param *LCList atoms in LinkedCell list
@@ -306 +307 @@
  * \return *TesselStruct is filled with convex boundary and tesselation is stored under \a *filename.
  */
-void FindConvexBorder(const molecule* mol, Tesselation *&TesselStruct, const LinkedCell *LCList, const char *filename)
+void FindConvexBorder(const molecule* mol, Boundaries *BoundaryPts, Tesselation *&TesselStruct, const LinkedCell *LCList, const char *filename)
 {
         Info FunctionInfo(__func__);
@@ -317 +318 @@
 
   // 1. Find all points on the boundary
-  if (BoundaryPoints == NULL) {
-      BoundaryFreeFlag = true;
-      BoundaryPoints = GetBoundaryPoints(mol, TesselStruct);
+  if (BoundaryPts == NULL) {
+    BoundaryFreeFlag = true;
+    BoundaryPoints = GetBoundaryPoints(mol, TesselStruct);
   } else {
-      DoLog(0) && (Log() << Verbose(0) << "Using given boundary points set." << endl);
+    BoundaryPoints = BoundaryPts;
+    DoLog(0) && (Log() << Verbose(0) << "Using given boundary points set." << endl);
   }
 
 // printing all inserted for debugging
-  for (int axis=0; axis < NDIM; axis++)
-    {
-      DoLog(1) && (Log() << Verbose(1) << "Printing list of candidates for axis " << axis << " which we have inserted so far." << endl);
-      int i=0;
-      for(Boundaries::iterator runner = BoundaryPoints[axis].begin(); runner != BoundaryPoints[axis].end(); runner++) {
-        if (runner != BoundaryPoints[axis].begin())
-          DoLog(0) && (Log() << Verbose(0) << ", " << i << ": " << *runner->second.second);
-        else
-          DoLog(0) && (Log() << Verbose(0) << i << ": " << *runner->second.second);
-        i++;
-      }
-      DoLog(0) && (Log() << Verbose(0) << endl);
-    }
+  for (int axis=0; axis < NDIM; axis++) {
+    DoLog(1) && (Log() << Verbose(1) << "Printing list of candidates for axis " << axis << " which we have inserted so far." << endl);
+    int i=0;
+    for(Boundaries::iterator runner = BoundaryPoints[axis].begin(); runner != BoundaryPoints[axis].end(); runner++) {
+      if (runner != BoundaryPoints[axis].begin())
+        DoLog(0) && (Log() << Verbose(0) << ", " << i << ": " << *runner->second.second);
+      else
+        DoLog(0) && (Log() << Verbose(0) << i << ": " << *runner->second.second);
+      i++;
+    }
+    DoLog(0) && (Log() << Verbose(0) << endl);
+  }
 
   // 2. fill the boundary point list
@@ -343 +344 @@
     for (Boundaries::iterator runner = BoundaryPoints[axis].begin(); runner != BoundaryPoints[axis].end(); runner++)
         if (!TesselStruct->AddBoundaryPoint(runner->second.second, 0))
-          DoeLog(2) && (eLog()<< Verbose(2) << "Point " << *(runner->second.second) << " is already present!" << endl);
+          DoLog(2) && (Log()<< Verbose(2) << "Point " << *(runner->second.second) << " is already present." << endl);
 
   DoLog(0) && (Log() << Verbose(0) << "I found " << TesselStruct->PointsOnBoundaryCount << " points on the convex boundary." << endl);
@@ -677 +678 @@
 
   IsAngstroem = configuration->GetIsAngstroem();
+  BoundaryPoints = GetBoundaryPoints(mol, TesselStruct);
   GreatestDiameter = GetDiametersOfCluster(BoundaryPoints, mol, TesselStruct, IsAngstroem);
-  BoundaryPoints = GetBoundaryPoints(mol, TesselStruct);
   LinkedCell *LCList = new LinkedCell(mol, 10.);
-  FindConvexBorder(mol, TesselStruct, (const LinkedCell *&)LCList, NULL);
+  FindConvexBorder(mol, BoundaryPoints, TesselStruct, (const LinkedCell *&)LCList, NULL);
   delete (LCList);
+  delete[] BoundaryPoints;
 
 
@@ -689 +691 @@
   else
     clustervolume = ClusterVolume;
+
+  delete TesselStruct;
 
   for (int i = 0; i < NDIM; i++)
@@ -741 +745 @@
     mol->CenterInBox();
   }
+  delete GreatestDiameter;
   // update Box of atoms by boundary
   mol->SetBoxDimension(&BoxLengths);
@@ -895 +900 @@
             }
       }
+  for (MoleculeList::iterator ListRunner = List->ListOfMolecules.begin(); ListRunner != List->ListOfMolecules.end(); ListRunner++) {
+    delete LCList[*ListRunner];
+    delete TesselStruct[(*ListRunner)];
+  }
 
   return Filling;