[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