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