0026937: Eliminate NO_CXX_EXCEPTION macro support
[occt.git] / src / TopOpeBRepDS / TopOpeBRepDS_Edge3dInterferenceTool.cxx
1 // Created on: 1997-12-24
2 // Created by: Prestataire Xuan PHAM PHU
3 // Copyright (c) 1997-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
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.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17
18 #include <BRep_Tool.hxx>
19 #include <BRepAdaptor_Curve.hxx>
20 #include <BRepAdaptor_Curve2d.hxx>
21 #include <gp_Vec.hxx>
22 #include <Precision.hxx>
23 #include <Standard_ProgramError.hxx>
24 #include <TopoDS.hxx>
25 #include <TopoDS_Shape.hxx>
26 #include <TopOpeBRepDS_CurvePointInterference.hxx>
27 #include <TopOpeBRepDS_define.hxx>
28 #include <TopOpeBRepDS_Edge3dInterferenceTool.hxx>
29 #include <TopOpeBRepDS_EdgeVertexInterference.hxx>
30 #include <TopOpeBRepDS_Interference.hxx>
31 #include <TopOpeBRepTool_EXPORT.hxx>
32 #include <TopOpeBRepTool_TOOL.hxx>
33
34 #define M_FORWARD(st) (st == TopAbs_FORWARD)
35 #define M_REVERSED(st) (st == TopAbs_REVERSED)
36
37 static void FUN_Raise()
38 {
39   throw Standard_ProgramError("Edge3dInterferenceTool");
40 }
41
42 #define POINT        (0)
43 #define VERTEXonref  (1)
44 #define VERTEXonOO   (2)
45 #define VERTEXonOref (3)
46
47 // myIsVertex :
48 // ------------
49 // POINT :<Eref> interfers with <E> at a point
50 // <Eref> interfers with <E> at a vertex V,
51 //   VERTEXonref  : V is on shape of <Eref>
52 //   VERTEXonOO   : V is on shape of <E>
53 //   VERTEXonOref : V is on 2 shapes.
54
55 // myVonOO : only for VERTEXonOO || VERTEXonOref
56 // --------
57
58 // myP3d : only for POINT || VERTEXonref
59 // -------
60
61 static Standard_Boolean FUN_hasparam(const Handle(TopOpeBRepDS_Interference)& I, Standard_Real& paronE)
62
63   // prequesitory : shapes <SIX> -> edge <E>
64   // ? <paronE> = parameter of <G> on <E>
65   TopOpeBRepDS_Kind GT = I->GeometryType();  
66   Standard_Boolean point =  (GT == TopOpeBRepDS_POINT);
67   Standard_Boolean vertex = (GT == TopOpeBRepDS_VERTEX);
68   if (point) {
69     Handle(TopOpeBRepDS_CurvePointInterference) CPI = 
70       Handle(TopOpeBRepDS_CurvePointInterference)::DownCast(I);
71     if (CPI.IsNull()) return Standard_False;
72     paronE = CPI->Parameter();
73   }
74   if (vertex) {
75     Handle(TopOpeBRepDS_EdgeVertexInterference) EVI = 
76       Handle(TopOpeBRepDS_EdgeVertexInterference)::DownCast(I);    
77     if (EVI.IsNull()) return Standard_False;
78     paronE = EVI->Parameter(); 
79   }
80   return Standard_True;
81 }
82
83 static Standard_Boolean FUN_paronOOE
84 (const TopoDS_Edge& OOE,const Standard_Integer IsVertex, const TopoDS_Shape& VonOO, const gp_Pnt& P3d, 
85  Standard_Real& paronOOE)
86 {  
87   Standard_Boolean ok = Standard_False;
88   Standard_Boolean hasVonOO = (IsVertex > 1);
89   if (hasVonOO) ok = FUN_tool_parVonE(TopoDS::Vertex(VonOO),OOE,paronOOE);
90   else {
91     Standard_Real dist; ok = FUN_tool_projPonE(P3d,OOE,paronOOE,dist);
92     Standard_Real tol1 = BRep_Tool::Tolerance(OOE);
93     Standard_Real tol2 = tol1*1.e3;
94     Standard_Real tol = tol2; if (tol > 1.e-2) tol = 1.e-2;
95     if (ok) ok = (dist <= tol);
96   }
97   return ok;
98 }
99
100 static Standard_Boolean FUN_keepIonF
101 (const gp_Vec& tgref, const Standard_Real& parE, const TopoDS_Edge& E, const TopoDS_Face& F,
102  const Standard_Real& tola)
103 // returns true if an interference I=(TonF,G=point/vertex,S=<E>)
104 // is to add to the Edge3dInterferenceTool resolving 3d complex transitions
105 // on edge E
106 {  
107   gp_Vec tmp; Standard_Boolean ok = TopOpeBRepTool_TOOL::TggeomE(parE,E,tmp);
108   if (!ok) return Standard_False;    
109   gp_Dir tgE = gp_Dir(tmp);
110   Standard_Real prod = Abs(tgref.Dot(tgE));
111   if (Abs(1-prod) < tola) return Standard_False; // <Eref> & <E> are tangent edges 
112   gp_Vec dd; ok = FUN_tool_nggeomF(parE,E,F,dd); gp_Dir ngF(dd);
113   if (!ok) return Standard_False;  
114   prod = Abs((tgref^tgE).Dot(ngF));
115   if (Abs(1-prod) < tola) return Standard_False;
116   return Standard_True;
117 }
118
119 // ----------------------------------------------------------------------
120 //             EDGE/FACE interferences reducing :
121 //
122 // ----------------------------------------------------------------------
123
124 //=======================================================================
125 //function : TopOpeBRepDS_Edge3dInterferenceTool
126 //purpose  : 
127 //=======================================================================
128
129 TopOpeBRepDS_Edge3dInterferenceTool::TopOpeBRepDS_Edge3dInterferenceTool()
130 :  myFaceOriented(0),myrefdef(Standard_False)
131
132 }
133
134
135 //=======================================================================
136 //function : InitPointVertex
137 //purpose  : Initializes reference data for edge/face complex transition           
138 //=======================================================================
139 // I = (TonF, G=POINT/VERTEX, S=<E>) interference on <Eref>
140 // G has parameter <paronEref> on <Eref>, <paronE> on <E>
141
142 void TopOpeBRepDS_Edge3dInterferenceTool::InitPointVertex
143 (const Standard_Integer IsVertex, const TopoDS_Shape& VonOO)
144 {
145   myIsVertex = IsVertex;
146   if (IsVertex > 1) myVonOO = VonOO;
147 }
148
149 //=======================================================================
150 //function : Init
151 //purpose  : Initializes reference data for edge/face complex transition           
152 //=======================================================================
153 // I = (T on <F>, G=POINT/VERTEX, S=<E>) interference on <Eref>
154 // G has parameter <paronEref> on <Eref>, <paronE> on <E>
155 // -- <E> is edge of <F> --
156
157 void TopOpeBRepDS_Edge3dInterferenceTool::Init
158 (const TopoDS_Shape& Eref, 
159  const TopoDS_Shape& E, const TopoDS_Shape& F, 
160  const Handle(TopOpeBRepDS_Interference)& I)
161 {  
162   const TopoDS_Edge& EEref = TopoDS::Edge(Eref);
163   const TopoDS_Edge& EE = TopoDS::Edge(E);
164   const TopoDS_Face& FF = TopoDS::Face(F);
165   myrefdef = Standard_False;
166
167   myTole = Precision::Angular(); // NYI
168
169   Standard_Real pref=0.0; Standard_Boolean ok = ::FUN_hasparam(I, pref);
170   if (!ok) {FUN_Raise(); return;}  
171   // <myP3d> :
172   {
173     BRepAdaptor_Curve BC(EEref);
174     myP3d = BC.Value(pref);
175   }
176   gp_Vec tmp; ok = TopOpeBRepTool_TOOL::TggeomE(pref,EEref,tmp); 
177   if (!ok) {FUN_Raise(); return;}  
178   gp_Dir tgref(tmp);
179   
180   Standard_Real pOO; ok = ::FUN_paronOOE(EE,myIsVertex,myVonOO,myP3d,pOO);
181   if (!ok) {FUN_Raise(); return;}
182   ok = TopOpeBRepTool_TOOL::TggeomE(pOO,EE,tmp); 
183   if (!ok) {FUN_Raise(); return;}  
184   gp_Dir tgOO(tmp);
185
186   Standard_Real dot = tgref.Dot(tgOO);
187   dot = 1 - Abs(dot);
188   Standard_Real tola = Precision::Confusion();
189   Standard_Boolean Esdm = (Abs(dot) < tola);
190   if (Esdm) return;
191   // NYI : il faut rejeter les interf I = (T,G,S=E) / E sdm with Eref
192
193   // <myrefdef> :
194   ok = ::FUN_keepIonF(tgref,pOO,EE,FF,myTole);
195   if (!ok) {
196     // <Eref> is tangent to <F>,
197     // If the transition is FORWARD or REVERSED, it describes a 2d 
198     // transition (while crossing <E> on <F>), we do not keep it.
199     const TopAbs_Orientation& O = I->Transition().Orientation(TopAbs_IN);
200     Standard_Boolean is2d = M_FORWARD(O) || M_REVERSED(O);
201     if (is2d) return;
202   }
203   myrefdef = Standard_True;
204
205   // <myFaceOriented> :
206   myFaceOriented = I->Transition().Index();
207
208   // <myTgtref>
209   myTgtref = tgref;
210
211   gp_Dir Norm = tgOO^tgref;  
212   myTool.Reset(tgOO, Norm);
213 }
214
215
216 //=======================================================================
217 //function : Add
218 //purpose  : 
219 //=======================================================================
220 // I = (T on <F>, G=POINT/VERTEX, S=<E>) interference on <Eref>
221 void TopOpeBRepDS_Edge3dInterferenceTool::Add
222 (const TopoDS_Shape& Eref,const TopoDS_Shape& E, const TopoDS_Shape& F,const Handle(TopOpeBRepDS_Interference)& I)
223 {
224   if (!myrefdef) {
225     Init(Eref,E,F,I);
226 //    return;
227   }
228
229   if (!myrefdef) return;
230
231   const TopoDS_Edge& EE = TopoDS::Edge(E);
232   const TopoDS_Face& FF = TopoDS::Face(F);
233
234   Standard_Real pOO; Standard_Boolean ok = ::FUN_paronOOE(EE,myIsVertex,myVonOO,myP3d,pOO);
235   if (!ok) return;
236   gp_Pnt2d uv;
237   {
238     BRepAdaptor_Curve2d BC2d(EE,FF);
239     uv = BC2d.Value(pOO);
240   }
241
242   ok = ::FUN_keepIonF(myTgtref,pOO,EE,FF,myTole);
243   if (!ok) {
244     const TopAbs_Orientation& O = I->Transition().Orientation(TopAbs_IN);
245     Standard_Boolean is2d = M_FORWARD(O) || M_REVERSED(O);
246     if (is2d) return;
247   }
248
249   TopAbs_Orientation oriloc = I->Transition().Orientation(TopAbs_IN);
250   TopAbs_Orientation oritan; ok = FUN_tool_orientEinFFORWARD(EE, FF, oritan); // xpu : 30/12/97
251   if (!ok) return;
252
253   gp_Dir Norm(FUN_tool_nggeomF(uv,FF)); 
254   myTool.Compare(myTole, Norm, oriloc, oritan);
255 }
256
257
258 //=======================================================================
259 //function : Transition
260 //purpose  : 
261 //=======================================================================
262
263 void TopOpeBRepDS_Edge3dInterferenceTool::Transition(const Handle(TopOpeBRepDS_Interference)& I) const 
264 {
265   TopOpeBRepDS_Transition& T = I->ChangeTransition();
266   I->Support(myFaceOriented);
267
268   TopAbs_State stb = myTool.StateBefore();
269   TopAbs_State sta = myTool.StateAfter();
270   T.Set(stb,sta);
271 }