-
- if (!isCircleFound && anEdgeIter != aProcEdges.end() &&
- anEdgeIter->IsSame(aProcEdges.back())) { // The loop on the same edge taken twice
- aProcVertexes.pop_back();
- aProcEdges.pop_back();
- aCurVertex = aProcVertexes.back();
- aCurDir = getOuterEdgeDirection(aProcEdges.back(), aCurVertex);
- aCurNorm = aNorm.Reversed();
- continue;
- }
-
- // When the orientation is correct or the edges looped through
- // the first element, create new face and remove unnecessary edges.
- TopoDS_Face aPatch;
- createFace(*aVertIter, anEdgeIter, aProcEdges.end(), aPlane, aPatch);
- if (!aPatch.IsNull()) {
- boost::shared_ptr<GeomAPI_Shape> aFace(new GeomAPI_Shape);
- aFace->setImpl(new TopoDS_Face(aPatch));
- theResultFaces.push_back(aFace);
- }
- // push the edges used in the loop to the map
- std::list<TopoDS_Edge>::iterator anIter;
- for (anIter = anEdgeIter; anIter != aProcEdges.end(); anIter++)
- anEdgesInLoops.insert(&(*anIter));
- // remove unnecessary edges
- std::list<TopoDS_Vertex>::iterator aCopyVLoop = aVertIter;
- std::list<TopoDS_Edge>::iterator aCopyELoop = anEdgeIter;
- removeWasteEdges(aVertIter, anEdgeIter, aProcVertexes.end(), aProcEdges.end(), aMapVE);
-
- // revert the list of remaining edges
- std::list<TopoDS_Vertex> aRemainVertexes;
- for (; aVertIter != aProcVertexes.end(); aVertIter++)
- aRemainVertexes.push_front(*aVertIter);
- std::list<TopoDS_Edge> aRemainEdges;
- for (; anEdgeIter != aProcEdges.end(); anEdgeIter++)
- aRemainEdges.push_front(*anEdgeIter);
- // remove edges and vertexes used in the loop and add remaining ones
- if (aCopyVLoop != aProcVertexes.begin()) {
- aVertIter = aCopyVLoop;
- aVertIter--;
- } else
- aVertIter = aProcVertexes.end();
- aProcVertexes.erase(aCopyVLoop, aProcVertexes.end());
- aProcVertexes.insert(aProcVertexes.end(), aRemainVertexes.begin(), aRemainVertexes.end());
- if (aCopyELoop != aProcEdges.begin()) {
- anEdgeIter = aCopyELoop;
- anEdgeIter--;
- } else
- anEdgeIter = aProcEdges.end();
- aProcEdges.erase(aCopyELoop, aProcEdges.end());
- aProcEdges.insert(aProcEdges.end(), aRemainEdges.begin(), aRemainEdges.end());
-
- if (aVertIter == aProcVertexes.end())
- aVertIter = aProcVertexes.begin();
- else
- aVertIter++;
- if (anEdgeIter == aProcEdges.end())
- anEdgeIter = aProcEdges.begin();
- else
- anEdgeIter++;
- aCopyVLoop = aVertIter;
- aCopyELoop = anEdgeIter;
-
- if (aVertIter != aProcVertexes.end() &&
- aMapVE.Contains(*aVertIter) && aMapVE.FindFromKey(*aVertIter).Extent() <= 1)
- removeWasteEdges(aVertIter, anEdgeIter, aProcVertexes.end(), aProcEdges.end(), aMapVE);
- if (aCopyVLoop != aVertIter)
- aProcVertexes.erase(aCopyVLoop, aVertIter);
- if (aCopyELoop != anEdgeIter)
- aProcEdges.erase(aCopyELoop, anEdgeIter);
-
- // Check whether the next vertex already exists
- if (aVertIter != aProcVertexes.end())
- aVertIter++;
- if (aVertIter != aProcVertexes.end() && anEdgeIter != aProcEdges.end()) {
- aNextVertex = *aVertIter;
- aNextDir = getOuterEdgeDirection(*anEdgeIter, aNextVertex);
- aProcVertexes.erase(++aVertIter, aProcVertexes.end());
- aProcEdges.erase(++anEdgeIter, aProcEdges.end());
- } else {
- // Recalculate current vertex and current direction
- aProcEdges.clear();
- aProcVertexes.clear();
- if (aMapVE.Extent() > 0) {
- aNextVertex = findStartVertex(aMapVE, aDirX, aDirY);
- aProcVertexes.push_back(aNextVertex);
- }
- aNextDir = aDirY.Reversed();
- aCurNorm = aNorm.Reversed();
- }
- }
-
- // if next vertex connected only to alone edge, this is a part of wire (not a closed loop),
- // we need to go back through the list of already checked edges to find a branching vertex
- if (!aMapVE.IsEmpty() && aMapVE.Contains(aNextVertex)
- && aMapVE.FindFromKey(aNextVertex).Size() == 1) {
- std::list<TopoDS_Vertex>::reverse_iterator aVRIter = aProcVertexes.rbegin();
- std::list<TopoDS_Edge>::reverse_iterator aERIter = aProcEdges.rbegin();
- if (aVRIter != aProcVertexes.rend())
- aVRIter++;
- if (aERIter != aProcEdges.rend())
- aERIter++;
-
- for (; aERIter != aProcEdges.rend(); aERIter++, aVRIter++)
- if (aMapVE.FindFromKey(*aVRIter).Size() > 2)
- break;
- if (aERIter != aProcEdges.rend()
- || (aVRIter != aProcVertexes.rend() && aMapVE.FindFromKey(*aVRIter).Size() == 1)) { // the branching vertex was found or current list of edges is a wire without branches
- std::list<TopoDS_Edge>::iterator aEIter;
- TopoDS_Edge aCurEdge;
- if (aERIter != aProcEdges.rend()) {
- aEIter = aERIter.base();
- aCurEdge = *aERIter;
- } else
- aEIter = aProcEdges.begin();
- std::list<TopoDS_Wire> aTail;
- createWireList(*aVRIter, aEIter, aProcEdges.end(), anEdgesInLoops, aTail);
- std::list<TopoDS_Wire>::const_iterator aTailIter = aTail.begin();
- for (; aTailIter != aTail.end(); aTailIter++)
- if (!aTailIter->IsNull()) {
- boost::shared_ptr<GeomAPI_Shape> aWire(new GeomAPI_Shape);
- aWire->setImpl(new TopoDS_Shape(*aTailIter));
- theResultWires.push_back(aWire);
- }
- std::list<TopoDS_Vertex>::iterator aVIter = aVRIter.base();
- std::list<TopoDS_Edge>::iterator aEItCopy = aEIter;
- removeWasteEdges(--aVIter, aEItCopy, aProcVertexes.end(), aProcEdges.end(), aMapVE);
-
- aProcEdges.erase(aEIter, aProcEdges.end());
- aVIter = aVRIter.base();
- aProcVertexes.erase(aVIter, aProcVertexes.end());
-
- if (!aProcVertexes.empty()) {
- aNextVertex = aProcVertexes.back();
- if (!aCurEdge.IsNull())
- aNextDir = getOuterEdgeDirection(aCurEdge, aNextVertex);
- }
- } else { // there is no branching vertex in the list of proceeded,
- // so we should revert the list and go the opposite way
- aProcVertexes.reverse();
- aProcEdges.reverse();
- aNextVertex = aProcVertexes.back();
- aNextDir =
- aProcEdges.empty() ? aDirY : getOuterEdgeDirection(aProcEdges.back(), aNextVertex);
- }
- }
-
- // When all edges of current component of connectivity are processed,
- // we should repeat procedure for finding initial vertex in the remaining
- if (!aMapVE.IsEmpty() && !aMapVE.Contains(aNextVertex)) {
- aProcVertexes.clear();
- aProcEdges.clear();
-
- aStartVertex = findStartVertex(aMapVE, aDirX, aDirY);
- aProcVertexes.push_back(aStartVertex);
- aNextVertex = aStartVertex;
- aNextDir = aDirY.Reversed();
- aCurNorm = aNorm.Reversed();
- }
-
- aCurVertex = aNextVertex;
- aCurDir = aNextDir;
- }
-
- if (theResultFaces.size() > 1)
- fixIntersections(theResultFaces);
-}
-
-void GeomAlgoAPI_SketchBuilder::fixIntersections(
- std::list<boost::shared_ptr<GeomAPI_Shape> >& theFaces)
-{
- BRepClass_FaceClassifier aClassifier;
-
- std::list<boost::shared_ptr<GeomAPI_Shape> >::iterator anIter1 = theFaces.begin();
- std::list<boost::shared_ptr<GeomAPI_Shape> >::iterator anIter2;
- for (; anIter1 != theFaces.end(); anIter1++) {
- anIter2 = anIter1;
- for (++anIter2; anIter2 != theFaces.end(); anIter2++) {
- const TopoDS_Face& aF1 = (*anIter1)->impl<TopoDS_Face>();
- assert(aF1.ShapeType() == TopAbs_FACE); // all items in result list should be faces
- TopExp_Explorer aVert2((*anIter2)->impl<TopoDS_Shape>(), TopAbs_VERTEX);
- for (; aVert2.More(); aVert2.Next()) {
- const TopoDS_Vertex& aV = (const TopoDS_Vertex&)aVert2.Current();
- aClassifier.Perform(aF1, BRep_Tool::Pnt(aV), tolerance);
- TopAbs_State aState = aClassifier.State();
- if (aState != TopAbs_IN && aState != TopAbs_ON)
- break;
- }
- if (aVert2.More()) { // second shape is not inside first, change the shapes order and repeat comparision
- const TopoDS_Face& aF2 = (*anIter2)->impl<TopoDS_Face>();
- assert(aF2.ShapeType() == TopAbs_FACE); // all items in result list should be faces
- TopExp_Explorer aVert1((*anIter1)->impl<TopoDS_Shape>(), TopAbs_VERTEX);
- for (; aVert1.More(); aVert1.Next()) {
- const TopoDS_Vertex& aV = (const TopoDS_Vertex&)aVert2.Current();
- aClassifier.Perform(aF2, BRep_Tool::Pnt(aV), tolerance);
- TopAbs_State aState = aClassifier.State();
- if (aState != TopAbs_IN && aState != TopAbs_ON)
- break;
- }
- if (!aVert1.More()) { // first shape should be cut from the second
- BRepAlgoAPI_Cut aCut((*anIter2)->impl<TopoDS_Shape>(), (*anIter1)->impl<TopoDS_Shape>());
- aCut.Build();
- TopExp_Explorer anExp(aCut.Shape(), TopAbs_FACE);
- bool isFirstFace = true;
- for (; anExp.More(); anExp.Next()) {
- if (anExp.Current().ShapeType() != TopAbs_FACE) continue;
- if (isFirstFace) {
- (*anIter2)->setImpl(new TopoDS_Shape(anExp.Current()));
- isFirstFace = false;
- } else {
- boost::shared_ptr<GeomAPI_Shape> aShape(new GeomAPI_Shape);
- aShape->setImpl(new TopoDS_Shape(anExp.Current()));
- theFaces.push_back(aShape);
- }
- }
- }
- } else { // second shape should be cut from the first
- BRepAlgoAPI_Cut aCut((*anIter1)->impl<TopoDS_Shape>(), (*anIter2)->impl<TopoDS_Shape>());
- aCut.Build();
- TopExp_Explorer anExp(aCut.Shape(), TopAbs_FACE);
- bool isFirstFace = true;
- for (; anExp.More(); anExp.Next()) {
- if (anExp.Current().ShapeType() != TopAbs_FACE) continue;
- if (isFirstFace) {
- (*anIter1)->setImpl(new TopoDS_Shape(anExp.Current()));
- isFirstFace = false;
- } else {
- boost::shared_ptr<GeomAPI_Shape> aShape(new GeomAPI_Shape);
- aShape->setImpl(new TopoDS_Shape(anExp.Current()));
- theFaces.push_back(aShape);
- }
- }
- }