Merge OCC22357 and replacing french comments by english one
[occt.git] / src / Select3D / Select3D_SensitiveFace.cxx
CommitLineData
7fd59977 1// Copyright: Matra-Datavision 1995
2// File: Select3D_SensitiveFace.cxx
3// Created: Mon Mar 27 10:15:15 1995
4// Author: Robert COUBLANC
5// <rob>
6//Modif on jun-24-97 : introduction de CSLib_Class2d de LBR
7// pour teste si on est dedans ou dehors...
8//Modif on jul-21-97 : changement en harray1 pour eventuelles connexions ...
9
10#include <Select3D_SensitiveFace.ixx>
4bf18dff 11#include <Select3D_Projector.hxx>
7fd59977 12#include <SelectBasics_BasicTool.hxx>
13#include <gp_Pnt2d.hxx>
14#include <gp_Pnt.hxx>
15#include <Precision.hxx>
16#include <ElCLib.hxx>
17
18#include <CSLib_Class2d.hxx>
19
20
21#define AutoInterMask 0x01
22#define AutoComputeMask 0x02
23// Standard_True if the flag is one
24#define AutoInterFlag(aflag) ( aflag & AutoInterMask )
25#define AutoComputeFlag(aflag) ( aflag & AutoComputeMask )
26// set the flag to one
27#define SetAutoInterFlag(aflag) ( aflag = aflag & AutoInterMask)
28#define SetAutoComputeFlag(aflag) ( aflag = aflag & AutoComputeMask)
29// Initialize flags
30#define AutoInitFlags(aflag) (aflag = 0)
31
32//==================================================
81bba717 33// Function: Hide this constructor to the next version...
34// Purpose : simply avoid interfering with the version update
7fd59977 35//==================================================
36
37Select3D_SensitiveFace::
38Select3D_SensitiveFace(const Handle(SelectBasics_EntityOwner)& OwnerId,
39 const TColgp_Array1OfPnt& ThePoints,
40 const Select3D_TypeOfSensitivity aType):
41Select3D_SensitivePoly(OwnerId, ThePoints),
42mytype (aType),
43myDetectedIndex(-1)
44{
45 AutoInitFlags(myautointer);
46}
47
48//==================================================
49// Function:
50// Purpose :
51//==================================================
52
53Select3D_SensitiveFace::
54Select3D_SensitiveFace(const Handle(SelectBasics_EntityOwner)& OwnerId,
55 const Handle(TColgp_HArray1OfPnt)& ThePoints,
56 const Select3D_TypeOfSensitivity aType):
57Select3D_SensitivePoly(OwnerId, ThePoints),
58mytype (aType),
59myDetectedIndex(-1)
60{
61 AutoInitFlags(myautointer);
62}
63
64//==================================================
65// Function:
66// Purpose :
67//==================================================
68
69Standard_Boolean Select3D_SensitiveFace::
70Matches(const Standard_Real X,
71 const Standard_Real Y,
72 const Standard_Real aTol,
73 Standard_Real& DMin)
74{
75#ifndef DEB
76 Standard_Real DMin2 = 0.;
77#else
78 Standard_Real DMin2;
79#endif
80 Standard_Real Xmin,Ymin,Xmax,Ymax;
81 if(!Bnd_Box2d(mybox2d).IsVoid()){
82 Bnd_Box2d(mybox2d).Get(Xmin,Ymin,Xmax,Ymax);
83 DMin2 = gp_XY(Xmax-Xmin,Ymax-Ymin).SquareModulus();
84 }
81bba717 85 // calculation of a criterion of minimum distance...
86 // from start Dmin = size of the bounding box 2D,
87 // then min. distance of the polyhedron or cdg...
7fd59977 88
89 Standard_Real aTol2 = aTol*aTol;
90 gp_XY CDG(0.,0.);
91// for(Standard_Integer I=1;I<=Nbp-1;I++){
92 Standard_Integer I;
93 for(I=1;I<mynbpoints-1;I++){
94 CDG+=((Select3D_Pnt2d*)mypolyg2d)[I-1];
95 }
96
97 if(mynbpoints>1){
98 CDG/= (mynbpoints-1);
99 }
100 DMin2=Min(DMin2,gp_XY(CDG.X()-X,CDG.Y()-Y).SquareModulus());
101 DMin = Sqrt(DMin2);
102
103
104 Standard_Boolean isplane2d(Standard_True);
105
106 for( I=1;I<mynbpoints-1;I++){
107 gp_XY V1(((Select3D_Pnt2d*)mypolyg2d)[I]),V(X,Y);
108 V1-=((Select3D_Pnt2d*)mypolyg2d)[I-1];
109 V-=((Select3D_Pnt2d*)mypolyg2d)[I-1];
110 Standard_Real Vector = V1^V;
111 Standard_Real V1V1 = V1.SquareModulus();
112 DMin2 =
113 (V1V1 <=aTol2) ?
81bba717 114 Min(DMin2,V.SquareModulus()): // if the segment is too small...
7fd59977 115 Min(DMin2,Vector*Vector/V1V1);
116 //cdg ...
117 gp_XY PlaneTest(CDG);PlaneTest-=((Select3D_Pnt2d*)mypolyg2d)[I-1];
118 Standard_Real valtst = PlaneTest^V1;
119 if(isplane2d && Abs(valtst)>aTol) isplane2d=Standard_False;
120 }
4bf18dff 121 if (isplane2d)
122 {
123 return Select3D_SensitiveEntity::Matches(X,Y,aTol,DMin);
7fd59977 124 }
125 //detection d'une auto - intersection dans le polygon 2D; si oui on sort
126// if (!AutoComputeFlag(myautointer)) {
127// if(mypolyg2d.Length()>4) {
128// if (SelectBasics_BasicTool::AutoInter(mypolyg2d)) {
129// SetAutoInterFlag(myautointer);
130// }
131// }
132// SetAutoComputeFlag(myautointer);
133// }
134// if (AutoInterFlag(myautointer)) return Standard_True;
135// //
136
81bba717 137 //otherwise it is checked if the point is in the face...
7fd59977 138 TColgp_Array1OfPnt2d aArrayOf2dPnt(1, mynbpoints);
139 Points2D(aArrayOf2dPnt);
140 CSLib_Class2d TheInOutTool(aArrayOf2dPnt,aTol,aTol,Xmin,Ymin,Xmax,Ymax);
141 Standard_Integer TheStat = TheInOutTool.SiDans(gp_Pnt2d(X,Y));
142
143 Standard_Boolean res(Standard_False);
144 switch(TheStat){
145 case 0:
146 res = Standard_True;
147 case 1:
148 {
149 if(mytype!=Select3D_TOS_BOUNDARY)
150 res = Standard_True;
151 }
152 }
4bf18dff 153 if (res)
154 {
155 return Select3D_SensitiveEntity::Matches(X,Y,aTol,DMin);
156 }
157 return Standard_False;
7fd59977 158}
159
160//=======================================================================
161//function : Matches
162//purpose :
163//=======================================================================
164
165Standard_Boolean Select3D_SensitiveFace::
166Matches (const Standard_Real XMin,
167 const Standard_Real YMin,
168 const Standard_Real XMax,
169 const Standard_Real YMax,
170 const Standard_Real aTol)
171{
172 Bnd_Box2d BoundBox;
173 BoundBox.Update(XMin-aTol,YMin-aTol,XMax+aTol,YMax+aTol);
174
175 for(Standard_Integer j=1;j<=mynbpoints-1;j++){
176 if(BoundBox.IsOut(((Select3D_Pnt2d*)mypolyg2d)[j-1])) return Standard_False;
177 }
178 return Standard_True;
179}
180
181//=======================================================================
182//function : Matches
183//purpose :
184//=======================================================================
185
186Standard_Boolean Select3D_SensitiveFace::
187Matches (const TColgp_Array1OfPnt2d& aPoly,
188 const Bnd_Box2d& aBox,
189 const Standard_Real aTol)
190{
191 Standard_Real Umin,Vmin,Umax,Vmax;
192 aBox.Get(Umin,Vmin,Umax,Vmax);
193 Standard_Real Tolu,Tolv;
194 Tolu = 1e-7;
195 Tolv = 1e-7;
196 CSLib_Class2d aClassifier2d(aPoly,aTol,aTol,Umin,Vmin,Umax,Vmax);
197
198 for(Standard_Integer j=1;j<=mynbpoints;j++){
199 Standard_Integer RES = aClassifier2d.SiDans(((Select3D_Pnt2d*)mypolyg2d)[j-1]);
200 if(RES!=1) return Standard_False;
201 }
202 return Standard_True;
203}
204
205
206//=======================================================================
207//function : Dump
208//purpose :
209//=======================================================================
210void Select3D_SensitiveFace::Dump(Standard_OStream& S,const Standard_Boolean FullDump) const
211{
212 S<<"\tSensitiveFace 3D :"<<endl;;
213 if(HasLocation())
214 S<<"\t\tExisting Location"<<endl;
215
216 if(mytype==Select3D_TOS_BOUNDARY)
217 S<<"\t\tSelection Of Bounding Polyline Only"<<endl;
218
219 if(FullDump){
220 S<<"\t\tNumber Of Points :"<<mynbpoints<<endl;
221
222// S<<"\t\t\tOwner:"<<myOwnerId<<endl;
223 Select3D_SensitiveEntity::DumpBox(S,mybox2d);
224 }
225}
226
227//=======================================================================
228//function : ComputeDepth
229//purpose :
230//=======================================================================
231Standard_Real Select3D_SensitiveFace::ComputeDepth(const gp_Lin& EyeLine) const
232{
4bf18dff 233 Standard_Real aDepth = Precision::Infinite();
234 Standard_Real aDepthMin = !mylastprj.IsNull() ? mylastprj->DepthMin() : -Precision::Infinite();
235 Standard_Real aDepthMax = !mylastprj.IsNull() ? mylastprj->DepthMax() : Precision::Infinite();
236 Standard_Real aDepthTest;
237 for (Standard_Integer i = 0; i < mynbpoints - 1; i++)
238 {
239 aDepthTest = ElCLib::Parameter (EyeLine, ((Select3D_Pnt* )mypolyg3d)[i]);
240 if (aDepthTest < aDepth && (aDepthTest > aDepthMin) && (aDepthTest < aDepthMax))
241 {
242 aDepth = aDepthTest;
243 }
244 }
245 return aDepth;
7fd59977 246}