[MOAB-dev] r4812 - in MOAB/trunk: src src/moab test

[iulian at mcs.anl.gov]

Author: iulian
Date: 2011-05-10 16:06:34 -0500 (Tue, 10 May 2011)
New Revision: 4812

Modified:
   MOAB/trunk/src/FBEngine.cpp
   MOAB/trunk/src/moab/FBEngine.hpp
   MOAB/trunk/test/mbfacet_test.cpp
Log:
fix bugs in splitting facet-based geometry faces
set neumann tags (mostly for debugging)
start filling for the 2 faces from the splitting line
first range is to the left (positive sense) , second range to the right
(negative sense)



Modified: MOAB/trunk/src/FBEngine.cpp
===================================================================
--- MOAB/trunk/src/FBEngine.cpp	2011-05-10 17:06:06 UTC (rev 4811)
+++ MOAB/trunk/src/FBEngine.cpp	2011-05-10 21:06:34 UTC (rev 4812)
@@ -29,7 +29,7 @@
 double tolerance = 0.01; // TODO: how is this used ????
 double tolerance_segment = 0.000001; // for segments intersection, points collapse
 const bool Debug_surf_eval = false;
-bool debug_splits = false;// false?
+bool debug_splits = false;
 
 // will compute intersection between a segment and slice of a plane
 // output is the intersection point
@@ -1336,53 +1336,10 @@
   _my_geomTopoTool->add_geo_set(new_geo_edge, 1);
   for (unsigned int k=0; k< geo_vertices.size(); k++)
     _my_geomTopoTool->add_geo_set(geo_vertices[k], 0);
-  // new face will have the new edge oriented positively
-  // get a mesh edge from geo edge;
 
-  Range mesh_edges;
-  rval = _mbImpl->get_entities_by_type(new_geo_edge, MBEDGE, mesh_edges);
-  MBERRORR(rval, "can't get new polyline edges");
-
-  // get a positive triangle adjacent to mesh_edge[0]
-  Range adj_tri;
-  EntityHandle firstMeshEdge = mesh_edges[0];
-  rval =  _mbImpl->get_adjacencies(&firstMeshEdge, 1,
-          2, false, adj_tri, Interface::UNION);
-  MBERRORR(rval, "can't get adj_tris to first edge");
-  EntityHandle pTriangle;
-  Range::iterator it=adj_tri.begin();
-  for ( ; it!=adj_tri.end(); it++)
-  {
-    EntityHandle tr=*it;
-    int num1, sense, offset;
-    rval = _mbImpl->side_number(tr, mesh_edges[0], num1, sense, offset);
-    MBERRORR(rval, "edge not adjacent");
-    if (sense==1)
-    {
-      pTriangle = tr;
-      break;
-    }
-  }
-  if (it==adj_tri.end())
-  {
-    MBERRORR(MB_FAILURE, "can't get a positive triangle");
-  }
-  bool firstIsNewFace = true;
-  if ( first.find(pTriangle)!=first.end())
-  {