Changes in src/tesselationhelpers.cpp [d74077:952f38]

File:

• src/tesselationhelpers.cpp (modified) (18 diffs)

src/tesselationhelpers.cpp

@@ -10 +10 @@
 #include <fstream>
 
-#include "BoundaryLineSet.hpp"
-#include "BoundaryPointSet.hpp"
-#include "BoundaryPolygonSet.hpp"
-#include "BoundaryTriangleSet.hpp"
-#include "CandidateForTesselation.hpp"
-#include "info.hpp"
+#include "Helpers/Info.hpp"
 #include "linkedcell.hpp"
-#include "linearsystemofequations.hpp"
-#include "log.hpp"
+#include "LinearAlgebra/linearsystemofequations.hpp"
+#include "Helpers/Log.hpp"
 #include "tesselation.hpp"
 #include "tesselationhelpers.hpp"
-#include "vector.hpp"
-#include "Line.hpp"
+#include "LinearAlgebra/Vector.hpp"
+#include "LinearAlgebra/Line.hpp"
 #include "vector_ops.hpp"
-#include "verbose.hpp"
-#include "Plane.hpp"
-#include "Matrix.hpp"
+#include "Helpers/Verbose.hpp"
+#include "LinearAlgebra/Plane.hpp"
+#include "LinearAlgebra/Matrix.hpp"
 
 void GetSphere(Vector * const center, const Vector &a, const Vector &b, const Vector &c, const double RADIUS)
@@ -149 +144 @@
  * \param *c third point
  */
-void GetCenterofCircumcircle(Vector &Center, const Vector &a, const Vector &b, const Vector &c)
+void GetCenterofCircumcircle(Vector * const Center, const Vector &a, const Vector &b, const Vector &c)
 {
         Info FunctionInfo(__func__);
@@ -161 +156 @@
   helper[2] = SideC.NormSquared()*(SideA.NormSquared()+SideB.NormSquared() - SideC.NormSquared());
 
-  Center.Zero();
-  Center += helper[0] * a;
-  Center += helper[1] * b;
-  Center += helper[2] * c;
-  Center.Scale(1./(helper[0]+helper[1]+helper[2]));
-  Log() << Verbose(1) << "INFO: Center (2nd algo) is at " << Center << "." << endl;
+  Center->Zero();
+  *Center += helper[0] * a;
+  *Center += helper[1] * b;
+  *Center += helper[2] * c;
+  Center->Scale(1./(helper[0]+helper[1]+helper[2]));
+  Log() << Verbose(1) << "INFO: Center (2nd algo) is at " << *Center << "." << endl;
 };
 
@@ -393 +388 @@
  * @return point which is second closest to the provided one
  */
-TesselPoint* FindSecondClosestTesselPoint(const Vector& Point, const LinkedCell* const LC)
+TesselPoint* FindSecondClosestTesselPoint(const Vector* Point, const LinkedCell* const LC)
 {
         Info FunctionInfo(__func__);
@@ -417 +412 @@
         if (List != NULL) {
           for (LinkedCell::LinkedNodes::const_iterator Runner = List->begin(); Runner != List->end(); Runner++) {
-            helper = (Point) - ((*Runner)->getPosition());
+            helper = (*Point) - (*(*Runner)->node);
             double currentNorm = helper. Norm();
             if (currentNorm < distance) {
@@ -446 +441 @@
  * @return point which is closest to the provided one, NULL if none found
  */
-TesselPoint* FindClosestTesselPoint(const Vector& Point, TesselPoint *&SecondPoint, const LinkedCell* const LC)
+TesselPoint* FindClosestTesselPoint(const Vector* Point, TesselPoint *&SecondPoint, const LinkedCell* const LC)
 {
         Info FunctionInfo(__func__);
@@ -470 +465 @@
         if (List != NULL) {
           for (LinkedCell::LinkedNodes::const_iterator Runner = List->begin(); Runner != List->end(); Runner++) {
-            helper = (Point) - ((*Runner)->getPosition());
+            helper = (*Point) - (*(*Runner)->node);
             double currentNorm = helper.NormSquared();
             if (currentNorm < distance) {
@@ -508 +503 @@
         Info FunctionInfo(__func__);
   // construct the plane of the two baselines (i.e. take both their directional vectors)
-  Vector Baseline = (Base->endpoints[1]->node->getPosition()) - (Base->endpoints[0]->node->getPosition());
-  Vector OtherBaseline = (OtherBase->endpoints[1]->node->getPosition()) - (OtherBase->endpoints[0]->node->getPosition());
+  Vector Baseline = (*Base->endpoints[1]->node->node) - (*Base->endpoints[0]->node->node);
+  Vector OtherBaseline = (*OtherBase->endpoints[1]->node->node) - (*OtherBase->endpoints[0]->node->node);
   Vector Normal = Baseline;
   Normal.VectorProduct(OtherBaseline);
@@ -516 +511 @@
 
   // project one offset point of OtherBase onto this plane (and add plane offset vector)
-  Vector NewOffset = (OtherBase->endpoints[0]->node->getPosition()) - (Base->endpoints[0]->node->getPosition());
+  Vector NewOffset = (*OtherBase->endpoints[0]->node->node) - (*Base->endpoints[0]->node->node);
   NewOffset.ProjectOntoPlane(Normal);
-  NewOffset += (Base->endpoints[0]->node->getPosition());
+  NewOffset += (*Base->endpoints[0]->node->node);
   Vector NewDirection = NewOffset + OtherBaseline;
 
   // calculate the intersection between this projected baseline and Base
   Vector *Intersection = new Vector;
-  Line line1 = makeLineThrough((Base->endpoints[0]->node->getPosition()),(Base->endpoints[1]->node->getPosition()));
+  Line line1 = makeLineThrough(*(Base->endpoints[0]->node->node),*(Base->endpoints[1]->node->node));
   Line line2 = makeLineThrough(NewOffset, NewDirection);
   *Intersection = line1.getIntersection(line2);
-  Normal = (*Intersection) - (Base->endpoints[0]->node->getPosition());
+  Normal = (*Intersection) - (*Base->endpoints[0]->node->node);
   DoLog(1) && (Log() << Verbose(1) << "Found closest point on " << *Base << " at " << *Intersection << ", factor in line is " << fabs(Normal.ScalarProduct(Baseline)/Baseline.NormSquared()) << "." << endl);
 
@@ -571 +566 @@
       *vrmlfile << "Sphere {" << endl << "  "; // 2 is sphere type
       for (i=0;i<NDIM;i++)
-        *vrmlfile << Walker->at(i)-center->at(i) << " ";
+        *vrmlfile << Walker->node->at(i)-center->at(i) << " ";
       *vrmlfile << "\t0.1\t1. 1. 1." << endl; // radius 0.05 and white as colour
       cloud->GoToNext();
@@ -581 +576 @@
       for (i=0;i<3;i++) { // print each node
         for (int j=0;j<NDIM;j++)  // and for each node all NDIM coordinates
-          *vrmlfile << TriangleRunner->second->endpoints[i]->node->at(j)-center->at(j) << " ";
+          *vrmlfile << TriangleRunner->second->endpoints[i]->node->node->at(j)-center->at(j) << " ";
         *vrmlfile << "\t";
       }
@@ -608 +603 @@
     Vector *center = cloud->GetCenter();
     // make the circumsphere's center absolute again
-    Vector helper = (1./3.) * ((Tess->LastTriangle->endpoints[0]->node->getPosition()) +
-                               (Tess->LastTriangle->endpoints[1]->node->getPosition()) +
-                               (Tess->LastTriangle->endpoints[2]->node->getPosition()));
+    Vector helper = (1./3.) * ((*Tess->LastTriangle->endpoints[0]->node->node) +
+                               (*Tess->LastTriangle->endpoints[1]->node->node) +
+                               (*Tess->LastTriangle->endpoints[2]->node->node));
     helper -= (*center);
     // and add to file plus translucency object
@@ -643 +638 @@
       *rasterfile << "2" << endl << "  ";  // 2 is sphere type
       for (int j=0;j<NDIM;j++) { // and for each node all NDIM coordinates
-        const double tmp = Walker->at(j)-center->at(j);
+        const double tmp = Walker->node->at(j)-center->at(j);
         *rasterfile << ((fabs(tmp) < MYEPSILON) ? 0 : tmp) << " ";
       }
@@ -656 +651 @@
       for (i=0;i<3;i++) {  // print each node
         for (int j=0;j<NDIM;j++) { // and for each node all NDIM coordinates
-          const double tmp = TriangleRunner->second->endpoints[i]->node->at(j)-center->at(j);
+          const double tmp = TriangleRunner->second->endpoints[i]->node->node->at(j)-center->at(j);
           *rasterfile << ((fabs(tmp) < MYEPSILON) ? 0 : tmp) << " ";
         }
@@ -710 +705 @@
       LookupList[Walker->nr] = Counter++;
       for (int i=0;i<NDIM;i++) {
-        const double tmp = Walker->at(i);
+        const double tmp = Walker->node->at(i);
         *tecplot << ((fabs(tmp) < MYEPSILON) ? 0 : tmp) << " ";
       }
@@ -759 +754 @@
     TriangleSet *triangles = TesselStruct->GetAllTriangles(PointRunner->second);
     for (TriangleSet::iterator TriangleRunner = triangles->begin(); TriangleRunner != triangles->end(); ++TriangleRunner) {
-      totalarea += CalculateAreaofGeneralTriangle((*TriangleRunner)->endpoints[0]->node->getPosition() , (*TriangleRunner)->endpoints[1]->node->getPosition() , (*TriangleRunner)->endpoints[2]->node->getPosition());
+      totalarea += CalculateAreaofGeneralTriangle(*(*TriangleRunner)->endpoints[0]->node->node, *(*TriangleRunner)->endpoints[1]->node->node, *(*TriangleRunner)->endpoints[2]->node->node);
     }
     ConcavityPerLine *= totalarea;
@@ -768 +763 @@
       line = (*TriangleRunner)->GetThirdLine(PointRunner->second);
       triangle = line->GetOtherTriangle(*TriangleRunner);
-      area = CalculateAreaofGeneralTriangle(triangle->endpoints[0]->node->getPosition() , triangle->endpoints[1]->node->getPosition() , triangle->endpoints[2]->node->getPosition());
-      area += CalculateAreaofGeneralTriangle((*TriangleRunner)->endpoints[0]->node->getPosition() , (*TriangleRunner)->endpoints[1]->node->getPosition() , (*TriangleRunner)->endpoints[2]->node->getPosition());
+      area = CalculateAreaofGeneralTriangle(*triangle->endpoints[0]->node->node, *triangle->endpoints[1]->node->node, *triangle->endpoints[2]->node->node);
+      area += CalculateAreaofGeneralTriangle(*(*TriangleRunner)->endpoints[0]->node->node, *(*TriangleRunner)->endpoints[1]->node->node, *(*TriangleRunner)->endpoints[2]->node->node);
       area *= -line->CalculateConvexity();
       if (area > 0)
@@ -802 +797 @@
     distance = 0.;
     for (TriangleMap::const_iterator TriangleRunner = Convex->TrianglesOnBoundary.begin(); TriangleRunner != Convex->TrianglesOnBoundary.end(); TriangleRunner++) {
-      const double CurrentDistance = Convex->GetDistanceSquaredToTriangle(PointRunner->second->node->getPosition() , TriangleRunner->second);
+      const double CurrentDistance = Convex->GetDistanceSquaredToTriangle(*PointRunner->second->node->node, TriangleRunner->second);
       if (CurrentDistance < distance)
         distance = CurrentDistance;