Changeset cc9225 for molecuilder/src/tesselationhelpers.cpp

Timestamp:

Nov 24, 2009, 3:34:49 PM (16 years ago)

Author:

Frederik Heber <heber@…>

Children:

1a3354

Parents:

7d1ad9

Message:

Fixes and naming of final Tecplot output file is now molecule name.

• FIXES to builder.cpp:
  • case 'p' would not dissect the molecule into connected subgraphs
  • if done so, the BondGraph was not yet initialised
  • if done so, we need to check whether BondGraphFileName has been set
  • BondGraphFileName.empty() could not be used as a check, as was set to "" instead of "\n"
  • if (finally) done so, we have to remove the empty molecule that we parsed in
  • ... and then pick the newly added molecule for mol to point at
  • SaveConfig() did not set the merged molecule name correctly.
  • if empty config is given, the empty molecule now receives the ConfigFileName as name

• FIXES to config.cpp
  • Load() did not set the name of the molecule

• changes to tesselationhelper.cpp and tesselation.cpp
  • WriteTecplotFile() uses N == -1 to make clear, that we are writing the final not a temporary tecplot file and then uses the molecule's name instead of the last triangle's endpoints
  • FindNextSuitableTriangle(): HasOtherBaselineBetterCandidate() was commented out for the moment (as it is not yet final and makes all the test fail)
  • IncludeSphereinRaster3D() did not check whether LastTriangle was NULL or not

• MoleculeListClass::DissectMoleculeIntoConnectedSubgraphs() names the subgraphs with suffixes "..-1", "..-2", ...

• changes to PointCloud and molecule
  • new virtual function PointCloud::GetName() returns "unknown"
  • new function molecule::GetName() returns pointer to name of molecule
  • IsEmpty() and IsEnd() now return true by default

• all .dat files in the "Tesselations" tests were changed accordingly (i.e. have now the name in the second line)
• benzene was added as a new test for an absolutely flat molecule. So far, it fails horribly.
• 13 of 17 tests run fine

Signed-off-by: Frederik Heber <heber@…>

File:

• molecuilder/src/tesselationhelpers.cpp (modified) (2 diffs)

molecuilder/src/tesselationhelpers.cpp

@@ -756 +756 @@
 {
   Vector helper;
-  // include the current position of the virtual sphere in the temporary raster3d file
-  Vector *center = cloud->GetCenter();
-  // make the circumsphere's center absolute again
-  helper.CopyVector(Tess->LastTriangle->endpoints[0]->node->node);
-  helper.AddVector(Tess->LastTriangle->endpoints[1]->node->node);
-  helper.AddVector(Tess->LastTriangle->endpoints[2]->node->node);
-  helper.Scale(1./3.);
-  helper.SubtractVector(center);
-  // and add to file plus translucency object
-  *rasterfile << "# current virtual sphere\n";
-  *rasterfile << "8\n  25.0    0.6     -1.0 -1.0 -1.0     0.2        0 0 0 0\n";
-  *rasterfile << "2\n  " << helper.x[0] << " " << helper.x[1] << " " << helper.x[2] << "\t" << 5. << "\t1 0 0\n";
-  *rasterfile << "9\n  terminating special property\n";
-  delete(center);
+
+  if (Tess->LastTriangle != NULL) {
+    // include the current position of the virtual sphere in the temporary raster3d file
+    Vector *center = cloud->GetCenter();
+    // make the circumsphere's center absolute again
+    helper.CopyVector(Tess->LastTriangle->endpoints[0]->node->node);
+    helper.AddVector(Tess->LastTriangle->endpoints[1]->node->node);
+    helper.AddVector(Tess->LastTriangle->endpoints[2]->node->node);
+    helper.Scale(1./3.);
+    helper.SubtractVector(center);
+    // and add to file plus translucency object
+    *rasterfile << "# current virtual sphere\n";
+    *rasterfile << "8\n  25.0    0.6     -1.0 -1.0 -1.0     0.2        0 0 0 0\n";
+    *rasterfile << "2\n  " << helper.x[0] << " " << helper.x[1] << " " << helper.x[2] << "\t" << 5. << "\t1 0 0\n";
+    *rasterfile << "9\n  terminating special property\n";
+    delete(center);
+  }
 };
 
@@ -829 +832 @@
     *tecplot << "TITLE = \"3D CONVEX SHELL\"" << endl;
     *tecplot << "VARIABLES = \"X\" \"Y\" \"Z\" \"U\"" << endl;
-    *tecplot << "ZONE T=\"" << N << "-";
-    for (int i=0;i<3;i++)
-      *tecplot << (i==0 ? "" : "_") << TesselStruct->LastTriangle->endpoints[i]->node->Name;
+    *tecplot << "ZONE T=\"";
+    if (N < 0) {
+      *tecplot << cloud->GetName();
+    } else {
+      *tecplot << N << "-";
+      for (int i=0;i<3;i++)
+        *tecplot << (i==0 ? "" : "_") << TesselStruct->LastTriangle->endpoints[i]->node->Name;
+    }
     *tecplot << "\", N=" << TesselStruct->PointsOnBoundary.size() << ", E=" << TesselStruct->TrianglesOnBoundary.size() << ", DATAPACKING=POINT, ZONETYPE=FETRIANGLE" << endl;
     int i=0;