[occt.git] / src / TopOpeBRepBuild / TopOpeBRepBuild_BuildEdges.cxx

// Created on: 1993-06-14
// Created by: Jean Yves LEBEY
// Copyright (c) 1993-1999 Matra Datavision
// Copyright (c) 1999-2014 OPEN CASCADE SAS
//
// This file is part of Open CASCADE Technology software library.
//
// This library is free software; you can redistribute it and/or modify it under
// the terms of the GNU Lesser General Public License version 2.1 as published
// by the Free Software Foundation, with special exception defined in the file
// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
// distribution for complete text of the license and disclaimer of any warranty.
//
// Alternatively, this file may be used under the terms of Open CASCADE
// commercial license or contractual agreement.


#include <gp_Pnt.hxx>
#include <Standard_NoSuchObject.hxx>
#include <TCollection_AsciiString.hxx>
#include <TColStd_HArray1OfInteger.hxx>
#include <TopoDS.hxx>
#include <TopoDS_Edge.hxx>
#include <TopoDS_Face.hxx>
#include <TopoDS_Shape.hxx>
#include <TopoDS_Vertex.hxx>
#include <TopOpeBRepBuild_Builder.hxx>
#include <TopOpeBRepBuild_define.hxx>
#include <TopOpeBRepBuild_EdgeBuilder.hxx>
#include <TopOpeBRepBuild_FaceBuilder.hxx>
#include <TopOpeBRepBuild_GTopo.hxx>
#include <TopOpeBRepBuild_HBuilder.hxx>
#include <TopOpeBRepBuild_PaveSet.hxx>
#include <TopOpeBRepBuild_ShapeSet.hxx>
#include <TopOpeBRepBuild_ShellFaceSet.hxx>
#include <TopOpeBRepBuild_SolidBuilder.hxx>
#include <TopOpeBRepBuild_WireEdgeSet.hxx>
#include <TopOpeBRepDS_BuildTool.hxx>
#include <TopOpeBRepDS_Curve.hxx>
#include <TopOpeBRepDS_CurveExplorer.hxx>
#include <TopOpeBRepDS_CurveIterator.hxx>
#include <TopOpeBRepDS_HDataStructure.hxx>
#include <TopOpeBRepDS_PointIterator.hxx>
#include <TopOpeBRepDS_SurfaceIterator.hxx>
#include <TopOpeBRepTool_ShapeExplorer.hxx>

#ifdef OCCT_DEBUG
extern Standard_Boolean TopOpeBRepBuild_GettraceCU();
#endif

//=======================================================================
//function : BuildEdges
//purpose  : 
//=======================================================================
void TopOpeBRepBuild_Builder::BuildEdges(const Standard_Integer iC,const Handle(TopOpeBRepDS_HDataStructure)& HDS)
{
#ifdef OCCT_DEBUG
  if (TopOpeBRepBuild_GettraceCU()) cout<<"\nBuildEdges on C "<<iC<<endl;
#endif
  const TopOpeBRepDS_Curve& C = HDS->Curve(iC);
  const Handle(Geom_Curve)& C3D = C.Curve(); 
  const Handle(TopOpeBRepDS_Interference)& I1 = C.GetSCI1();
  const Handle(TopOpeBRepDS_Interference)& I2 = C.GetSCI2();
  Standard_Boolean nnn = C3D.IsNull() && I1.IsNull() && I2.IsNull();
  if (nnn) {
    return;
  }
  
  TopoDS_Shape anEdge;
  const TopOpeBRepDS_Curve& curC = HDS->Curve(iC);
  myBuildTool.MakeEdge(anEdge,curC,HDS->DS());
  TopOpeBRepBuild_PaveSet PVS(anEdge);
  TopOpeBRepDS_PointIterator CPIT(HDS->CurvePoints(iC));
  FillVertexSet(CPIT,TopAbs_IN,PVS);
  TopOpeBRepBuild_PaveClassifier VCL(anEdge);
  Standard_Boolean equalpar = PVS.HasEqualParameters();
  if (equalpar) VCL.SetFirstParameter(PVS.EqualParameters());
  Standard_Boolean closvert = PVS.ClosedVertices();
  VCL.ClosedVertices(closvert);  
  PVS.InitLoop();
  if ( !PVS.MoreLoop() ) {
    return;
  }
  TopOpeBRepBuild_EdgeBuilder EDBU(PVS,VCL);
  TopTools_ListOfShape& EL = ChangeNewEdges(iC);
  MakeEdges(anEdge,EDBU,EL);
  TopTools_ListIteratorOfListOfShape It(EL);
  Standard_Integer inewC = -1;
  for (; It.More(); It.Next()) {
    TopoDS_Edge& newEdge = TopoDS::Edge(It.Value());
    myBuildTool.RecomputeCurves(curC,TopoDS::Edge(anEdge),newEdge,inewC,HDS);
    if (inewC != -1) ChangeNewEdges(inewC).Append(newEdge);
  }
  if (inewC != -1) {
    HDS->RemoveCurve(iC);
  }
  else {
    for (It.Initialize(EL);It.More();It.Next()) {
      TopoDS_Edge& newEdge = TopoDS::Edge(It.Value());
      myBuildTool.UpdateEdge(anEdge,newEdge);
    }
  }
}

//=======================================================================
//function : BuildEdges
//purpose  : 
//=======================================================================
void TopOpeBRepBuild_Builder::BuildEdges(const Handle(TopOpeBRepDS_HDataStructure)& HDS)
{
  TopOpeBRepDS_DataStructure& BDS = HDS->ChangeDS();

  myNewEdges.Clear();
  TopOpeBRepDS_CurveExplorer cex;

  Standard_Integer ick = 0;
  for (cex.Init(BDS,Standard_False); cex.More(); cex.Next()) {
    Standard_Integer ic = cex.Index();
    Standard_Boolean ck = cex.IsCurveKeep(ic);
    Standard_Integer im = cex.Curve(ic).Mother();
    if (ck == 1 && im != 0 && ick == 0) { ick = ic; break; }
  }
  if (ick) {
    for (cex.Init(BDS,Standard_True); cex.More(); cex.Next()) {
      Standard_Integer ic = cex.Index();
      BDS.RemoveCurve(ic);
    }
    BDS.ChangeNbCurves(ick-1);
  }

  for (cex.Init(BDS,Standard_False); cex.More(); cex.Next()) {
    Standard_Integer ic = cex.Index();
    Standard_Integer im = cex.Curve(ic).Mother();
    if (im != 0 ) continue;
    BuildEdges(ic,HDS);
  }

  Standard_Integer ip, np = HDS->NbPoints();
  TColStd_HArray1OfInteger tp(0,np,0);
  for (cex.Init(BDS); cex.More(); cex.Next()) {
#ifdef OCCT_DEBUG
//    const TopOpeBRepDS_Curve& C = cex.Curve();
#endif
    Standard_Integer ic = cex.Index();
    TopOpeBRepDS_ListIteratorOfListOfInterference it(BDS.CurveInterferences(ic));
    for(;it.More();it.Next()) {
      const Handle(TopOpeBRepDS_Interference)& I = it.Value();
      {
        Standard_Integer ig = I->Geometry();
        TopOpeBRepDS_Kind kg = I->GeometryType();
        if (kg == TopOpeBRepDS_POINT && ig <= np) tp.ChangeValue(ig) = tp.Value(ig) + 1;
      }
      {
        Standard_Integer is = I->Support();
        TopOpeBRepDS_Kind ks = I->SupportType();
        if (ks == TopOpeBRepDS_POINT) tp.ChangeValue(is) = tp.Value(is) + 1;
      }
    }
  }
  Standard_Integer is, ns = BDS.NbShapes();
  for (is = 1; is <= ns; is++ ) {
    const TopoDS_Shape& S = BDS.Shape(is);
    if(S.IsNull()) continue;
    Standard_Boolean test = (S.ShapeType() == TopAbs_EDGE); 
    if ( !test ) continue;
    TopOpeBRepDS_ListIteratorOfListOfInterference it(BDS.ShapeInterferences(is));
    for(;it.More();it.Next()) {
      const Handle(TopOpeBRepDS_Interference)& I = it.Value(); 
      {
        Standard_Integer ig = I->Geometry();
        TopOpeBRepDS_Kind kg = I->GeometryType();
        if (kg == TopOpeBRepDS_POINT) tp.ChangeValue(ig) = tp.Value(ig) + 1;
      }
      {
        Standard_Integer is1 = I->Support();
        TopOpeBRepDS_Kind ks = I->SupportType();
       if (ks == TopOpeBRepDS_POINT) tp.ChangeValue(is1) = tp.Value(is1) + 1;
      }
    }
  }
  for (ip = 1; ip <= np; ip++) if (tp.Value(ip) == 0) BDS.RemovePoint(ip);
}
Commit	Line	Data
b311480e	1	// Created on: 1993-06-14
	2	// Created by: Jean Yves LEBEY
	3	// Copyright (c) 1993-1999 Matra Datavision
973c2be1	4	// Copyright (c) 1999-2014 OPEN CASCADE SAS
b311480e	5	//
973c2be1	6	// This file is part of Open CASCADE Technology software library.
b311480e	7	//
d5f74e42	8	// This library is free software; you can redistribute it and/or modify it under
d5f74e42	9	// the terms of the GNU Lesser General Public License version 2.1 as published
973c2be1	10	// by the Free Software Foundation, with special exception defined in the file
	11	// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
	12	// distribution for complete text of the license and disclaimer of any warranty.
b311480e	13	//
973c2be1	14	// Alternatively, this file may be used under the terms of Open CASCADE
973c2be1	15	// commercial license or contractual agreement.
7fd59977	16
42cf5bc1	17
	18	#include <gp_Pnt.hxx>
	19	#include <Standard_NoSuchObject.hxx>
	20	#include <TCollection_AsciiString.hxx>
7fd59977	21	#include <TColStd_HArray1OfInteger.hxx>
7fd59977	22	#include <TopoDS.hxx>
42cf5bc1	23	#include <TopoDS_Edge.hxx>
	24	#include <TopoDS_Face.hxx>
	25	#include <TopoDS_Shape.hxx>
	26	#include <TopoDS_Vertex.hxx>
	27	#include <TopOpeBRepBuild_Builder.hxx>
7fd59977	28	#include <TopOpeBRepBuild_define.hxx>
42cf5bc1	29	#include <TopOpeBRepBuild_EdgeBuilder.hxx>
	30	#include <TopOpeBRepBuild_FaceBuilder.hxx>
	31	#include <TopOpeBRepBuild_GTopo.hxx>
	32	#include <TopOpeBRepBuild_HBuilder.hxx>
	33	#include <TopOpeBRepBuild_PaveSet.hxx>
	34	#include <TopOpeBRepBuild_ShapeSet.hxx>
	35	#include <TopOpeBRepBuild_ShellFaceSet.hxx>
	36	#include <TopOpeBRepBuild_SolidBuilder.hxx>
	37	#include <TopOpeBRepBuild_WireEdgeSet.hxx>
	38	#include <TopOpeBRepDS_BuildTool.hxx>
	39	#include <TopOpeBRepDS_Curve.hxx>
	40	#include <TopOpeBRepDS_CurveExplorer.hxx>
	41	#include <TopOpeBRepDS_CurveIterator.hxx>
	42	#include <TopOpeBRepDS_HDataStructure.hxx>
	43	#include <TopOpeBRepDS_PointIterator.hxx>
	44	#include <TopOpeBRepDS_SurfaceIterator.hxx>
	45	#include <TopOpeBRepTool_ShapeExplorer.hxx>
7fd59977	46
0797d9d3	47	#ifdef OCCT_DEBUG
1d0a9d4d	48	extern Standard_Boolean TopOpeBRepBuild_GettraceCU();
7fd59977	49	#endif
	50
	51	//=======================================================================
	52	//function : BuildEdges
	53	//purpose :
	54	//=======================================================================
	55	void TopOpeBRepBuild_Builder::BuildEdges(const Standard_Integer iC,const Handle(TopOpeBRepDS_HDataStructure)& HDS)
	56	{
0797d9d3	57	#ifdef OCCT_DEBUG
7fd59977	58	if (TopOpeBRepBuild_GettraceCU()) cout<<"\nBuildEdges on C "<<iC<<endl;
	59	#endif
	60	const TopOpeBRepDS_Curve& C = HDS->Curve(iC);
	61	const Handle(Geom_Curve)& C3D = C.Curve();
	62	const Handle(TopOpeBRepDS_Interference)& I1 = C.GetSCI1();
	63	const Handle(TopOpeBRepDS_Interference)& I2 = C.GetSCI2();
	64	Standard_Boolean nnn = C3D.IsNull() && I1.IsNull() && I2.IsNull();
	65	if (nnn) {
	66	return;
	67	}
	68
	69	TopoDS_Shape anEdge;
	70	const TopOpeBRepDS_Curve& curC = HDS->Curve(iC);
	71	myBuildTool.MakeEdge(anEdge,curC,HDS->DS());
	72	TopOpeBRepBuild_PaveSet PVS(anEdge);
	73	TopOpeBRepDS_PointIterator CPIT(HDS->CurvePoints(iC));
	74	FillVertexSet(CPIT,TopAbs_IN,PVS);
	75	TopOpeBRepBuild_PaveClassifier VCL(anEdge);
	76	Standard_Boolean equalpar = PVS.HasEqualParameters();
	77	if (equalpar) VCL.SetFirstParameter(PVS.EqualParameters());
	78	Standard_Boolean closvert = PVS.ClosedVertices();
	79	VCL.ClosedVertices(closvert);
	80	PVS.InitLoop();
	81	if ( !PVS.MoreLoop() ) {
	82	return;
	83	}
	84	TopOpeBRepBuild_EdgeBuilder EDBU(PVS,VCL);
	85	TopTools_ListOfShape& EL = ChangeNewEdges(iC);
	86	MakeEdges(anEdge,EDBU,EL);
	87	TopTools_ListIteratorOfListOfShape It(EL);
	88	Standard_Integer inewC = -1;
	89	for (; It.More(); It.Next()) {
	90	TopoDS_Edge& newEdge = TopoDS::Edge(It.Value());
	91	myBuildTool.RecomputeCurves(curC,TopoDS::Edge(anEdge),newEdge,inewC,HDS);
	92	if (inewC != -1) ChangeNewEdges(inewC).Append(newEdge);
	93	}
	94	if (inewC != -1) {
	95	HDS->RemoveCurve(iC);
	96	}
	97	else {
	98	for (It.Initialize(EL);It.More();It.Next()) {
	99	TopoDS_Edge& newEdge = TopoDS::Edge(It.Value());
	100	myBuildTool.UpdateEdge(anEdge,newEdge);
	101	}
	102	}
	103	}
	104
	105	//=======================================================================
	106	//function : BuildEdges
	107	//purpose :
	108	//=======================================================================
	109	void TopOpeBRepBuild_Builder::BuildEdges(const Handle(TopOpeBRepDS_HDataStructure)& HDS)
	110	{
	111	TopOpeBRepDS_DataStructure& BDS = HDS->ChangeDS();
	112
	113	myNewEdges.Clear();
7fd59977	114	TopOpeBRepDS_CurveExplorer cex;
	115
	116	Standard_Integer ick = 0;
	117	for (cex.Init(BDS,Standard_False); cex.More(); cex.Next()) {
	118	Standard_Integer ic = cex.Index();
	119	Standard_Boolean ck = cex.IsCurveKeep(ic);
	120	Standard_Integer im = cex.Curve(ic).Mother();
	121	if (ck == 1 && im != 0 && ick == 0) { ick = ic; break; }
	122	}
	123	if (ick) {
	124	for (cex.Init(BDS,Standard_True); cex.More(); cex.Next()) {
	125	Standard_Integer ic = cex.Index();
	126	BDS.RemoveCurve(ic);
	127	}
	128	BDS.ChangeNbCurves(ick-1);
	129	}
	130
7fd59977	131	for (cex.Init(BDS,Standard_False); cex.More(); cex.Next()) {
	132	Standard_Integer ic = cex.Index();
	133	Standard_Integer im = cex.Curve(ic).Mother();
	134	if (im != 0 ) continue;
	135	BuildEdges(ic,HDS);
	136	}
	137
	138	Standard_Integer ip, np = HDS->NbPoints();
	139	TColStd_HArray1OfInteger tp(0,np,0);
	140	for (cex.Init(BDS); cex.More(); cex.Next()) {
0797d9d3	141	#ifdef OCCT_DEBUG
7fd59977	142	// const TopOpeBRepDS_Curve& C = cex.Curve();
	143	#endif
	144	Standard_Integer ic = cex.Index();
	145	TopOpeBRepDS_ListIteratorOfListOfInterference it(BDS.CurveInterferences(ic));
	146	for(;it.More();it.Next()) {
	147	const Handle(TopOpeBRepDS_Interference)& I = it.Value();
	148	{
	149	Standard_Integer ig = I->Geometry();
	150	TopOpeBRepDS_Kind kg = I->GeometryType();
	151	if (kg == TopOpeBRepDS_POINT && ig <= np) tp.ChangeValue(ig) = tp.Value(ig) + 1;
	152	}
	153	{
	154	Standard_Integer is = I->Support();
	155	TopOpeBRepDS_Kind ks = I->SupportType();
	156	if (ks == TopOpeBRepDS_POINT) tp.ChangeValue(is) = tp.Value(is) + 1;
	157	}
	158	}
	159	}
	160	Standard_Integer is, ns = BDS.NbShapes();
	161	for (is = 1; is <= ns; is++ ) {
	162	const TopoDS_Shape& S = BDS.Shape(is);
	163	if(S.IsNull()) continue;
	164	Standard_Boolean test = (S.ShapeType() == TopAbs_EDGE);
	165	if ( !test ) continue;
	166	TopOpeBRepDS_ListIteratorOfListOfInterference it(BDS.ShapeInterferences(is));
	167	for(;it.More();it.Next()) {
	168	const Handle(TopOpeBRepDS_Interference)& I = it.Value();
	169	{
	170	Standard_Integer ig = I->Geometry();
	171	TopOpeBRepDS_Kind kg = I->GeometryType();
	172	if (kg == TopOpeBRepDS_POINT) tp.ChangeValue(ig) = tp.Value(ig) + 1;
	173	}
	174	{
	175	Standard_Integer is1 = I->Support();
	176	TopOpeBRepDS_Kind ks = I->SupportType();
	177	if (ks == TopOpeBRepDS_POINT) tp.ChangeValue(is1) = tp.Value(is1) + 1;
	178	}
	179	}
	180	}
	181	for (ip = 1; ip <= np; ip++) if (tp.Value(ip) == 0) BDS.RemovePoint(ip);
	182	}