0024002: Overall code and build procedure refactoring -- automatic
[occt.git] / src / AIS / AIS_ParallelRelation.cxx
1 // Created on: 1996-12-05
2 // Created by: Jean-Pierre COMBE/Odile Olivier
3 // Copyright (c) 1996-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 <AIS.hxx>
19 #include <AIS_ParallelRelation.hxx>
20 #include <BRep_Tool.hxx>
21 #include <BRepAdaptor_Curve.hxx>
22 #include <BRepAdaptor_Surface.hxx>
23 #include <DsgPrs_LengthPresentation.hxx>
24 #include <ElCLib.hxx>
25 #include <gce_MakeLin.hxx>
26 #include <Geom_Ellipse.hxx>
27 #include <Geom_Line.hxx>
28 #include <Geom_Plane.hxx>
29 #include <Geom_Transformation.hxx>
30 #include <gp_Ax1.hxx>
31 #include <gp_Ax2.hxx>
32 #include <gp_Dir.hxx>
33 #include <gp_Lin.hxx>
34 #include <gp_Pln.hxx>
35 #include <gp_Pnt.hxx>
36 #include <Precision.hxx>
37 #include <Prs3d_ArrowAspect.hxx>
38 #include <Prs3d_DimensionAspect.hxx>
39 #include <Prs3d_Drawer.hxx>
40 #include <Prs3d_Presentation.hxx>
41 #include <Prs3d_Projector.hxx>
42 #include <Select3D_SensitiveBox.hxx>
43 #include <Select3D_SensitiveSegment.hxx>
44 #include <SelectMgr_EntityOwner.hxx>
45 #include <SelectMgr_Selection.hxx>
46 #include <Standard_DomainError.hxx>
47 #include <Standard_NotImplemented.hxx>
48 #include <Standard_Type.hxx>
49 #include <TCollection_AsciiString.hxx>
50 #include <TCollection_ExtendedString.hxx>
51 #include <TopExp_Explorer.hxx>
52 #include <TopoDS.hxx>
53 #include <TopoDS_Shape.hxx>
54
55 //=======================================================================
56 //function : Constructor
57 //purpose  : 
58 //=======================================================================
59 AIS_ParallelRelation::AIS_ParallelRelation(const TopoDS_Shape& aFShape, 
60                                            const TopoDS_Shape& aSShape, 
61                                            const Handle(Geom_Plane)& aPlane)
62 {
63   myFShape = aFShape;
64   mySShape = aSShape;
65   myPlane = aPlane;
66   myAutomaticPosition = Standard_True;
67   myArrowSize = 0.01;
68   mySymbolPrs = DsgPrs_AS_BOTHAR;
69 }
70
71 //=======================================================================
72 //function : Constructor
73 //purpose  : 
74 //=======================================================================
75 AIS_ParallelRelation::AIS_ParallelRelation(const TopoDS_Shape& aFShape, 
76                                            const TopoDS_Shape& aSShape, 
77                                            const Handle(Geom_Plane)& aPlane, 
78                                            const gp_Pnt& aPosition,
79                                            const DsgPrs_ArrowSide aSymbolPrs, 
80                                            const Standard_Real anArrowSize)
81 {
82   myFShape = aFShape;
83   mySShape = aSShape;
84   myPlane = aPlane;
85   myAutomaticPosition = Standard_False;
86   SetArrowSize( anArrowSize );
87   myPosition = aPosition;
88   mySymbolPrs = aSymbolPrs;
89 }
90
91 //=======================================================================
92 //function : Compute
93 //purpose  : 
94 //=======================================================================
95 void AIS_ParallelRelation::Compute(const Handle(PrsMgr_PresentationManager3d)&, 
96                                    const Handle(Prs3d_Presentation)& aPresentation, 
97                                    const Standard_Integer)
98 {
99   aPresentation->Clear();
100
101   switch (myFShape.ShapeType())
102     {
103     case TopAbs_FACE :
104       {
105         // cas longueur entre deux faces
106         ComputeTwoFacesParallel(aPresentation);
107       }
108       break;
109     case TopAbs_EDGE :
110       {
111         // cas longueur entre deux edges
112         ComputeTwoEdgesParallel(aPresentation);
113       }
114       break;
115     default:
116       break;
117     }
118 }
119
120 //=======================================================================
121 //function : Compute
122 //purpose  : to avoid warning
123 //=======================================================================
124 void AIS_ParallelRelation::Compute(const Handle(Prs3d_Projector)& aProjector,
125                                    const Handle(Prs3d_Presentation)& aPresentation)
126 {
127 // Standard_NotImplemented::Raise("AIS_ParallelRelation::Compute(const Handle(Prs3d_Projector)&,const Handle(Prs3d_Presentation)&)");
128   PrsMgr_PresentableObject::Compute( aProjector , aPresentation ) ;
129 }
130
131 void AIS_ParallelRelation::Compute(const Handle(Prs3d_Projector)& aProjector, const Handle(Geom_Transformation)& aTransformation, const Handle(Prs3d_Presentation)& aPresentation)
132 {
133 // Standard_NotImplemented::Raise("AIS_ParallelRelation::Compute(const Handle(Prs3d_Projector)&, const Handle(Geom_Transformation)&, const Handle(Prs3d_Presentation)&)");
134   PrsMgr_PresentableObject::Compute( aProjector , aTransformation , aPresentation ) ;
135 }
136
137 //=======================================================================
138 //function : ComputeSelection
139 //purpose  : 
140 //=======================================================================
141 void AIS_ParallelRelation::ComputeSelection(const Handle(SelectMgr_Selection)& aSelection, 
142                                             const Standard_Integer)
143 {
144   gp_Lin L1 (myFAttach,myDirAttach);
145   gp_Lin L2 (mySAttach,myDirAttach);
146   gp_Pnt Proj1 = ElCLib::Value(ElCLib::Parameter(L1,myPosition),L1);
147   gp_Pnt Proj2 = ElCLib::Value(ElCLib::Parameter(L2,myPosition),L2);
148   
149   gp_Lin L3;
150   Handle(SelectMgr_EntityOwner) own = new SelectMgr_EntityOwner(this,7);
151
152   if (!Proj1.IsEqual(Proj2,Precision::Confusion()))
153     {
154       L3 = gce_MakeLin(Proj1,Proj2);
155     }
156   else
157     {
158       L3 = gce_MakeLin(Proj1,myDirAttach);
159       Standard_Real size(Min(myVal/100.+1.e-6,myArrowSize+1.e-6));
160       Handle(Select3D_SensitiveBox) box =
161         new Select3D_SensitiveBox(own,
162                                   myPosition.X(),
163                                   myPosition.Y(),
164                                   myPosition.Z(),
165                                   myPosition.X()+size,
166                                   myPosition.Y()+size,
167                                   myPosition.Z()+size);
168       aSelection->Add(box);
169     }
170   Standard_Real parmin,parmax,parcur;
171   parmin = ElCLib::Parameter(L3,Proj1);
172   parmax = parmin;
173   
174   parcur = ElCLib::Parameter(L3,Proj2);
175   parmin = Min(parmin,parcur);
176   parmax = Max(parmax,parcur);
177   
178   parcur = ElCLib::Parameter(L3,myPosition);
179   parmin = Min(parmin,parcur);
180   parmax = Max(parmax,parcur);
181   
182   gp_Pnt PointMin = ElCLib::Value(parmin,L3);
183   gp_Pnt PointMax = ElCLib::Value(parmax,L3);
184
185   Handle(Select3D_SensitiveSegment) seg;
186   
187   if (!PointMin.IsEqual(PointMax,Precision::Confusion()))
188     {
189       seg = new Select3D_SensitiveSegment(own,
190                                           PointMin,
191                                           PointMax);
192       aSelection->Add(seg);
193     }
194   if (!myFAttach.IsEqual(Proj1,Precision::Confusion()))
195     {
196       seg = new Select3D_SensitiveSegment(own, myFAttach, Proj1);
197       aSelection->Add(seg);
198     }
199   if (!mySAttach.IsEqual(Proj2,Precision::Confusion()))
200     {
201       seg = new Select3D_SensitiveSegment(own, mySAttach, Proj2);
202       aSelection->Add(seg);
203     }
204 }
205
206 //=======================================================================
207 //function : ComputeTwoFacesParallel
208 //purpose  : 
209 //=======================================================================
210 void AIS_ParallelRelation::ComputeTwoFacesParallel(const Handle(Prs3d_Presentation)&)
211 {
212   Standard_NotImplemented::Raise("AIS_ParallelRelation::ComputeTwoFacesParallel not implemented");
213 }
214
215 //=======================================================================
216 //function : ComputeTwoEdgesParallel
217 //purpose  : 
218 //=======================================================================
219 void AIS_ParallelRelation::ComputeTwoEdgesParallel(const Handle(Prs3d_Presentation)& aPresentation)
220 {
221   TopoDS_Edge E1 = TopoDS::Edge(myFShape);
222   TopoDS_Edge E2 = TopoDS::Edge(mySShape);
223
224   gp_Pnt ptat11,ptat12,ptat21,ptat22;//,pint3d;
225   Handle(Geom_Curve) geom1,geom2;
226   Standard_Boolean isInfinite1,isInfinite2;
227   Handle(Geom_Curve) extCurv;
228   if (!AIS::ComputeGeometry(E1,E2,myExtShape,
229                             geom1,geom2,
230                             ptat11,ptat12,ptat21,ptat22,
231                             extCurv,
232                             isInfinite1,isInfinite2,
233                             myPlane))
234     {
235       return;
236     }
237
238   aPresentation->SetInfiniteState((isInfinite1 || isInfinite2) && (myExtShape != 0));
239
240   gp_Lin l1;
241   gp_Lin l2;
242   Standard_Boolean isEl1 = Standard_False, isEl2 = Standard_False;
243
244   if (geom1->IsInstance(STANDARD_TYPE(Geom_Ellipse)))
245     {
246       Handle(Geom_Ellipse) geom_el1 (Handle(Geom_Ellipse)::DownCast (geom1));
247       // construct lines through focuses
248       gp_Ax1 elAx = geom_el1->XAxis();
249       l1 = gp_Lin(elAx);
250       Standard_Real focex = geom_el1->MajorRadius() - geom_el1->Focal()/2.0;
251       gp_Vec transvec = gp_Vec(elAx.Direction())*focex;
252       ptat11 = geom_el1->Focus1().Translated(transvec);
253       ptat12 = geom_el1->Focus2().Translated(-transvec);
254       isEl1 = Standard_True;
255     }
256   else if (geom1->IsInstance(STANDARD_TYPE(Geom_Line)))
257     {
258       Handle(Geom_Line) geom_lin1 (Handle(Geom_Line)::DownCast (geom1));
259       l1 = geom_lin1->Lin();
260     }
261   else return;
262
263   if (geom2->IsInstance(STANDARD_TYPE(Geom_Ellipse)))
264     {
265       Handle(Geom_Ellipse) geom_el2 (Handle(Geom_Ellipse)::DownCast (geom2));
266       // construct lines through focuses
267       gp_Ax1 elAx = geom_el2->XAxis();
268       l2 = gp_Lin(elAx);
269       Standard_Real focex = geom_el2->MajorRadius() - geom_el2->Focal()/2.0;
270       gp_Vec transvec = gp_Vec(elAx.Direction())*focex;
271       ptat21 = geom_el2->Focus1().Translated(transvec);
272       ptat22 = geom_el2->Focus2().Translated(-transvec);
273       isEl2 = Standard_True;
274     }
275   else if (geom2->IsInstance(STANDARD_TYPE(Geom_Line)))
276     {
277       Handle(Geom_Line) geom_lin2 (Handle(Geom_Line)::DownCast (geom2));
278       l2 = geom_lin2->Lin();
279     }
280   else return;
281
282   const Handle(Geom_Line)& geom_lin1 = new Geom_Line(l1);
283   const Handle(Geom_Line)& geom_lin2 = new Geom_Line(l2);
284
285   myDirAttach = l1.Direction();
286   // size
287   if( !myArrowSizeIsDefined ) {
288     Standard_Real arrSize1 (myArrowSize), arrSize2 (myArrowSize);
289     if (!isInfinite1) arrSize1 = ptat11.Distance(ptat12)/50.;
290     if (!isInfinite2) arrSize2 = ptat21.Distance(ptat22)/50.;
291     myArrowSize = Max(myArrowSize,Max(arrSize1,arrSize2));
292 //  myArrowSize = Min(myArrowSize,Min(arrSize1,arrSize2));
293   }
294
295   if ( myAutomaticPosition )
296     {    
297       gp_Pnt curpos;
298       if ( !isInfinite1 )
299         {
300           gp_Pnt p2 = ElCLib::Value(ElCLib::Parameter(l2,ptat11),l2);
301           curpos.SetXYZ((ptat11.XYZ() + p2.XYZ())/2.);
302         }
303       else if ( !isInfinite2 )
304         {
305           gp_Pnt p2 = ElCLib::Value(ElCLib::Parameter(l1,ptat21),l1);
306           curpos.SetXYZ((ptat21.XYZ()+p2.XYZ())/2.);
307         }
308       else
309         {
310           curpos.SetXYZ((l1.Location().XYZ()+l2.Location().XYZ())/2.);
311         }
312       // offset pour eviter confusion Edge et Dimension
313       gp_Vec offset (myDirAttach);
314       offset = offset*myArrowSize*(-10.);
315       curpos.Translate(offset);
316       myPosition = curpos;
317     }
318
319   // recherche points attache
320   if (!isInfinite1)
321     {
322       if ( isEl1 )
323         {
324           if (myPosition.Distance(ptat11) < myPosition.Distance(ptat12)) myFAttach = ptat12;
325           else myFAttach = ptat11;
326         }
327       else
328         {
329           if (myPosition.Distance(ptat11) > myPosition.Distance(ptat12)) myFAttach = ptat12;
330           else myFAttach = ptat11;
331         }
332     }
333   else
334     {
335       myFAttach = ElCLib::Value(ElCLib::Parameter(l1,myPosition),l1);
336     }
337   
338   if (!isInfinite2)
339     {
340       if ( isEl2 )
341         {
342           if (myPosition.Distance(ptat21) < myPosition.Distance(ptat22)) mySAttach = ptat22;
343           else mySAttach = ptat21;
344         }
345       else
346         {
347           if (myPosition.Distance(ptat21) > myPosition.Distance(ptat22)) mySAttach = ptat22;
348           else mySAttach = ptat21;
349         }
350     }
351   else
352     {
353       mySAttach = ElCLib::Value(ElCLib::Parameter(l2,myPosition),l2);
354     }
355   TCollection_ExtendedString aText (" //");
356   
357   if (l1.Distance(l2) <= Precision::Confusion())
358     {
359       myArrowSize = 0.;
360     }
361   Handle(Prs3d_DimensionAspect) la = myDrawer->DimensionAspect();
362   Handle(Prs3d_ArrowAspect) arr = la->ArrowAspect();
363   arr->SetLength(myArrowSize);
364   arr = la->ArrowAspect();
365   arr->SetLength(myArrowSize);
366   if ( myExtShape == 1)
367     mySymbolPrs = DsgPrs_AS_FIRSTPT_LASTAR;
368   else if ( myExtShape == 2)
369     mySymbolPrs = DsgPrs_AS_FIRSTAR_LASTPT;
370
371   DsgPrs_LengthPresentation::Add(aPresentation,
372                                  myDrawer,
373                                  aText,
374                                  myFAttach,
375                                  mySAttach,
376                                  myDirAttach,
377                                  myPosition,
378                                  mySymbolPrs);
379   if ( (myExtShape != 0) &&  !extCurv.IsNull())
380     {
381       gp_Pnt pf, pl;
382       if ( myExtShape == 1 )
383         {
384           if (!isInfinite1)
385             {
386               pf = ptat11; 
387               pl = ptat12;
388             }
389           ComputeProjEdgePresentation(aPresentation,E1,geom_lin1,pf,pl);
390         }
391       else
392         {
393           if (!isInfinite2)
394             {
395               pf = ptat21; 
396               pl = ptat22;
397             }
398           ComputeProjEdgePresentation(aPresentation,E2,geom_lin2,pf,pl);
399         }
400     }
401 }