Changeset 27459a for molecuilder/src/boundary.cpp

Timestamp:

Aug 18, 2009, 8:42:39 AM (16 years ago)

Author:

Frederik Heber <heber@…>

Children:

daf5d6

Parents:

3de1d2

Message:

Further bugfixes to three-step-procedure in convex envelope, not yet working.

• new functions StoreTrianglesinFile() and CalculateConcavityPerBoundaryPoint
• new functions Tesselation::PrintAll...() to print all points, lines and triangles currently on the boundary (simply to remove commented-out code)
• FIX: Tesselation::IsConvexRectangle(): just needs one closest point and simply checks the sign of the scalar product with both endpoints to determine in-/outside
• BUGFIX: Tesselation::PickFarthestofTwoBaselines(): NewOffset was wrong by the offset to the plane

File:

• molecuilder/src/boundary.cpp (modified) (7 diffs)

molecuilder/src/boundary.cpp

@@ -620 +620 @@
   }
 
+  //CalculateConcavityPerBoundaryPoint(out, TesselStruct);
+  StoreTrianglesinFile(out, mol, filename, "-first");
+
   // First step: RemovePointFromTesselatedSurface
   PointRunner = TesselStruct->PointsOnBoundary.begin();
@@ -641 +644 @@
   }
 
-
-//  // print all lines
-//  LineRunner = TesselStruct->LinesOnBoundary.begin();
-//  LineAdvance = LineRunner;  // we need an advanced line, as the LineRunner might get removed
-//  *out << Verbose(1) << "Printing all boundary lines for debugging." << endl;
-//  while (LineRunner != TesselStruct->LinesOnBoundary.end()) {
-//    LineAdvance++;
-//    line = LineRunner->second;
-//    *out << Verbose(2) << "INFO: Current line is " << *line << "." << endl;
-//    if (LineAdvance != TesselStruct->LinesOnBoundary.end())
-//      *out << Verbose(2) << "INFO: Next line will be " << *(LineAdvance->second) << "." << endl;
-//    LineRunner = LineAdvance;
-//  }
-//
-//  // print all triangles
-//  TriangleRunner = TesselStruct->TrianglesOnBoundary.begin();
-//  TriangleAdvance = TriangleRunner;  // we need an advanced line, as the LineRunner might get removed
-//  *out << Verbose(1) << "Printing all boundary triangles for debugging." << endl;
-//  while (TriangleRunner != TesselStruct->TrianglesOnBoundary.end()) {
-//    TriangleAdvance++;
-//    *out << Verbose(2) << "INFO: Current triangle is " << *(TriangleRunner->second) << "." << endl;
-//    if (TriangleAdvance != TesselStruct->TrianglesOnBoundary.end())
-//      *out << Verbose(2) << "INFO: Next triangle will be " << *(TriangleAdvance->second) << "." << endl;
-//    TriangleRunner = TriangleAdvance;
-//  }
+  //CalculateConcavityPerBoundaryPoint(out, TesselStruct);
+  StoreTrianglesinFile(out, mol, filename, "-second");
 
   // second step: PickFarthestofTwoBaselines
@@ -674 +654 @@
     line = LineRunner->second;
     *out << Verbose(1) << "INFO: Picking farthest baseline for line is " << *line << "." << endl;
-    if (LineAdvance != TesselStruct->LinesOnBoundary.end())
     // take highest of both lines
-    TesselStruct->PickFarthestofTwoBaselines(out, line);
+    if (TesselStruct->IsConvexRectangle(out, line) == NULL)
+      TesselStruct->PickFarthestofTwoBaselines(out, line);
     LineRunner = LineAdvance;
   }
 
-  // calculate remaining concavity
-  for (PointRunner = TesselStruct->PointsOnBoundary.begin(); PointRunner != TesselStruct->PointsOnBoundary.end(); PointRunner++) {
-    point = PointRunner->second;
-    *out << Verbose(1) << "INFO: Current point is " << *point << "." << endl;
-    point->value = 0;
-    for (LineMap::iterator LineRunner = point->lines.begin(); LineRunner != point->lines.end(); LineRunner++) {
-      line = LineRunner->second;
-      *out << Verbose(2) << "INFO: Current line of point " << *point << " is " << *line << "." << endl;
-      if (!line->CheckConvexityCriterion(out))
-        point->value += 1;
-    }
-  }
-
-  // 4. Store triangles in tecplot file
-  if (filename != NULL) {
-    if (DoTecplotOutput) {
-      string OutputName(filename);
-      OutputName.append("_intermed");
-      OutputName.append(TecplotSuffix);
-      ofstream *tecplot = new ofstream(OutputName.c_str());
-      write_tecplot_file(out, tecplot, mol->TesselStruct, mol, 0);
-      tecplot->close();
-      delete(tecplot);
-    }
-    if (DoRaster3DOutput) {
-      string OutputName(filename);
-      OutputName.append("_intermed");
-      OutputName.append(Raster3DSuffix);
-      ofstream *rasterplot = new ofstream(OutputName.c_str());
-      write_raster3d_file(out, rasterplot, mol->TesselStruct, mol);
-      rasterplot->close();
-      delete(rasterplot);
-    }
-  }
+  //CalculateConcavityPerBoundaryPoint(out, TesselStruct);
+  StoreTrianglesinFile(out, mol, filename, "-third");
 
   // third step: IsConvexRectangle
@@ -735 +683 @@
   }
 
+  CalculateConcavityPerBoundaryPoint(out, TesselStruct);
+
+  // end
+  *out << Verbose(0) << "End of ConvexizeNonconvexEnvelope" << endl;
+  return volume;
+};
+
+/** Calculates the concavity for each of the BoundaryPointSet's in a Tesselation.
+ * Sets BoundaryPointSet::value equal to the number of connected lines that are not convex.
+ * \param *out output stream for debugging
+ * \param *TesselStruct pointer to Tesselation structure
+ */
+void CalculateConcavityPerBoundaryPoint(ofstream *out, class Tesselation *TesselStruct)
+{
+  class BoundaryPointSet *point = NULL;
+  class BoundaryLineSet *line = NULL;
   // calculate remaining concavity
-  for (PointRunner = TesselStruct->PointsOnBoundary.begin(); PointRunner != TesselStruct->PointsOnBoundary.end(); PointRunner++) {
+  for (PointMap::iterator PointRunner = TesselStruct->PointsOnBoundary.begin(); PointRunner != TesselStruct->PointsOnBoundary.end(); PointRunner++) {
     point = PointRunner->second;
     *out << Verbose(1) << "INFO: Current point is " << *point << "." << endl;
@@ -747 +711 @@
     }
   }
-
+};
+
+/** Stores triangles to file.
+ * \param *out output stream for debugging
+ * \param *mol molecule with atoms and bonds
+ * \param *filename prefix of filename
+ * \param *extraSuffix intermediate suffix
+ */
+void StoreTrianglesinFile(ofstream *out, molecule *mol, const char *filename, const char *extraSuffix)
+{
   // 4. Store triangles in tecplot file
   if (filename != NULL) {
     if (DoTecplotOutput) {
       string OutputName(filename);
+      OutputName.append(extraSuffix);
       OutputName.append(TecplotSuffix);
       ofstream *tecplot = new ofstream(OutputName.c_str());
@@ -760 +734 @@
     if (DoRaster3DOutput) {
       string OutputName(filename);
+      OutputName.append(extraSuffix);
       OutputName.append(Raster3DSuffix);
       ofstream *rasterplot = new ofstream(OutputName.c_str());
@@ -767 +742 @@
     }
   }
-
-  // end
-  *out << Verbose(0) << "End of ConvexizeNonconvexEnvelope" << endl;
-  return volume;
 };