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

[wilsonp at engr.wisc.edu]

Author: wilsonp
Date: 2010-06-01 15:13:17 -0500 (Tue, 01 Jun 2010)
New Revision: 3986

Modified:
   MOAB/trunk/tools/dagmc/DagMC.cpp
Log:
Revert boundary case tolerance in point membership to be the same as
discard distance tolerance.  This tolerance should not be greater than
the discard distance tolerance because it thus allows boundary cases
to be resolved as inside a volume and ray firing to be resolved as
outside the volume, revering parity on the ray's position.

Also added topology connections between the implicit complement's
"volume" entity handle and the list of "surface" entity handles.
Note: that this requires an update in the newcel function in MCNP's
connection code.



Modified: MOAB/trunk/tools/dagmc/DagMC.cpp
===================================================================
--- MOAB/trunk/tools/dagmc/DagMC.cpp	2010-06-01 19:04:21 UTC (rev 3985)
+++ MOAB/trunk/tools/dagmc/DagMC.cpp	2010-06-01 20:13:17 UTC (rev 3986)
@@ -591,7 +591,7 @@
 {
   ErrorCode rval;
   const double epsilon = discardDistTol;
-  const double boundary_epsilon = facetingTolerance;
+  const double boundary_epsilon = epsilon;
 
     // Get OBB Tree for volume
   assert(volume - setOffset < rootSets.size());
@@ -1806,6 +1806,11 @@
       
         // add obb root to list of obb roots
       comp_tree.insert( surf_obb_root );
+
+      // add this surf to the topology of the implicit complement volume
+      rval = MBI->add_parent_child(impl_compl_handle,*surf_i);
+      if (MB_SUCCESS != rval)
+	return rval;
     }  
   }