Implemented function for determining the two endpoints of a given midpoint for the 3D cases.

Author: MarkoSterbentz

src/axom/mir/CellGenerator.cpp

@@ -110,7 +110,7 @@ void  CellGenerator::generateVertexVolumeFractions(const mir::Shape shapeType,
   {
     int vID = itr->first;
 
-    for (int matID = 0; matID < vertexVF.size(); ++matID)
+    for (unsigned long matID = 0; matID < vertexVF.size(); ++matID)
     {
       if ( vID < mir::utilities::numVerts(shapeType) )
       {

src/axom/mir/MIRUtilities.hpp

@@ -177,6 +177,80 @@ namespace utilities
         if ( midpointVertexID == 7 && isFromVertex  ) { return 3; }
         if ( midpointVertexID == 7 && !isFromVertex ) { return 0; }
         break;
+      case mir::Shape::Tetrahedron:
+        if ( midpointVertexID == 4 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID == 4 && !isFromVertex ) { return 1; }
+        if ( midpointVertexID == 5 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID == 5 && !isFromVertex ) { return 2; } 
+        if ( midpointVertexID == 6 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID == 6 && !isFromVertex ) { return 3; }
+        if ( midpointVertexID == 7 && isFromVertex  ) { return 1; }
+        if ( midpointVertexID == 7 && !isFromVertex ) { return 2; }
+        if ( midpointVertexID == 8 && isFromVertex  ) { return 2; }
+        if ( midpointVertexID == 8 && !isFromVertex ) { return 3; }
+        if ( midpointVertexID == 9 && isFromVertex  ) { return 3; }
+        if ( midpointVertexID == 9 && !isFromVertex ) { return 1; }
+      case mir::Shape::Pyramid:
+        if ( midpointVertexID ==  5 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID ==  5 && !isFromVertex ) { return 1; } 
+        if ( midpointVertexID ==  6 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID ==  6 && !isFromVertex ) { return 2; }
+        if ( midpointVertexID ==  7 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID ==  7 && !isFromVertex ) { return 3; }
+        if ( midpointVertexID ==  8 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID ==  8 && !isFromVertex ) { return 4; }
+        if ( midpointVertexID ==  9 && isFromVertex  ) { return 1; }
+        if ( midpointVertexID ==  9 && !isFromVertex ) { return 2; }
+        if ( midpointVertexID == 10 && isFromVertex  ) { return 2; }
+        if ( midpointVertexID == 10 && !isFromVertex ) { return 3; }
+        if ( midpointVertexID == 11 && isFromVertex  ) { return 3; }
+        if ( midpointVertexID == 11 && !isFromVertex ) { return 4; }
+        if ( midpointVertexID == 12 && isFromVertex  ) { return 4; }
+        if ( midpointVertexID == 12 && !isFromVertex ) { return 1; }
+      case mir::Shape::Triangular_Prism:
+        if ( midpointVertexID ==  6 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID ==  6 && !isFromVertex ) { return 1; }
+        if ( midpointVertexID ==  7 && isFromVertex  ) { return 1; }
+        if ( midpointVertexID ==  7 && !isFromVertex ) { return 2; }
+        if ( midpointVertexID ==  8 && isFromVertex  ) { return 2; }
+        if ( midpointVertexID ==  8 && !isFromVertex ) { return 0; }
+        if ( midpointVertexID ==  9 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID ==  9 && !isFromVertex ) { return 3; }
+        if ( midpointVertexID == 10 && isFromVertex  ) { return 1; }
+        if ( midpointVertexID == 10 && !isFromVertex ) { return 4; }
+        if ( midpointVertexID == 11 && isFromVertex  ) { return 2; }
+        if ( midpointVertexID == 11 && !isFromVertex ) { return 5; }
+        if ( midpointVertexID == 12 && isFromVertex  ) { return 3; }
+        if ( midpointVertexID == 12 && !isFromVertex ) { return 4; }
+        if ( midpointVertexID == 13 && isFromVertex  ) { return 4; }
+        if ( midpointVertexID == 13 && !isFromVertex ) { return 5; }
+        if ( midpointVertexID == 14 && isFromVertex  ) { return 5; }
+        if ( midpointVertexID == 14 && !isFromVertex ) { return 3; }
+      case mir::Shape::Hexahedron:
+        if ( midpointVertexID ==  8 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID ==  8 && !isFromVertex ) { return 1; }
+        if ( midpointVertexID ==  9 && isFromVertex  ) { return 1; }
+        if ( midpointVertexID ==  9 && !isFromVertex ) { return 2; }
+        if ( midpointVertexID == 10 && isFromVertex  ) { return 2; }
+        if ( midpointVertexID == 10 && !isFromVertex ) { return 3; }
+        if ( midpointVertexID == 11 && isFromVertex  ) { return 3; }
+        if ( midpointVertexID == 11 && !isFromVertex ) { return 0; }
+        if ( midpointVertexID == 12 && isFromVertex  ) { return 0; }
+        if ( midpointVertexID == 12 && !isFromVertex ) { return 4; }
+        if ( midpointVertexID == 13 && isFromVertex  ) { return 1; }
+        if ( midpointVertexID == 13 && !isFromVertex ) { return 5; }
+        if ( midpointVertexID == 14 && isFromVertex  ) { return 2; }
+        if ( midpointVertexID == 14 && !isFromVertex ) { return 6; }
+        if ( midpointVertexID == 15 && isFromVertex  ) { return 3; }
+        if ( midpointVertexID == 15 && !isFromVertex ) { return 7; }
+        if ( midpointVertexID == 16 && isFromVertex  ) { return 4; }
+        if ( midpointVertexID == 16 && !isFromVertex ) { return 5; }
+        if ( midpointVertexID == 17 && isFromVertex  ) { return 5; }
+        if ( midpointVertexID == 17 && !isFromVertex ) { return 6; }
+        if ( midpointVertexID == 18 && isFromVertex  ) { return 6; }
+        if ( midpointVertexID == 18 && !isFromVertex ) { return 7; }
+        if ( midpointVertexID == 19 && isFromVertex  ) { return 7; }
+        if ( midpointVertexID == 19 && !isFromVertex ) { return 4; }
       default:
         printf("Edge endpoint case not implemented.\n");
         return -1;

src/axom/mir/examples/mir_tutorial_simple.cpp

@@ -48,7 +48,7 @@ int main( int AXOM_NOT_USED(argc), char** AXOM_NOT_USED(argv) )
   std::cout << "Material interface reconstruction time: " << std::chrono::duration_cast<std::chrono::milliseconds>(endTime - startTime).count() << " ms" << std::endl;
 
   // Output results
-  processedMesh.writeMeshToFile("/Users/sterbentz3/Desktop/testOutputNewMARKONEW4.vtk");
+  processedMesh.writeMeshToFile("/Users/sterbentz3/Desktop/processedMesh.vtk");
 
   std::vector<std::vector<axom::float64> > materialVolumeFractionsElement = processedMesh.computeOriginalElementVolumeFractions();