Warnings on vc14 were eliminated
[occt.git] / src / IntImpParGen / IntImpParGen_Tool.cxx
1 // Created on: 1992-06-10
2 // Created by: Laurent BUCHARD
3 // Copyright (c) 1992-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 #include <IntImpParGen_Tool.hxx>
18 #include <gp.hxx>
19
20
21 #define TOLERANCE_ANGULAIRE 0.00000001
22
23 //----------------------------------------------------------------------
24 Standard_Real NormalizeOnDomain(Standard_Real& Param,const IntRes2d_Domain& TheDomain) {
25   Standard_Real modParam = Param;
26   if(TheDomain.IsClosed()) {
27     Standard_Real Periode,t;
28     TheDomain.EquivalentParameters(t,Periode);
29     Periode-=t;
30     if(TheDomain.HasFirstPoint()) {
31       while(modParam<TheDomain.FirstParameter()) {
32         modParam+=Periode;
33       }
34     }
35     if(TheDomain.HasLastPoint()) {
36       if(modParam>TheDomain.LastParameter()) {
37         modParam-=Periode;
38       }
39     }
40   }
41   return(modParam);
42 }
43 //----------------------------------------------------------------------
44 void Determine_Position(IntRes2d_Position& Pos1,
45                         const IntRes2d_Domain& TheDomain,
46                         const gp_Pnt2d& Pnt1,
47                         const Standard_Real Param1) {
48   
49   Pos1=IntRes2d_Middle;
50
51   if(TheDomain.HasFirstPoint()) { 
52     if(Pnt1.Distance(TheDomain.FirstPoint()) 
53        <= TheDomain.FirstTolerance()) {
54       Pos1=IntRes2d_Head;       
55     }
56   }
57    
58   if(TheDomain.HasLastPoint()) {
59     if(Pnt1.Distance(TheDomain.LastPoint()) 
60        <= TheDomain.LastTolerance()) {
61       if(Pos1==IntRes2d_Head) {
62         if(Abs(Param1-TheDomain.LastParameter())
63            < Abs(Param1-TheDomain.FirstParameter()))
64           Pos1=IntRes2d_End;    
65       }
66       else {
67         Pos1=IntRes2d_End; 
68       }
69     } 
70   }
71 }   
72 //----------------------------------------------------------------------
73 void Determine_Transition(const IntRes2d_Position    Pos1,
74                           gp_Vec2d&                  Tan1,
75                           const gp_Vec2d&            Norm1,
76                           IntRes2d_Transition&       T1,
77                           const IntRes2d_Position    Pos2,
78                           gp_Vec2d&                  Tan2,
79                           const gp_Vec2d&            Norm2,
80                           IntRes2d_Transition&       T2,
81 //                        const Standard_Real        Tolerance_Angulaire) {
82                           const Standard_Real        ) {
83   
84   Standard_Boolean courbure1=Standard_True;
85   Standard_Boolean courbure2=Standard_True;
86   Standard_Boolean decide=Standard_True;
87   if (Tan1.Magnitude()<=gp::Resolution()) {
88     Tan1=Norm1;
89     courbure1=Standard_False;
90     if (Tan1.Magnitude()<=gp::Resolution()) {    // transition undecided
91       decide=Standard_False;
92     }
93   }
94   
95   if (Tan2.Magnitude()<=gp::Resolution()) {
96     Tan2=Norm2;
97     courbure2=Standard_False;
98     if (Tan2.Magnitude()<=gp::Resolution()) {    // transition undecided
99       decide=Standard_False;
100     }
101   }
102   
103   if (!decide) {
104     T1.SetValue(Pos1);
105     T2.SetValue(Pos2);
106   }
107   else {
108     Standard_Real sgn=Tan1.Crossed(Tan2);
109     Standard_Real norm=Tan1.Magnitude()*Tan2.Magnitude();
110
111     if (Abs(sgn)<=TOLERANCE_ANGULAIRE*norm) {   // Transition TOUCH #########
112       Standard_Boolean opos=(Tan1.Dot(Tan2))<0;
113       if (!(courbure1||courbure2)) {
114         T1.SetValue(Standard_True,Pos1,IntRes2d_Unknown,opos);
115         T2.SetValue(Standard_True,Pos2,IntRes2d_Unknown,opos);
116       }
117       else {
118         gp_Vec2d Norm;
119         Tan1.Normalized();
120         Norm.SetCoord(-Tan1.Y(),Tan1.X());
121         Standard_Real Val1,Val2;
122         if (!courbure1) {
123           Val1=0.0;
124         }
125         else {
126           Val1=Norm.Dot(Norm1);
127         }
128         if (!courbure2) {
129           Val2=0.0;
130         }
131         else {
132           Val2=Norm.Dot(Norm2);
133         }
134         
135         if (Abs(Val1-Val2) <= gp::Resolution()) {
136           T1.SetValue(Standard_True,Pos1,IntRes2d_Unknown,opos);
137           T2.SetValue(Standard_True,Pos2,IntRes2d_Unknown,opos);
138         }
139         else if (Val2 > Val1) {
140           T2.SetValue(Standard_True,Pos2,IntRes2d_Inside,opos);
141           if (opos) {
142             T1.SetValue(Standard_True,Pos1,IntRes2d_Inside,opos);
143           }
144           else {
145             T1.SetValue(Standard_True,Pos1,IntRes2d_Outside,opos);
146           }
147         }
148         else {         // Val1 > Val2
149           T2.SetValue(Standard_True,Pos2,IntRes2d_Outside,opos);
150           if (opos) {
151             T1.SetValue(Standard_True,Pos1,IntRes2d_Outside,opos);
152           }
153           else {
154             T1.SetValue(Standard_True,Pos1,IntRes2d_Inside,opos);
155           }
156         }
157       }
158     }
159     else if (sgn<0) {
160       T1.SetValue(Standard_False,Pos1,IntRes2d_In);
161       T2.SetValue(Standard_False,Pos2,IntRes2d_Out);
162     }
163     else {     // sgn>0
164       T1.SetValue(Standard_False,Pos1,IntRes2d_Out);
165       T2.SetValue(Standard_False,Pos2,IntRes2d_In);
166     }
167   }
168 }
169