0026936: Drawbacks of inlining in new type system in OCCT 7.0 -- automatic
[occt.git] / src / AIS / AIS_PerpendicularRelation.cxx
CommitLineData
b311480e 1// Created on: 1996-12-05
2// Created by: Jean-Pierre COMBE/Odile Olivier
3// Copyright (c) 1996-1999 Matra Datavision
973c2be1 4// Copyright (c) 1999-2014 OPEN CASCADE SAS
b311480e 5//
973c2be1 6// This file is part of Open CASCADE Technology software library.
b311480e 7//
d5f74e42 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
973c2be1 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.
b311480e 13//
973c2be1 14// Alternatively, this file may be used under the terms of Open CASCADE
15// commercial license or contractual agreement.
7fd59977 16
7fd59977 17
18#include <AIS.hxx>
42cf5bc1 19#include <AIS_PerpendicularRelation.hxx>
7fd59977 20#include <BRep_Tool.hxx>
42cf5bc1 21#include <BRepAdaptor_Surface.hxx>
22#include <BRepBuilderAPI_MakeFace.hxx>
7fd59977 23#include <DsgPrs_PerpenPresentation.hxx>
7fd59977 24#include <ElCLib.hxx>
42cf5bc1 25#include <gce_MakeDir.hxx>
7fd59977 26#include <Geom2d_Line.hxx>
42cf5bc1 27#include <Geom_Ellipse.hxx>
7fd59977 28#include <Geom_Line.hxx>
29#include <Geom_Plane.hxx>
42cf5bc1 30#include <Geom_Transformation.hxx>
31#include <GeomAPI.hxx>
32#include <gp_Pln.hxx>
33#include <gp_Pnt.hxx>
34#include <gp_Pnt2d.hxx>
35#include <gp_Trsf.hxx>
36#include <gp_Vec.hxx>
7fd59977 37#include <IntAna2d_AnaIntersection.hxx>
38#include <IntAna2d_IntPoint.hxx>
7fd59977 39#include <Precision.hxx>
42cf5bc1 40#include <Prs3d_Presentation.hxx>
41#include <Prs3d_Projector.hxx>
7fd59977 42#include <Select3D_SensitiveSegment.hxx>
43#include <SelectMgr_EntityOwner.hxx>
f751596e 44#include <SelectMgr_Selection.hxx>
42cf5bc1 45#include <Standard_NotImplemented.hxx>
46#include <Standard_Type.hxx>
7fd59977 47#include <TopoDS.hxx>
48#include <TopoDS_Edge.hxx>
49#include <TopoDS_Face.hxx>
42cf5bc1 50#include <TopoDS_Shape.hxx>
7fd59977 51#include <TopoDS_Vertex.hxx>
52
92efcf78 53IMPLEMENT_STANDARD_RTTIEXT(AIS_PerpendicularRelation,AIS_Relation)
54
7fd59977 55//=======================================================================
56//function : Constructor
57//purpose : TwoEdgesPerpendicular
58//=======================================================================
59AIS_PerpendicularRelation::AIS_PerpendicularRelation(const TopoDS_Shape& aFShape,
60 const TopoDS_Shape& aSShape,
61 const Handle(Geom_Plane)& aPlane)
62:AIS_Relation()
63{
64 myFShape = aFShape;
65 mySShape = aSShape;
66 myPlane = aPlane;
67}
68
69//=======================================================================
70//function : Constructor
71//purpose : TwoFacesPerpendicular
72//=======================================================================
73AIS_PerpendicularRelation::AIS_PerpendicularRelation(const TopoDS_Shape& aFShape,
74 const TopoDS_Shape& aSShape)
75:AIS_Relation()
76{
77 myFShape = aFShape;
78 mySShape = aSShape;
79}
80
81//=======================================================================
82//function : Compute
83//purpose :
84//=======================================================================
85void AIS_PerpendicularRelation::Compute(const Handle(PrsMgr_PresentationManager3d)&,
86 const Handle(Prs3d_Presentation)& aPresentation,
87 const Standard_Integer)
88{
89 aPresentation->Clear();
90
91 if (myFShape.ShapeType() == mySShape.ShapeType()) {
92 switch (myFShape.ShapeType()) {
93 case TopAbs_FACE :
94 {
95 // cas perpendiculaire entre deux faces
96 ComputeTwoFacesPerpendicular(aPresentation);
97 }
98 break;
99 case TopAbs_EDGE :
100 {
101 // cas perpendiculaire entre deux edges
102 ComputeTwoEdgesPerpendicular(aPresentation);
103 }
104 break;
105 default:
106 break;
107 }
108 }
109 // Cas pas traite - Edge/Face
110}
111
112//=======================================================================
113//function : Compute
114//purpose : to avoid warning
115//=======================================================================
116void AIS_PerpendicularRelation::Compute(const Handle(Prs3d_Projector)& aProjector,
117 const Handle(Prs3d_Presentation)& aPresentation)
118{
119// Standard_NotImplemented::Raise("AIS_PerpendicularRelation::Compute(const Handle(Prs3d_Projector)&,const Handle(Prs3d_Presentation)&)");
120 PrsMgr_PresentableObject::Compute( aProjector , aPresentation ) ;
121}
122
857ffd5e 123void AIS_PerpendicularRelation::Compute(const Handle(Prs3d_Projector)& aProjector, const Handle(Geom_Transformation)& aTransformation, const Handle(Prs3d_Presentation)& aPresentation)
7fd59977 124{
857ffd5e 125// Standard_NotImplemented::Raise("AIS_PerpendicularRelation::Compute(const Handle(Prs3d_Projector)&, const Handle(Geom_Transformation)&, const Handle(Prs3d_Presentation)&)");
7fd59977 126 PrsMgr_PresentableObject::Compute( aProjector , aTransformation , aPresentation ) ;
127}
128
129//=======================================================================
130//function : ComputeSelection
131//purpose :
132//=======================================================================
133void AIS_PerpendicularRelation::ComputeSelection(const Handle(SelectMgr_Selection)& aSelection,
134 const Standard_Integer)
135{
136 Handle(SelectMgr_EntityOwner) own = new SelectMgr_EntityOwner(this,7);
137 const gp_Pnt& pos = myPosition;
138 Handle(Select3D_SensitiveSegment) seg;
139 Standard_Boolean ok1(Standard_False),ok2(Standard_False);
140
141 if (!myFAttach.IsEqual(pos,Precision::Confusion())) {
142 seg = new Select3D_SensitiveSegment(own,
143 myFAttach,
144 pos);
145 aSelection->Add(seg);
146 ok1 = Standard_True;
147 }
148 if (!mySAttach.IsEqual(myPosition,Precision::Confusion())) {
149 seg = new Select3D_SensitiveSegment(own,
150 mySAttach,
151 pos);
152 aSelection->Add(seg);
153 ok2 = Standard_True;
154 }
155
156 if (ok1 && ok2) {
157 gp_Vec vec1(gce_MakeDir(pos,myFAttach));
158 gp_Vec vec2(gce_MakeDir(pos,mySAttach));
159 Standard_Real dist1(pos.Distance(myFAttach));
160 Standard_Real dist2(pos.Distance(mySAttach));
161 vec1 *= dist1;
162 vec1 *= .2;
163 vec2 *= dist2;
164 vec2 *= .2;
165
166 gp_Pnt pAx11 = pos.Translated(vec1);
167 gp_Pnt pAx22 = pos.Translated(vec2);
168 gp_Pnt p_symb = pAx22.Translated(vec1);
169 seg = new Select3D_SensitiveSegment(own,pAx11,p_symb);
170 aSelection->Add(seg);
171 seg = new Select3D_SensitiveSegment(own,p_symb,pAx22);
172 aSelection->Add(seg);
173 }
174}
175
176//=======================================================================
177//function : ComputeTwoFacesPerpendicular
178//purpose :
179//=======================================================================
180void AIS_PerpendicularRelation::ComputeTwoFacesPerpendicular
181 (const Handle(Prs3d_Presentation)& /*aPresentation*/)
182{
183}
184
185//=======================================================================
186//function : ComputeTwoEdgesPerpendicular
187//purpose :
188//=======================================================================
189void AIS_PerpendicularRelation::ComputeTwoEdgesPerpendicular(const Handle(Prs3d_Presentation)& aPresentation)
190{
191 // 3d lines
192 Handle(Geom_Curve) geom1,geom2;
193 gp_Pnt pint3d,p1,p2,pAx1,pAx2,ptat11,ptat12,ptat21,ptat22;
194 Standard_Boolean isInfinite1,isInfinite2;
195 Handle(Geom_Curve) extCurv;
196 if ( !AIS::ComputeGeometry(TopoDS::Edge(myFShape),TopoDS::Edge(mySShape),
197 myExtShape,
198 geom1,geom2,
199 ptat11,ptat12,ptat21,ptat22,
200 extCurv,
201 isInfinite1,isInfinite2,
202 myPlane) ) return;
203
204 Standard_Boolean interOut1(Standard_False),interOut2(Standard_False);
205
206 Handle(Geom_Line) geom_lin1;
207 Handle(Geom_Line) geom_lin2;
208 if ( geom1->IsInstance(STANDARD_TYPE(Geom_Ellipse)) )
209 {
c5f3a425 210 Handle(Geom_Ellipse) geom_el (Handle(Geom_Ellipse)::DownCast (geom1));
7fd59977 211 // construct lines through focuses
212 gp_Ax1 elAx = geom_el->XAxis();
213 gp_Lin ll (elAx);
214 geom_lin1 = new Geom_Line(ll);
215 Standard_Real focex = geom_el->MajorRadius() - geom_el->Focal()/2.0;
216 gp_Vec transvec = gp_Vec(elAx.Direction())*focex;
217 ptat11 = geom_el->Focus1().Translated(transvec);
218 ptat12 = geom_el->Focus2().Translated(-transvec);
219 interOut1 = Standard_True;
220 }
221 else if ( geom1->IsInstance(STANDARD_TYPE(Geom_Line)) )
222 {
c5f3a425 223 geom_lin1 = Handle(Geom_Line)::DownCast (geom1);
7fd59977 224 }
225 else return;
226
227 if (geom2->IsInstance(STANDARD_TYPE(Geom_Ellipse)))
228 {
c5f3a425 229 Handle(Geom_Ellipse) geom_el (Handle(Geom_Ellipse)::DownCast (geom2));
7fd59977 230 // construct lines through focuses
231 gp_Ax1 elAx = geom_el->XAxis();
232 gp_Lin ll (elAx);
233 geom_lin2 = new Geom_Line(ll);
234 Standard_Real focex = geom_el->MajorRadius() - geom_el->Focal()/2.0;
235 gp_Vec transvec = gp_Vec(elAx.Direction())*focex;
236 ptat21 = geom_el->Focus1().Translated(transvec);
237 ptat22 = geom_el->Focus2().Translated(-transvec);
238 interOut2 = Standard_True;
239 }
240 else if ( geom2->IsInstance(STANDARD_TYPE(Geom_Line)) )
241 {
c5f3a425 242 geom_lin2 = Handle(Geom_Line)::DownCast (geom2);
7fd59977 243 }
244 else return;
245
246 // current face
247 BRepBuilderAPI_MakeFace makeface (myPlane->Pln());
248 TopoDS_Face face (makeface.Face());
249 BRepAdaptor_Surface adp (makeface.Face());
250
251 // 2d lines => projection of 3d on current plane
252 Handle(Geom2d_Curve) aGeom2dCurve = GeomAPI::To2d(geom_lin1,myPlane->Pln());
c5f3a425 253 Handle(Geom2d_Line) lin1_2d = Handle(Geom2d_Line)::DownCast (aGeom2dCurve) ;
7fd59977 254 aGeom2dCurve = GeomAPI::To2d(geom_lin2,myPlane->Pln());
c5f3a425 255 Handle(Geom2d_Line) lin2_2d = Handle(Geom2d_Line)::DownCast (aGeom2dCurve) ;
7fd59977 256 IntAna2d_AnaIntersection inter(lin1_2d->Lin2d(),lin2_2d->Lin2d());
257 if (!inter.IsDone()) return;
258 if (!inter.NbPoints()) return;
259
260 gp_Pnt2d pint(inter.Point(1).Value());
261 pint3d = adp.Value(pint.X(),pint.Y());
262
263 myPosition = pint3d;
264 // recherche points attache
265 Standard_Real par1,par2,curpar,pmin,pmax;//,dist,sign;
266 Standard_Real length(0.);
267
268 if ( isInfinite1 && isInfinite2 )
269 {
270 Standard_Real curpar1 = ElCLib::Parameter(geom_lin1->Lin(),pint3d);
271 Standard_Real curpar2 = ElCLib::Parameter(geom_lin2->Lin(),pint3d);
272 par1 = par2 = 50.;
273 p1 = p2 = pint3d;
274 myFAttach = ElCLib::Value(curpar1+par1,geom_lin1->Lin());
275 mySAttach = ElCLib::Value(curpar2+par2,geom_lin2->Lin());
276 }
277 else
278 {
279 Standard_Boolean lengthComputed (Standard_False);
280 if ( !isInfinite1 )
281 {
282 curpar = ElCLib::Parameter(geom_lin1->Lin(),pint3d);
283 par1 = ElCLib::Parameter(geom_lin1->Lin(),ptat11);
284 par2 = ElCLib::Parameter(geom_lin1->Lin(),ptat12);
285 pmin = Min(par1,par2);
286 pmax = Max(par1,par2);
287
288 if ( myPosition.SquareDistance(ptat11) > myPosition.SquareDistance(ptat12) )
289 p1 = ptat11;
290 else
291 p1 = ptat12;
292 if ( (curpar < pmin) || (curpar > pmax) )
293 {
294 interOut1 = Standard_True;
295 }
296 if ( !isInfinite2 ) length = 2.*Min(ptat11.Distance(ptat12),ptat21.Distance(ptat22))/5.;
297 else length = 2.*ptat11.Distance(ptat12)/5.;
298 lengthComputed = Standard_True;
299 gp_Vec vec1 (gce_MakeDir(myPosition,p1));
300 vec1.Multiply(length);
301 pAx1 = myPosition.Translated(vec1);
302 myFAttach = pAx1;
303 }
304 if ( !isInfinite2 )
305 {
306 curpar = ElCLib::Parameter(geom_lin2->Lin(),pint3d);
307 par1 = ElCLib::Parameter(geom_lin2->Lin(),ptat21);
308 par2 = ElCLib::Parameter(geom_lin2->Lin(),ptat22);
309 pmin = Min(par1,par2);
310 pmax = Max(par1,par2);
311
312 if ( myPosition.SquareDistance(ptat21) > myPosition.SquareDistance(ptat22) ) p2 = ptat21;
313 else p2 = ptat22;
314 if ( (curpar < pmin) || (curpar > pmax) )
315 {
316 interOut2 = Standard_True;
317 }
318 gp_Vec vec2 (gce_MakeDir(myPosition,p2));
319 if ( !lengthComputed )
320 {
321 if ( !isInfinite1 ) length = 2.*Min(ptat11.Distance(ptat12),ptat21.Distance(ptat22))/5.;
322 else length = 2.*ptat21.Distance(ptat22)/5.;
323 }
324 vec2.Multiply(length);
325 pAx2 = myPosition.Translated(vec2);
326 mySAttach = pAx2;
327 }
328 if ( isInfinite1 )
329 {
330 p1 = myPosition;
331 gp_Vec vec1(geom_lin1->Lin().Direction());
332 vec1.Multiply(length);
333 myFAttach = myPosition.Translated(vec1);
334 }
335 if ( isInfinite2 )
336 {
337 p2 = myPosition;
338 gp_Vec vec2(geom_lin2->Lin().Direction());
339 vec2.Multiply(length);
340 mySAttach = myPosition.Translated(vec2);
341 }
342 }
343 DsgPrs_PerpenPresentation::Add(aPresentation,myDrawer,
344 myFAttach,mySAttach,
345 p1,p2,
346 myPosition,
347 interOut1,interOut2);
348
349 if ( (myExtShape != 0) && !extCurv.IsNull()) {
350 gp_Pnt pf,pl;
351 if ( myExtShape == 1 ) {
352 if (!isInfinite1) {
353 pf = ptat11;
354 pl = ptat12;
355 }
356 aPresentation->SetInfiniteState(isInfinite1);
357 ComputeProjEdgePresentation(aPresentation,TopoDS::Edge(myFShape),geom_lin1,pf,pl);
358 }
359 else {
360 if (!isInfinite2) {
361 pf = ptat21;
362 pl = ptat22;
363 }
364 aPresentation->SetInfiniteState(isInfinite2);
365 ComputeProjEdgePresentation(aPresentation,TopoDS::Edge(mySShape),geom_lin2,pf,pl);
366 }
367 }
368}