[MOAB-dev] r3997 - MOAB/trunk/tools/dagmc

[wilsonp at engr.wisc.edu]

Author: wilsonp
Date: 2010-06-04 12:56:14 -0500 (Fri, 04 Jun 2010)
New Revision: 3997

Modified:
   MOAB/trunk/tools/dagmc/DagMC.cpp
Log:
DagMC::point_in_volume update:
When searching for all the triangles within epsilon of a current point that was previously found to be IN a given triangle, ensure that that list includes that given triangle.



Modified: MOAB/trunk/tools/dagmc/DagMC.cpp
===================================================================
--- MOAB/trunk/tools/dagmc/DagMC.cpp	2010-06-04 16:58:37 UTC (rev 3996)
+++ MOAB/trunk/tools/dagmc/DagMC.cpp	2010-06-04 17:56:14 UTC (rev 3997)
@@ -617,6 +617,13 @@
   rval = obbTree.sphere_intersect_triangles( closest.array(), epsilon, root, tris, &surfs );
   if (MB_SUCCESS != rval) return rval;
   
+    // Make sure this list of triangles includes the facet of this closest point
+  if(tris.size() == 0 || 
+     find(tris.begin(),tris.end(),facet) == tris.end() ) {
+    tris.push_back(facet);
+    surfs.push_back(surface);
+  }
+
     // One-triangle case : determine if point is above or below triangle
   const EntityHandle* conn;
   int len, sense;