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
6 // This file is part of Open CASCADE Technology software library.
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.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
19 #include <gp_Pnt2d.hxx>
20 #include <gp_Vec2d.hxx>
21 #include <IntImpParGen.hxx>
22 #include <IntImpParGen_Tool.hxx>
23 #include <IntRes2d_Domain.hxx>
24 #include <IntRes2d_Position.hxx>
25 #include <IntRes2d_Transition.hxx>
27 #define TOLERANCE_ANGULAIRE 0.00000001
28 #define DERIVEE_PREMIERE_NULLE 0.000000000001
29 //----------------------------------------------------------------------
30 Standard_Real IntImpParGen::NormalizeOnDomain(Standard_Real& Param,
31 const IntRes2d_Domain& TheDomain) {
32 Standard_Real modParam = Param;
33 if(TheDomain.IsClosed()) {
34 Standard_Real Periode,t;
35 TheDomain.EquivalentParameters(t,Periode);
37 while( modParam<TheDomain.FirstParameter()
38 && modParam+Periode < TheDomain.LastParameter()) {
41 while( modParam>TheDomain.LastParameter()
42 && modParam-Periode > TheDomain.FirstParameter()) {
48 //----------------------------------------------------------------------
49 void IntImpParGen::DeterminePosition(IntRes2d_Position& Pos1,
50 const IntRes2d_Domain& TheDomain,
52 const Standard_Real Param1) {
56 if(TheDomain.HasFirstPoint()) {
57 if(Pnt1.Distance(TheDomain.FirstPoint())
58 <= TheDomain.FirstTolerance()) {
63 if(TheDomain.HasLastPoint()) {
64 if(Pnt1.Distance(TheDomain.LastPoint())
65 <= TheDomain.LastTolerance()) {
66 if(Pos1==IntRes2d_Head) {
67 if(Abs(Param1-TheDomain.LastParameter())
68 < Abs(Param1-TheDomain.FirstParameter()))
77 //----------------------------------------------------------------------
78 void IntImpParGen::DetermineTransition(const IntRes2d_Position Pos1,
80 const gp_Vec2d& Norm1,
81 IntRes2d_Transition& T1,
82 const IntRes2d_Position Pos2,
84 const gp_Vec2d& Norm2,
85 IntRes2d_Transition& T2,
86 const Standard_Real ) {
88 Standard_Boolean courbure1=Standard_True;
89 Standard_Boolean courbure2=Standard_True;
90 Standard_Boolean decide=Standard_True;
96 if (Tan1.SquareMagnitude()<=DERIVEE_PREMIERE_NULLE) {
98 courbure1=Standard_False;
99 if (Tan1.SquareMagnitude()<=DERIVEE_PREMIERE_NULLE) { // transition undecided
100 decide=Standard_False;
104 if (Tan2.SquareMagnitude()<=DERIVEE_PREMIERE_NULLE) {
106 courbure2=Standard_False;
107 if (Tan2.SquareMagnitude()<=DERIVEE_PREMIERE_NULLE) { // transition undecided
108 decide=Standard_False;
117 Standard_Real sgn=Tan1.Crossed(Tan2);
118 Standard_Real norm=Tan1.Magnitude()*Tan2.Magnitude();
120 if (Abs(sgn)<=TOLERANCE_ANGULAIRE*norm) { // Transition TOUCH #########
121 Standard_Boolean opos=(Tan1.Dot(Tan2))<0;
122 if (!(courbure1||courbure2)) {
123 T1.SetValue(Standard_True,Pos1,IntRes2d_Unknown,opos);
124 T2.SetValue(Standard_True,Pos2,IntRes2d_Unknown,opos);
129 Norm.SetCoord(-Tan1.Y(),Tan1.X());
130 Standard_Real Val1,Val2;
135 Val1=Norm.Dot(Norm1);
141 Val2=Norm.Dot(Norm2);
144 if (Abs(Val1-Val2) <= TOLERANCE_ANGULAIRE) {
145 T1.SetValue(Standard_True,Pos1,IntRes2d_Unknown,opos);
146 T2.SetValue(Standard_True,Pos2,IntRes2d_Unknown,opos);
148 else if (Val2 > Val1) {
149 T2.SetValue(Standard_True,Pos2,IntRes2d_Inside,opos);
151 T1.SetValue(Standard_True,Pos1,IntRes2d_Inside,opos);
154 T1.SetValue(Standard_True,Pos1,IntRes2d_Outside,opos);
157 else { // Val1 > Val2
158 T2.SetValue(Standard_True,Pos2,IntRes2d_Outside,opos);
160 T1.SetValue(Standard_True,Pos1,IntRes2d_Outside,opos);
163 T1.SetValue(Standard_True,Pos1,IntRes2d_Inside,opos);
169 T1.SetValue(Standard_False,Pos1,IntRes2d_In);
170 T2.SetValue(Standard_False,Pos2,IntRes2d_Out);
173 T1.SetValue(Standard_False,Pos1,IntRes2d_Out);
174 T2.SetValue(Standard_False,Pos2,IntRes2d_In);
179 //----------------------------------------------------------------------
180 Standard_Boolean IntImpParGen::DetermineTransition(const IntRes2d_Position Pos1,
182 IntRes2d_Transition& T1,
183 const IntRes2d_Position Pos2,
185 IntRes2d_Transition& T2,
186 const Standard_Real ) {
188 T1.SetPosition(Pos1);
189 T2.SetPosition(Pos2);
191 Standard_Real Tan1Magnitude = Tan1.Magnitude();
192 if (Tan1Magnitude<=DERIVEE_PREMIERE_NULLE) {
193 return(Standard_False);
196 Standard_Real Tan2Magnitude = Tan2.Magnitude();
197 if (Tan2Magnitude<=DERIVEE_PREMIERE_NULLE) {
198 return(Standard_False);
201 Standard_Real sgn=Tan1.Crossed(Tan2);
202 Standard_Real norm=Tan1Magnitude*Tan2Magnitude;
204 if (Abs(sgn)<=TOLERANCE_ANGULAIRE*norm) { // Transition TOUCH #########
205 return(Standard_False);
208 T1.SetValue(Standard_False,Pos1,IntRes2d_In);
209 T2.SetValue(Standard_False,Pos2,IntRes2d_Out);
212 T1.SetValue(Standard_False,Pos1,IntRes2d_Out);
213 T2.SetValue(Standard_False,Pos2,IntRes2d_In);
215 return(Standard_True);