[MOAB-dev] r4813 - MOAB/trunk/src

iulian at mcs.anl.gov iulian at mcs.anl.gov
Tue May 10 16:40:19 CDT 2011


Author: iulian
Date: 2011-05-10 16:40:18 -0500 (Tue, 10 May 2011)
New Revision: 4813

Modified:
   MOAB/trunk/src/FBEngine.cpp
Log:
when an old edge is split, need to also set the
sense of the new edge in the adjacent faces
still to do: when an edge is removed from a face, we need to clean the sense tags
too


Modified: MOAB/trunk/src/FBEngine.cpp
===================================================================
--- MOAB/trunk/src/FBEngine.cpp	2011-05-10 21:06:34 UTC (rev 4812)
+++ MOAB/trunk/src/FBEngine.cpp	2011-05-10 21:40:18 UTC (rev 4813)
@@ -2374,14 +2374,20 @@
   Range faceRange;
   rval = getAdjacentEntities(edge, 2, faceRange);
 
-  // these faces will be adjacent to the new edge!
+  // these faces will be adjacent to the new edge too!
+  // need to set the senses of new edge within faces
 
   for (Range::iterator it= faceRange.begin(); it!=faceRange.end(); it++)
   {
     EntityHandle face = *it;
-    // at this point, we just need to add to both edges the new node set (vertex)
     rval = _mbImpl->add_parent_child(face, new_edge);
     MBERRORR(rval, " can't add new edge - face parent relation");
+    int sense;
+    rval = _my_geomTopoTool->get_sense(edge, face, sense);
+    MBERRORR(rval, " can't get initial sense of edge in the adjacent face");
+    // keep the same sense for the new edge within the faces
+    rval = _my_geomTopoTool->set_sense(new_edge, face, sense);
+    MBERRORR(rval, " can't set sense of new edge in the adjacent face");
   }
 
   return MB_SUCCESS;


























More information about the moab-dev mailing list