[occt.git] / src / QANewModTopOpe / QANewModTopOpe_Glue_vertex.cxx

// Created on: 2001-01-16
// Created by: Michael SAZONOV
// Copyright (c) 2001-2012 OPEN CASCADE SAS
//
// The content of this file is subject to the Open CASCADE Technology Public
// License Version 6.5 (the "License"). You may not use the content of this file
// except in compliance with the License. Please obtain a copy of the License
// at http://www.opencascade.org and read it completely before using this file.
//
// The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
// main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
//
// The Original Code and all software distributed under the License is
// distributed on an "AS IS" basis, without warranty of any kind, and the
// Initial Developer hereby disclaims all such warranties, including without
// limitation, any warranties of merchantability, fitness for a particular
// purpose or non-infringement. Please see the License for the specific terms
// and conditions governing the rights and limitations under the License.

#include <QANewModTopOpe_Glue.ixx>
#include <TopoDS.hxx>
#include <TopoDS_Iterator.hxx>
#include <BRep_Builder.hxx>
#include <BRep_Tool.hxx>
#include <TopExp.hxx>
#include <TopTools_ListIteratorOfListOfShape.hxx>
#include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
#include <BRepExtrema_DistShapeShape.hxx>
#include <TopExp_Explorer.hxx>

//=======================================================================
//function : PerformVertex
//purpose  : 
//=======================================================================

void
QANewModTopOpe_Glue::PerformVertex() 
{
  BRepExtrema_DistShapeShape aExtrema (myS1, myS2);
  if (!aExtrema.IsDone()) return;
  if (aExtrema.InnerSolution()) {
    if (myS2.ShapeType() != TopAbs_VERTEX) return;
    TopoDS_Vertex aV = TopoDS::Vertex(myS2);
    BRep_Builder aBld;
    aBld.Add(myS1,aV.Oriented(TopAbs_INTERNAL));

    TopTools_ListOfShape aList;
    aList.Append (aV);
    myMapGener.Bind(myS1, aList);
        
    myShape = myS1;
    myShape.Orientation(myS1.Orientation());
    Done();
 
  }
  else {
    Standard_Integer nbSol = aExtrema.NbSolution(), i;
    Standard_Real aDist = aExtrema.Value();
    TopoDS_Shape aDummy = myS2.Oriented(TopAbs_FORWARD);
    TopoDS_Vertex aVer2 = TopoDS::Vertex (aDummy);
    Standard_Real aTol2 = BRep_Tool::Tolerance(aVer2);
    
    TopTools_MapOfShape aMapPassed;
    aMapPassed.Add (myS2);
    TopTools_IndexedDataMapOfShapeListOfShape aMapAnc;
    TopExp::MapShapesAndAncestors (myS1, TopAbs_VERTEX, TopAbs_EDGE, aMapAnc);
    TopExp::MapShapesAndAncestors (myS1, TopAbs_VERTEX, TopAbs_FACE, aMapAnc);
    TopExp::MapShapesAndAncestors (myS1, TopAbs_EDGE, TopAbs_FACE, aMapAnc);
    BRep_Builder aBld;
    
    // pass 1: process contacted vertices
    for (i=1; i <= nbSol; i++) {
      TopoDS_Shape aShape = aExtrema.SupportOnShape1(i);
      if (aShape.ShapeType() != TopAbs_VERTEX ||
	  aMapPassed.Contains(aShape)) continue;
      
      const TopoDS_Vertex& aVer = TopoDS::Vertex (aShape);
      
      if(aVer.IsSame(aVer2)) continue;
      
      Standard_Real aTol1 = BRep_Tool::Tolerance(aVer);
      if (aDist > aTol1 && aDist > aTol2) continue;
      
      aTol2 = Max (aTol2, aTol1 + aDist);
      aBld.UpdateVertex (aVer2, aTol2);
      
      // substitute aVer2 instead of aVer
      TopTools_ListOfShape aList;
      aList.Append (aVer2);
      mySubst.Substitute (aVer, aList);
      aMapPassed.Add(aVer);
      TopTools_ListIteratorOfListOfShape aIter(aMapAnc.FindFromKey(aVer));
      for (; aIter.More(); aIter.Next()) {
	aMapPassed.Add(aIter.Value());
      }
      myMapModif.Bind(aVer, aList);
    }
    
    // pass 2: process contacted edges
    for (i=1; i <= nbSol; i++) {
      TopoDS_Shape aShape = aExtrema.SupportOnShape1(i);
      if (aShape.ShapeType() != TopAbs_EDGE ||
	  aMapPassed.Contains(aShape)) continue;
      
      const TopoDS_Edge& aEdge = TopoDS::Edge (aShape);
      Standard_Real aTol1 = BRep_Tool::Tolerance(aEdge);
      if (aDist > aTol1 && aDist > aTol2) continue;
      Standard_Real aPar;
      aExtrema.ParOnEdgeS1(i, aPar);
      
      // construct new edge
      TopoDS_Edge aNewEdge;
      InsertVertexInEdge (aEdge, aVer2, aPar, aNewEdge);
      
      // substitute edge
      TopTools_ListOfShape aList;
      aList.Append (aNewEdge.Oriented(TopAbs_FORWARD));
      mySubst.Substitute (aEdge, aList);
      aMapPassed.Add(aEdge);
      TopTools_ListIteratorOfListOfShape aIter(aMapAnc.FindFromKey(aEdge));
      for (; aIter.More(); aIter.Next()) {
	aMapPassed.Add(aIter.Value());
      }
      // for Mandrake-10 - mkv,02.06.06 - myMapModif.Bind(aEdge, TopTools_ListOfShape());
      TopTools_ListOfShape aList1;
      myMapModif.Bind(aEdge, aList1);
      myMapModif(aEdge).Append (aNewEdge);
      // for Mandrake-10 - mkv,02.06.06 - myMapGener.Bind(aEdge, TopTools_ListOfShape());
      TopTools_ListOfShape aList2;
      myMapGener.Bind(aEdge, aList2);
      myMapGener(aEdge).Append (aVer2);
    }

    // pass 3: process contacted faces
    for (i=1; i <= nbSol; i++) {
      TopoDS_Shape aShape = aExtrema.SupportOnShape1(i);
      if (aShape.ShapeType() != TopAbs_FACE ||
	  aMapPassed.Contains(aShape)) continue;
      
      const TopoDS_Face& aFace = TopoDS::Face (aShape);
      Standard_Real aTol1 = BRep_Tool::Tolerance(aFace);
      if (aDist > aTol1 && aDist > aTol2) continue;
      Standard_Real aParU, aParV;
      aExtrema.ParOnFaceS1(i, aParU, aParV);
      
      // construct new face
      TopoDS_Face aNewFace = aFace;
      aNewFace.EmptyCopy();
      aNewFace.Orientation (TopAbs_FORWARD);
      aBld.NaturalRestriction (aNewFace, BRep_Tool::NaturalRestriction(aFace));
      // add old subshapes
      TopoDS_Iterator aIterF (aFace, Standard_False);
      for (; aIterF.More(); aIterF.Next()) {
	aBld.Add (aNewFace, aIterF.Value());
      }
      // add new internal vertex
      aTol2 = Max (aTol2, aTol1);
      aBld.Add (aNewFace, aVer2.Oriented(TopAbs_INTERNAL));
      aBld.UpdateVertex (aVer2, aParU, aParV, aNewFace, aTol2);
      
      // substitute face
      TopTools_ListOfShape aList;
      aList.Append (aNewFace);
      mySubst.Substitute (aFace, aList);
      aMapPassed.Add(aFace);
      // for Mandrake-10 - mkv,02.06.06 - myMapModif.Bind(aFace, TopTools_ListOfShape());
      TopTools_ListOfShape aList3;
      myMapModif.Bind(aFace, aList3);
      myMapModif(aFace).Append (aNewFace);
      // for Mandrake-10 - mkv,02.06.06 - myMapGener.Bind(aFace, TopTools_ListOfShape());
      TopTools_ListOfShape aList4;
      myMapGener.Bind(aFace, aList4);
      myMapGener(aFace).Append (aVer2);
    }
    
    mySubst.Build(myS1);
    
    TopExp_Explorer aExp (myS1, TopAbs_FACE);
    for (; aExp.More(); aExp.Next()) {
      const TopoDS_Shape& aFace = aExp.Current();
      if(myMapModif.IsBound(aFace)) continue;
      if (mySubst.IsCopied(aFace)) {
	if (!mySubst.Copy(aFace).IsEmpty()) {	
	  myMapModif.Bind(aFace,mySubst.Copy(aFace));
	}
      }
    }
    
    aExp.Init(myS1, TopAbs_EDGE);
    for (; aExp.More(); aExp.Next()) {
      const TopoDS_Shape& anEdge = aExp.Current();
      if(myMapModif.IsBound(anEdge)) continue;
      if (mySubst.IsCopied(anEdge)) {
	if (!mySubst.Copy(anEdge).IsEmpty()) {	
	  myMapModif.Bind(anEdge,mySubst.Copy(anEdge));
	}
      }
    }

    if (mySubst.IsCopied(myS1)) {
      myShape = mySubst.Copy(myS1).First();
      myShape.Orientation(myS1.Orientation());
      Done();
    }
  }
}
Commit	Line	Data
b311480e	1	// Created on: 2001-01-16
	2	// Created by: Michael SAZONOV
	3	// Copyright (c) 2001-2012 OPEN CASCADE SAS
	4	//
	5	// The content of this file is subject to the Open CASCADE Technology Public
	6	// License Version 6.5 (the "License"). You may not use the content of this file
	7	// except in compliance with the License. Please obtain a copy of the License
	8	// at http://www.opencascade.org and read it completely before using this file.
	9	//
	10	// The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
	11	// main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
	12	//
	13	// The Original Code and all software distributed under the License is
	14	// distributed on an "AS IS" basis, without warranty of any kind, and the
	15	// Initial Developer hereby disclaims all such warranties, including without
	16	// limitation, any warranties of merchantability, fitness for a particular
	17	// purpose or non-infringement. Please see the License for the specific terms
	18	// and conditions governing the rights and limitations under the License.
7fd59977	19
	20	#include <QANewModTopOpe_Glue.ixx>
	21	#include <TopoDS.hxx>
	22	#include <TopoDS_Iterator.hxx>
	23	#include <BRep_Builder.hxx>
	24	#include <BRep_Tool.hxx>
	25	#include <TopExp.hxx>
	26	#include <TopTools_ListIteratorOfListOfShape.hxx>
	27	#include <TopTools_IndexedDataMapOfShapeListOfShape.hxx>
	28	#include <BRepExtrema_DistShapeShape.hxx>
	29	#include <TopExp_Explorer.hxx>
	30
	31	//=======================================================================
	32	//function : PerformVertex
	33	//purpose :
	34	//=======================================================================
	35
	36	void
	37	QANewModTopOpe_Glue::PerformVertex()
	38	{
	39	BRepExtrema_DistShapeShape aExtrema (myS1, myS2);
	40	if (!aExtrema.IsDone()) return;
	41	if (aExtrema.InnerSolution()) {
	42	if (myS2.ShapeType() != TopAbs_VERTEX) return;
	43	TopoDS_Vertex aV = TopoDS::Vertex(myS2);
	44	BRep_Builder aBld;
	45	aBld.Add(myS1,aV.Oriented(TopAbs_INTERNAL));
	46
	47	TopTools_ListOfShape aList;
	48	aList.Append (aV);
	49	myMapGener.Bind(myS1, aList);
	50
	51	myShape = myS1;
	52	myShape.Orientation(myS1.Orientation());
	53	Done();
	54
	55	}
	56	else {
	57	Standard_Integer nbSol = aExtrema.NbSolution(), i;
	58	Standard_Real aDist = aExtrema.Value();
	59	TopoDS_Shape aDummy = myS2.Oriented(TopAbs_FORWARD);
	60	TopoDS_Vertex aVer2 = TopoDS::Vertex (aDummy);
	61	Standard_Real aTol2 = BRep_Tool::Tolerance(aVer2);
	62
	63	TopTools_MapOfShape aMapPassed;
	64	aMapPassed.Add (myS2);
	65	TopTools_IndexedDataMapOfShapeListOfShape aMapAnc;
	66	TopExp::MapShapesAndAncestors (myS1, TopAbs_VERTEX, TopAbs_EDGE, aMapAnc);
	67	TopExp::MapShapesAndAncestors (myS1, TopAbs_VERTEX, TopAbs_FACE, aMapAnc);
	68	TopExp::MapShapesAndAncestors (myS1, TopAbs_EDGE, TopAbs_FACE, aMapAnc);
	69	BRep_Builder aBld;
	70
	71	// pass 1: process contacted vertices
	72	for (i=1; i <= nbSol; i++) {
	73	TopoDS_Shape aShape = aExtrema.SupportOnShape1(i);
	74	if (aShape.ShapeType() != TopAbs_VERTEX \|\|
	75	aMapPassed.Contains(aShape)) continue;
	76
	77	const TopoDS_Vertex& aVer = TopoDS::Vertex (aShape);
	78
	79	if(aVer.IsSame(aVer2)) continue;
	80
	81	Standard_Real aTol1 = BRep_Tool::Tolerance(aVer);
	82	if (aDist > aTol1 && aDist > aTol2) continue;
83
84	aTol2 = Max (aTol2, aTol1 + aDist);
85	aBld.UpdateVertex (aVer2, aTol2);
86
87	// substitute aVer2 instead of aVer
88	TopTools_ListOfShape aList;
89	aList.Append (aVer2);
90	mySubst.Substitute (aVer, aList);
91	aMapPassed.Add(aVer);
92	TopTools_ListIteratorOfListOfShape aIter(aMapAnc.FindFromKey(aVer));
93	for (; aIter.More(); aIter.Next()) {
94	aMapPassed.Add(aIter.Value());
95	}
96	myMapModif.Bind(aVer, aList);
97	}
98
99	// pass 2: process contacted edges
100	for (i=1; i <= nbSol; i++) {
101	TopoDS_Shape aShape = aExtrema.SupportOnShape1(i);
102	if (aShape.ShapeType() != TopAbs_EDGE \|\|
103	aMapPassed.Contains(aShape)) continue;
104
105	const TopoDS_Edge& aEdge = TopoDS::Edge (aShape);
106	Standard_Real aTol1 = BRep_Tool::Tolerance(aEdge);
107	if (aDist > aTol1 && aDist > aTol2) continue;
108	Standard_Real aPar;
109	aExtrema.ParOnEdgeS1(i, aPar);
110
111	// construct new edge
112	TopoDS_Edge aNewEdge;
113	InsertVertexInEdge (aEdge, aVer2, aPar, aNewEdge);
114
115	// substitute edge
116	TopTools_ListOfShape aList;
117	aList.Append (aNewEdge.Oriented(TopAbs_FORWARD));
118	mySubst.Substitute (aEdge, aList);
119	aMapPassed.Add(aEdge);
120	TopTools_ListIteratorOfListOfShape aIter(aMapAnc.FindFromKey(aEdge));
121	for (; aIter.More(); aIter.Next()) {
122	aMapPassed.Add(aIter.Value());
123	}
124	// for Mandrake-10 - mkv,02.06.06 - myMapModif.Bind(aEdge, TopTools_ListOfShape());
125	TopTools_ListOfShape aList1;
126	myMapModif.Bind(aEdge, aList1);
127	myMapModif(aEdge).Append (aNewEdge);
128	// for Mandrake-10 - mkv,02.06.06 - myMapGener.Bind(aEdge, TopTools_ListOfShape());
129	TopTools_ListOfShape aList2;
130	myMapGener.Bind(aEdge, aList2);
131	myMapGener(aEdge).Append (aVer2);
132	}
133
134	// pass 3: process contacted faces
135	for (i=1; i <= nbSol; i++) {
136	TopoDS_Shape aShape = aExtrema.SupportOnShape1(i);
137	if (aShape.ShapeType() != TopAbs_FACE \|\|
138	aMapPassed.Contains(aShape)) continue;
139
140	const TopoDS_Face& aFace = TopoDS::Face (aShape);
141	Standard_Real aTol1 = BRep_Tool::Tolerance(aFace);
142	if (aDist > aTol1 && aDist > aTol2) continue;
143	Standard_Real aParU, aParV;
144	aExtrema.ParOnFaceS1(i, aParU, aParV);
145
146	// construct new face
147	TopoDS_Face aNewFace = aFace;
148	aNewFace.EmptyCopy();
149	aNewFace.Orientation (TopAbs_FORWARD);
150	aBld.NaturalRestriction (aNewFace, BRep_Tool::NaturalRestriction(aFace));
151	// add old subshapes
152	TopoDS_Iterator aIterF (aFace, Standard_False);
153	for (; aIterF.More(); aIterF.Next()) {
154	aBld.Add (aNewFace, aIterF.Value());
155	}
156	// add new internal vertex
157	aTol2 = Max (aTol2, aTol1);
158	aBld.Add (aNewFace, aVer2.Oriented(TopAbs_INTERNAL));
159	aBld.UpdateVertex (aVer2, aParU, aParV, aNewFace, aTol2);
160
161	// substitute face
162	TopTools_ListOfShape aList;
163	aList.Append (aNewFace);
164	mySubst.Substitute (aFace, aList);
165	aMapPassed.Add(aFace);
166	// for Mandrake-10 - mkv,02.06.06 - myMapModif.Bind(aFace, TopTools_ListOfShape());
167	TopTools_ListOfShape aList3;
168	myMapModif.Bind(aFace, aList3);
169	myMapModif(aFace).Append (aNewFace);
170	// for Mandrake-10 - mkv,02.06.06 - myMapGener.Bind(aFace, TopTools_ListOfShape());
171	TopTools_ListOfShape aList4;
172	myMapGener.Bind(aFace, aList4);
173	myMapGener(aFace).Append (aVer2);
174	}
175
176	mySubst.Build(myS1);
177
178	TopExp_Explorer aExp (myS1, TopAbs_FACE);
179	for (; aExp.More(); aExp.Next()) {
180	const TopoDS_Shape& aFace = aExp.Current();
181	if(myMapModif.IsBound(aFace)) continue;
182	if (mySubst.IsCopied(aFace)) {
183	if (!mySubst.Copy(aFace).IsEmpty()) {
184	myMapModif.Bind(aFace,mySubst.Copy(aFace));
185	}
186	}
187	}
188
189	aExp.Init(myS1, TopAbs_EDGE);
190	for (; aExp.More(); aExp.Next()) {
191	const TopoDS_Shape& anEdge = aExp.Current();
192	if(myMapModif.IsBound(anEdge)) continue;
193	if (mySubst.IsCopied(anEdge)) {
194	if (!mySubst.Copy(anEdge).IsEmpty()) {
195	myMapModif.Bind(anEdge,mySubst.Copy(anEdge));
196	}
197	}
198	}
199
200	if (mySubst.IsCopied(myS1)) {
201	myShape = mySubst.Copy(myS1).First();
202	myShape.Orientation(myS1.Orientation());
203	Done();
204	}
205	}
206	}