0031459: Visualization, AIS_TextLabel - add missing getters
[occt.git] / src / PrsDim / PrsDim_EqualDistanceRelation.cxx
1 // Created on: 1998-01-24
2 // Created by: Julia GERASIMOVA
3 // Copyright (c) 1998-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_EqualDistanceRelation.hxx>
18
19 #include <PrsDim.hxx>
20 #include <PrsDim_LengthDimension.hxx>
21 #include <Bnd_Box.hxx>
22 #include <BRep_Tool.hxx>
23 #include <BRepAdaptor_Curve.hxx>
24 #include <BRepBuilderAPI_MakeEdge.hxx>
25 #include <BRepBuilderAPI_MakeVertex.hxx>
26 #include <DsgPrs_EqualDistancePresentation.hxx>
27 #include <ElCLib.hxx>
28 #include <Geom_Circle.hxx>
29 #include <Geom_Line.hxx>
30 #include <Geom_Plane.hxx>
31 #include <GeomAPI_ProjectPointOnCurve.hxx>
32 #include <GeomAPI_ProjectPointOnSurf.hxx>
33 #include <gp_Lin.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_SensitiveCircle.hxx>
42 #include <Select3D_SensitiveSegment.hxx>
43 #include <SelectMgr_EntityOwner.hxx>
44 #include <Standard_NotImplemented.hxx>
45 #include <TopoDS_Edge.hxx>
46 #include <TopoDS_Shape.hxx>
47 #include <TopoDS_Vertex.hxx>
48
49 IMPLEMENT_STANDARD_RTTIEXT(PrsDim_EqualDistanceRelation, PrsDim_Relation)
50
51 //=======================================================================
52 //function : PrsDim_EqualDistanceRelation
53 //purpose  : 
54 //=======================================================================
55 PrsDim_EqualDistanceRelation::PrsDim_EqualDistanceRelation( const TopoDS_Shape& aShape1,
56                                                       const TopoDS_Shape& aShape2,
57                                                       const TopoDS_Shape& aShape3,
58                                                       const TopoDS_Shape& aShape4,
59                                                       const Handle( Geom_Plane )& aPlane )
60 : PrsDim_Relation()
61 {
62   myFShape = aShape1;
63   mySShape = aShape2;
64   myShape3 = aShape3;
65   myShape4 = aShape4;
66   myPlane  = aPlane;
67
68   //Temporary
69   myArrowSize = 3.0; //set the concrete value
70   mySymbolPrs = DsgPrs_AS_BOTHAR;
71 }
72
73 //=======================================================================
74 //function : Compute
75 //purpose  : 
76 //=======================================================================
77
78 void PrsDim_EqualDistanceRelation::Compute( const Handle( PrsMgr_PresentationManager3d )&,
79                                          const Handle( Prs3d_Presentation )& aPresentation,
80                                          const Standard_Integer ) 
81 {
82   gp_Pnt Position12 =  myPosition, Position34 = myPosition;
83
84   Handle(Prs3d_DimensionAspect) la = myDrawer->DimensionAspect();
85   Handle(Prs3d_ArrowAspect) arr = la->ArrowAspect();
86   arr->SetLength(myArrowSize);
87 // -- ota -- begin
88   if (!myAutomaticPosition ){
89     gp_Pnt aMiddle12 ((myPoint1.XYZ() + myPoint2.XYZ())*0.5);
90     gp_Pnt aMiddle34 ((myPoint3.XYZ() + myPoint4.XYZ())*0.5);
91
92     if (myPosition.Distance(aMiddle12) > myPosition.Distance(aMiddle34))
93       Position12.SetXYZ((myPoint1.XYZ() + myPoint2.XYZ())*0.5);
94     else        
95       Position34.SetXYZ((myPoint3.XYZ() + myPoint4.XYZ())*0.5);
96
97   }
98   
99   if (myFShape.ShapeType() == TopAbs_EDGE && mySShape.ShapeType() == TopAbs_EDGE)
100     PrsDim_EqualDistanceRelation::ComputeTwoEdgesLength(aPresentation,
101                                                      myDrawer,
102                                                      myArrowSize,
103                                                      TopoDS::Edge(myFShape),
104                                                      TopoDS::Edge(mySShape),
105                                                      myPlane,
106                                                      myAutomaticPosition,
107                                                      myIsSetBndBox,
108                                                      myBndBox,
109                                                      Position12,
110                                                      myAttachPoint1,
111                                                      myAttachPoint2,
112                                                      myPoint1,
113                                                      myPoint2,
114                                                      mySymbolPrs );
115          
116   
117   else if (myFShape.ShapeType() == TopAbs_VERTEX && mySShape.ShapeType() == TopAbs_VERTEX)
118     PrsDim_EqualDistanceRelation::ComputeTwoVerticesLength( aPresentation,
119                                                         myDrawer,
120                                                         myArrowSize,
121                                                         TopoDS::Vertex(myFShape),
122                                                         TopoDS::Vertex(mySShape),
123                                                         myPlane,
124                                                         myAutomaticPosition,
125                                                         myIsSetBndBox,
126                                                         myBndBox,
127                                                         PrsDim_TypeOfDist_Unknown,
128                                                         Position12,
129                                                         myAttachPoint1,
130                                                         myAttachPoint2,
131                                                         myPoint1,
132                                                         myPoint2,
133                                                         mySymbolPrs );
134   else 
135     PrsDim_EqualDistanceRelation::ComputeOneEdgeOneVertexLength( aPresentation,
136                                                              myDrawer,
137                                                              myArrowSize,
138                                                              myFShape,
139                                                              mySShape,
140                                                              myPlane,
141                                                              myAutomaticPosition,
142                                                              myIsSetBndBox,
143                                                              myBndBox,
144                                                              Position12,
145                                                              myAttachPoint1,
146                                                              myAttachPoint2,
147                                                              myPoint1,
148                                                              myPoint2,
149                                                              mySymbolPrs );
150
151   if (myShape3.ShapeType() == TopAbs_EDGE && myShape4.ShapeType() == TopAbs_EDGE)
152     PrsDim_EqualDistanceRelation::ComputeTwoEdgesLength(aPresentation,
153                                                      myDrawer,
154                                                      myArrowSize,
155                                                      TopoDS::Edge(myShape3),
156                                                      TopoDS::Edge(myShape4),
157                                                      myPlane,
158                                                      myAutomaticPosition,
159                                                      myIsSetBndBox,
160                                                      myBndBox,
161                                                      Position34,
162                                                      myAttachPoint3,
163                                                      myAttachPoint4,
164                                                      myPoint3,
165                                                      myPoint4,
166                                                      mySymbolPrs );
167   
168   else if (myShape3.ShapeType() == TopAbs_VERTEX && myShape4.ShapeType() == TopAbs_VERTEX)
169     PrsDim_EqualDistanceRelation::ComputeTwoVerticesLength( aPresentation,
170                                                         myDrawer,
171                                                         myArrowSize,
172                                                         TopoDS::Vertex(myShape3),
173                                                         TopoDS::Vertex(myShape4),
174                                                         myPlane,
175                                                         myAutomaticPosition,
176                                                         myIsSetBndBox,
177                                                         myBndBox,
178                                                         PrsDim_TypeOfDist_Unknown,
179                                                         Position34,
180                                                         myAttachPoint3,
181                                                         myAttachPoint4,
182                                                         myPoint3,
183                                                         myPoint4,
184                                                         mySymbolPrs );
185    
186   else
187     PrsDim_EqualDistanceRelation::ComputeOneEdgeOneVertexLength( aPresentation,
188                                                              myDrawer,
189                                                              myArrowSize,
190                                                              myShape3,
191                                                              myShape4,
192                                                              myPlane,
193                                                              myAutomaticPosition,
194                                                              myIsSetBndBox,
195                                                              myBndBox,
196                                                              Position34,
197                                                              myAttachPoint3,
198                                                              myAttachPoint4,
199                                                              myPoint3,
200                                                              myPoint4,
201                                                              mySymbolPrs );
202   
203   DsgPrs_EqualDistancePresentation::Add( aPresentation, myDrawer, 
204                                         myPoint1, myPoint2, myPoint3, myPoint4, myPlane );
205 }
206
207 //=======================================================================
208 //function : ComputeSelection
209 //purpose  : 
210 //=======================================================================
211
212 void PrsDim_EqualDistanceRelation::ComputeSelection( const Handle( SelectMgr_Selection )& aSelection,
213                                                   const Standard_Integer ) 
214 {
215   Handle( SelectMgr_EntityOwner ) own = new SelectMgr_EntityOwner( this, 7 );
216   Handle( Select3D_SensitiveSegment ) seg;
217
218   seg = new Select3D_SensitiveSegment( own, myPoint1, myPoint2 );
219   aSelection->Add( seg );
220
221   seg = new Select3D_SensitiveSegment( own, myPoint3, myPoint4 );
222   aSelection->Add( seg );
223
224   // Line between two middles
225   gp_Pnt Middle12( (myPoint1.XYZ() + myPoint2.XYZ()) * 0.5 ), 
226          Middle34( (myPoint3.XYZ() + myPoint4.XYZ()) *0.5 );
227   seg = new Select3D_SensitiveSegment( own, Middle12, Middle34 );
228   aSelection->Add( seg );
229
230   gp_Pnt Middle((Middle12.XYZ() +  Middle34.XYZ())*0.5);
231   Standard_Real SmallDist = .001;
232   Handle( Select3D_SensitiveBox ) box = new Select3D_SensitiveBox( own,
233                                                                   Middle.X() - SmallDist,
234                                                                   Middle.Y() - SmallDist,
235                                                                   Middle.Z() - SmallDist,
236                                                                   Middle.X() + SmallDist,
237                                                                   Middle.Y() + SmallDist,
238                                                                   Middle.Z() + SmallDist );
239   aSelection->Add(box);
240
241   if (myFShape.ShapeType() == TopAbs_EDGE){
242     BRepAdaptor_Curve aCurve(TopoDS::Edge(myFShape));
243     if (aCurve.GetType() == GeomAbs_Line){
244         //add sensetive element - line
245       seg = new Select3D_SensitiveSegment( own, myAttachPoint1, myPoint1);
246       aSelection->Add( seg );
247     }
248     else if (aCurve.GetType() == GeomAbs_Circle){
249       Handle(Geom_Circle) aCircle = Handle(Geom_Circle)::DownCast(aCurve.Curve().Curve());
250       Standard_Real FirstPar = ElCLib::Parameter(aCircle->Circ(), myAttachPoint1),
251                     LastPar  = ElCLib::Parameter(aCircle->Circ(), myPoint1);
252       if (LastPar < FirstPar ) LastPar+=M_PI*2;
253       //add sensetive arc
254       Handle(Select3D_SensitiveCircle) circ = 
255         new Select3D_SensitiveCircle( own, aCircle,  FirstPar, LastPar);
256       aSelection->Add( circ );
257     }
258   }
259   else {
260       seg = new Select3D_SensitiveSegment( own, myAttachPoint1, myPoint1);
261       aSelection->Add( seg );
262     }
263   
264   if (mySShape.ShapeType() == TopAbs_EDGE){
265     BRepAdaptor_Curve aCurve(TopoDS::Edge(mySShape));
266     if (aCurve.GetType() == GeomAbs_Line) {
267       //add sensetive element - line
268       seg = new Select3D_SensitiveSegment( own, myAttachPoint2, myPoint2);
269       aSelection->Add( seg );
270     }
271     else if (aCurve.GetType() == GeomAbs_Circle){
272       Handle(Geom_Circle) aCircle = Handle(Geom_Circle)::DownCast(aCurve.Curve().Curve());
273       Standard_Real FirstPar = ElCLib::Parameter(aCircle->Circ(), myAttachPoint2),
274       LastPar  = ElCLib::Parameter(aCircle->Circ(), myPoint2);
275       if (LastPar < FirstPar ) LastPar+=M_PI*2;
276       //add sensetive arc
277       Handle(Select3D_SensitiveCircle) circ = 
278         new Select3D_SensitiveCircle( own,aCircle,  FirstPar, LastPar);
279       aSelection->Add( circ );
280     }
281   }
282   else {
283     seg = new Select3D_SensitiveSegment( own, myAttachPoint2, myPoint2);
284     aSelection->Add( seg );
285   }
286     
287   if (myShape3.ShapeType() == TopAbs_EDGE){
288     BRepAdaptor_Curve aCurve(TopoDS::Edge(myShape3));
289     if (aCurve.GetType() == GeomAbs_Line) {
290       //add sensetive element - line
291       seg = new Select3D_SensitiveSegment( own, myAttachPoint3, myPoint3);
292       aSelection->Add( seg );
293     }
294     else if (aCurve.GetType() == GeomAbs_Circle){
295       Handle(Geom_Circle) aCircle = Handle(Geom_Circle)::DownCast(aCurve.Curve().Curve());
296       Standard_Real FirstPar = ElCLib::Parameter(aCircle->Circ(), myAttachPoint3),
297       LastPar  = ElCLib::Parameter(aCircle->Circ(), myPoint3);
298       if (LastPar < FirstPar ) LastPar+=M_PI*2;
299       Handle(Select3D_SensitiveCircle) circ = 
300         new Select3D_SensitiveCircle( own, aCircle,  FirstPar, LastPar);
301       aSelection->Add( circ );
302     }
303     else {
304       seg = new Select3D_SensitiveSegment( own, myAttachPoint3, myPoint3);
305       aSelection->Add( seg );
306     }
307   }
308   else {
309     seg = new Select3D_SensitiveSegment( own, myAttachPoint3, myPoint3);
310     aSelection->Add( seg );
311   }
312
313   if (myShape4.ShapeType() == TopAbs_EDGE){
314     BRepAdaptor_Curve aCurve(TopoDS::Edge(myShape4));
315     if (aCurve.GetType() == GeomAbs_Line) {
316       //add sensetive element - line
317       seg = new Select3D_SensitiveSegment( own, myAttachPoint4, myPoint4);
318       aSelection->Add( seg );
319     }
320     else if (aCurve.GetType() == GeomAbs_Circle){
321       Handle(Geom_Circle) aCircle = Handle(Geom_Circle)::DownCast(aCurve.Curve().Curve());
322       Standard_Real FirstPar = ElCLib::Parameter(aCircle->Circ(), myAttachPoint4),
323       LastPar  = ElCLib::Parameter(aCircle->Circ(), myPoint4);
324       if (LastPar < FirstPar ) LastPar+=M_PI*2;
325       //add sensetive arc
326       Handle(Select3D_SensitiveCircle) circ = 
327         new Select3D_SensitiveCircle( own,aCircle,  FirstPar, LastPar);
328       aSelection->Add( circ );
329     }
330   }
331   else {
332     seg = new Select3D_SensitiveSegment( own, myAttachPoint4, myPoint4);
333     aSelection->Add( seg );
334   }
335 }
336
337 //=======================================================================
338 //function : ComputeTwoEdgesLength
339 //purpose  : 
340 //=======================================================================
341 void PrsDim_EqualDistanceRelation::ComputeTwoEdgesLength( const Handle( Prs3d_Presentation )& aPresentation,
342                                                       const Handle( Prs3d_Drawer )& aDrawer,
343                                                       const Standard_Real ArrowSize,
344                                                       const TopoDS_Edge & FirstEdge,
345                                                       const TopoDS_Edge & SecondEdge,
346                                                       const Handle( Geom_Plane )& Plane,
347                                                       const Standard_Boolean AutomaticPos,
348                                                       const Standard_Boolean IsSetBndBox,
349                                                       const Bnd_Box & BndBox,
350                                                       gp_Pnt& Position,
351                                                       gp_Pnt& FirstAttach,
352                                                       gp_Pnt& SecondAttach,
353                                                       gp_Pnt& FirstExtreme,
354                                                       gp_Pnt& SecondExtreme, 
355                                                       DsgPrs_ArrowSide & SymbolPrs )
356
357   gp_Dir DirAttach;
358   BRepAdaptor_Curve cu1( FirstEdge );
359   BRepAdaptor_Curve cu2( SecondEdge );
360   
361   // 3d lines
362   Handle(Geom_Curve) geom1,geom2;
363   gp_Pnt ptat11,ptat12,ptat21,ptat22;
364   
365   Standard_Boolean isInfinite1(Standard_False),isInfinite2(Standard_False);
366   Handle(Geom_Curve) extCurv;
367   Standard_Real arrsize = ArrowSize;// size
368   Standard_Real Val=0.;
369   Standard_Boolean isInPlane1, isInPlane2;
370
371   if(!PrsDim::ComputeGeometry(FirstEdge,geom1, ptat11, ptat12,extCurv,isInfinite1,isInPlane1, Plane ))
372     return;
373   if(!PrsDim::ComputeGeometry(SecondEdge, geom2, ptat21, ptat22, extCurv, isInfinite2,isInPlane2, Plane))
374     return;
375   
376   aPresentation->SetInfiniteState(isInfinite1 || isInfinite2);
377   
378   if (cu1.GetType() == GeomAbs_Line && cu2.GetType() == GeomAbs_Line) 
379     {
380       Handle(Geom_Line) geom_lin1 (Handle(Geom_Line)::DownCast (geom1));
381       Handle(Geom_Line) geom_lin2 (Handle(Geom_Line)::DownCast (geom2));
382       const gp_Lin& l1 = geom_lin1->Lin();
383       const gp_Lin& l2 = geom_lin2->Lin();
384       
385       //Get Val value
386       Val = l1.Distance( l2 );
387       
388       DirAttach = l1.Direction();
389       
390       if (AutomaticPos) {
391         // compute position of offset point
392         gp_Pnt curpos;
393         Standard_Real par1=0., par2=0.;
394         if(!(isInfinite1 || isInfinite2))
395           {
396             par1 = ElCLib::Parameter(l1,ptat11);
397             par2 = ElCLib::Parameter(l1,ptat21);
398             if (par1 <par2){//project ptat11 on l2
399               gp_Pnt p2 = ElCLib::Value(ElCLib::Parameter(l2,ptat11),l2);
400               curpos.SetXYZ((ptat11.XYZ()+p2.XYZ())*0.5);
401             }
402             else {//project ptat21 on l1
403               gp_Pnt p2 = ElCLib::Value(par2, l1);
404               curpos.SetXYZ((ptat21.XYZ()+p2.XYZ())*0.5);
405             }
406           }
407         else if (!isInfinite1){
408           par2 = ElCLib::Parameter(l1,ptat21);
409           gp_Pnt p2 = ElCLib::Value(par2,l1);
410           curpos.SetXYZ((ptat21.XYZ()+p2.XYZ())/2.);
411         }
412         else if (!isInfinite2) {
413           gp_Pnt p2 = ElCLib::Value(ElCLib::Parameter(l2,ptat11),l2);
414           curpos.SetXYZ((ptat11.XYZ()+p2.XYZ())*0.5);
415         }
416         else   
417           curpos.SetXYZ((l1.Location().XYZ()+l2.Location().XYZ())*0.5);
418  
419     // compute  offset
420     gp_Vec offset(DirAttach);
421     offset = offset*ArrowSize*(-10.);
422     curpos.Translate(offset);
423     Position = curpos;
424   }
425   else {    // project point on the plane
426     Position = PrsDim::ProjectPointOnPlane( Position, Plane->Pln() );
427   }
428
429   // find attach points
430   if (!isInfinite1) {
431     if (Position.Distance(ptat11) > Position.Distance(ptat12)) FirstAttach = ptat12;
432     else FirstAttach = ptat11;
433   }
434   else {
435     FirstAttach = ElCLib::Value(ElCLib::Parameter(l1,Position),l1);
436   }
437   
438   if (!isInfinite2) {
439     if (Position.Distance(ptat21) > Position.Distance(ptat22)) SecondAttach = ptat22;
440     else SecondAttach = ptat21;
441   }
442   else {
443     SecondAttach = ElCLib::Value(ElCLib::Parameter(l2,Position),l2);
444   }
445
446   Standard_Real confusion(Precision::Confusion());
447   if (arrsize < confusion) arrsize = Val*0.1;
448   if (Abs(Val) <= confusion) {arrsize = 0.;}
449
450   Handle(Prs3d_DimensionAspect) la = aDrawer->DimensionAspect();
451   Handle(Prs3d_ArrowAspect) arr = la->ArrowAspect();  
452   arr->SetLength(arrsize);
453   arr = la->ArrowAspect();
454   arr->SetLength(arrsize);
455
456   if (AutomaticPos && IsSetBndBox)
457     Position = PrsDim::TranslatePointToBound( Position, DirAttach, BndBox );
458  
459    DsgPrs_EqualDistancePresentation::AddInterval(aPresentation,
460                                                 aDrawer,
461                                                 FirstAttach,
462                                                 SecondAttach,
463                                                 DirAttach,
464                                                 Position,
465                                                 SymbolPrs,
466                                                 FirstExtreme,
467                                                 SecondExtreme);
468                                                 
469       
470 }
471   if (cu1.GetType() == GeomAbs_Circle && cu2.GetType() == GeomAbs_Circle){
472     //Get first and last points of circles
473     Handle(Geom_Circle) aCir1 (Handle(Geom_Circle)::DownCast(geom1));
474     Handle(Geom_Circle) aCir2 (Handle(Geom_Circle)::DownCast(geom2));
475     gp_Circ aCirc1 = aCir1->Circ();
476     gp_Circ aCirc2 = aCir2->Circ();
477
478     //To avoid circles with different orientaion
479     Standard_Real aTol = Precision::Confusion();
480     if(aCirc2.Axis().IsOpposite(aCirc1.Axis(), aTol) ||
481        aCirc2.XAxis().IsOpposite(aCirc1.XAxis(), aTol) || 
482        aCirc2.YAxis().IsOpposite(aCirc1.YAxis(), aTol) )
483       {
484         aCirc2.SetPosition(aCirc1.Position());
485         aCirc2.SetAxis(aCirc1.Axis());
486       }
487     
488     if (AutomaticPos){ 
489       Standard_Real par1 = 0, par2 = 0;
490       gp_Pln aPln =  Plane->Pln();
491       //Project ptat12 and ptat22 on constraint plane
492       gp_Pnt PrPnt12 = PrsDim::ProjectPointOnPlane(ptat12, aPln);
493       gp_Pnt PrPnt22 = PrsDim::ProjectPointOnPlane(ptat22, aPln);
494       //Project circles center on constraint plane
495       gp_Pnt PrCenter = PrsDim::ProjectPointOnPlane(aCirc1.Location(), aPln);
496
497       gp_Dir XDir = aPln.XAxis().Direction();
498       gp_Dir YDir = aPln.YAxis().Direction();
499       
500       
501       if (PrPnt12.Distance(PrCenter) >Precision::Confusion())
502         {
503           gp_Dir aDir1(PrPnt12.XYZ() - PrCenter.XYZ());
504           Standard_Real anAngle = aDir1.Angle(XDir); //Get the angle in range [0, M_PI]
505           if (aDir1.Dot(YDir) < 0)
506             anAngle = 2 * M_PI - anAngle;
507           par1 = anAngle;
508         }
509       
510       if (PrPnt22.Distance(PrCenter) >Precision::Confusion())
511         {
512           gp_Dir aDir2(PrPnt22.XYZ() - PrCenter.XYZ());
513           Standard_Real anAngle = aDir2.Angle(XDir); //Get the angle in range [0, M_PI]
514           if (aDir2.Dot(YDir) < 0)
515             anAngle = 2 * M_PI - anAngle;
516           par2 = anAngle;
517         }
518       
519       
520       if(par1 > par2 ){
521         FirstExtreme = ptat12;
522         Standard_Real aPar1 = ElCLib::Parameter(aCirc2, ptat12);
523         SecondExtreme = ElCLib::Value(aPar1, aCirc2);
524       }
525       else {
526         Standard_Real aPar2 =  ElCLib::Parameter(aCirc1, ptat22);
527         FirstExtreme = ElCLib::Value(aPar2, aCirc1);
528         SecondExtreme = ptat22;
529       }
530     }
531     else {
532       Standard_Real pospar = ElCLib::Parameter(aCirc1, Position);
533       FirstExtreme  =  ElCLib::Value(pospar, aCirc1);
534       pospar = ElCLib::Parameter(aCirc2, Position);
535       SecondExtreme =  ElCLib::Value(pospar, aCirc2);
536     }
537
538     DsgPrs_EqualDistancePresentation::AddIntervalBetweenTwoArcs(aPresentation,
539                                                                 aDrawer,
540                                                                 aCirc1,
541                                                                 aCirc2,
542                                                                 ptat12,
543                                                                 FirstExtreme, 
544                                                                 ptat22,
545                                                                 SecondExtreme,
546                                                                 SymbolPrs);
547     
548     FirstAttach = ptat12; SecondAttach = ptat22; //assign attach points
549     Position.SetXYZ( (FirstAttach.XYZ() + SecondAttach.XYZ())*0.5);
550   }
551
552   if (arrsize < Precision::Confusion()) arrsize = Val*0.1;
553   if (Abs(Val) <=  Precision::Confusion()) {arrsize = 0.;}
554
555 //  gp_Pnt pf, pl;
556   if (!isInPlane1) {
557     PrsDim::ComputeProjEdgePresentation( aPresentation, aDrawer, FirstEdge, geom1, ptat11, ptat12);
558   }
559   if(!isInPlane2) {
560     PrsDim::ComputeProjEdgePresentation( aPresentation, aDrawer, SecondEdge, geom2, ptat21, ptat22);
561   }
562 }
563
564 //=======================================================================
565 //function : ComputeTwoVerticesLength
566 //purpose  : 
567 //=======================================================================
568
569 void PrsDim_EqualDistanceRelation::ComputeTwoVerticesLength( const Handle( Prs3d_Presentation )& aPresentation,
570                                                          const Handle( Prs3d_Drawer )& aDrawer,
571                                                          const Standard_Real ArrowSize,
572                                                          const TopoDS_Vertex& FirstVertex,
573                                                          const TopoDS_Vertex& SecondVertex,
574                                                          const Handle( Geom_Plane )& Plane,
575                                                          const Standard_Boolean AutomaticPos,
576                                                          const Standard_Boolean IsSetBndBox,
577                                                          const Bnd_Box& BndBox,
578                                                          const PrsDim_TypeOfDist TypeDist,
579                                                          gp_Pnt& Position,
580                                                          gp_Pnt& FirstAttach,
581                                                          gp_Pnt& SecondAttach,
582                                                          gp_Pnt& FirstExtreme,
583                                                          gp_Pnt& SecondExtreme, 
584                                                          DsgPrs_ArrowSide& SymbolPrs )
585 {
586   Standard_Boolean isOnPlane1, isOnPlane2;
587   gp_Dir DirAttach;
588   PrsDim::ComputeGeometry( FirstVertex, FirstAttach, Plane, isOnPlane1);
589   PrsDim::ComputeGeometry( SecondVertex, SecondAttach, Plane, isOnPlane2);
590
591   Standard_Real confusion(Precision::Confusion());
592   Standard_Boolean samePoint(FirstAttach.IsEqual(SecondAttach,confusion));
593
594   if (TypeDist == PrsDim_TypeOfDist_Vertical) DirAttach =  Plane->Pln().XAxis().Direction();
595   else if (TypeDist == PrsDim_TypeOfDist_Horizontal) DirAttach =  Plane->Pln().YAxis().Direction();
596   else {
597     if (!samePoint) {
598       DirAttach.SetXYZ(SecondAttach.XYZ() - FirstAttach.XYZ());
599       DirAttach.Rotate(Plane->Pln().Axis(),M_PI/2.);
600     }
601   }
602   
603   // size
604   if (AutomaticPos) {
605    if (!samePoint) {
606      gp_Pnt curpos((FirstAttach.XYZ()+SecondAttach.XYZ())*0.5);
607      // make offset of curpos
608      gp_Vec offset(DirAttach);
609      offset = offset*ArrowSize*(-10.);
610      curpos.Translate(offset);
611      Position = curpos;
612    }
613    else {
614      gp_Dir aDir = Plane->Pln().Axis().Direction();
615      gp_Vec aVec (aDir.XYZ()*10*ArrowSize);
616      //Position = gp_Pnt(FirstAttach.XYZ()+gp_XYZ(1.,1.,1.)); // not correct
617      Position = FirstAttach.Translated(aVec);
618      Position = PrsDim::ProjectPointOnPlane( Position, Plane->Pln() );//not needed really
619      DirAttach.SetXYZ(Position.XYZ() - FirstAttach.XYZ());
620    }
621   }
622   else {   
623     Position = PrsDim::ProjectPointOnPlane( Position, Plane->Pln() );
624   }
625
626  
627   Handle(Prs3d_DimensionAspect) la = aDrawer->DimensionAspect();
628   Handle(Prs3d_ArrowAspect) arr = la->ArrowAspect();  
629   arr->SetLength(ArrowSize);
630   arr = la->ArrowAspect();
631   arr->SetLength(ArrowSize);
632
633   if (AutomaticPos && IsSetBndBox)
634     Position = PrsDim::TranslatePointToBound( Position, DirAttach, BndBox );
635
636   DsgPrs_EqualDistancePresentation::AddInterval(aPresentation,
637                                                 aDrawer,
638                                                 FirstAttach,
639                                                 SecondAttach,
640                                                 DirAttach,
641                                                 Position,
642                                                 SymbolPrs,
643                                                 FirstExtreme, //returned
644                                                 SecondExtreme); //returned
645
646   // Compute projection
647   if ( !isOnPlane1)
648     PrsDim::ComputeProjVertexPresentation(aPresentation, aDrawer, FirstVertex, FirstAttach);
649   if ( !isOnPlane2)
650     PrsDim::ComputeProjVertexPresentation(aPresentation, aDrawer, SecondVertex, SecondAttach);
651
652
653
654
655 //=======================================================================
656 //function : ComputeOneEdgeOneVertexLength
657 //purpose  : 
658 //=======================================================================
659
660 void PrsDim_EqualDistanceRelation::ComputeOneEdgeOneVertexLength( const Handle( Prs3d_Presentation )& aPresentation,
661                                                               const Handle( Prs3d_Drawer )& aDrawer,
662                                                               const Standard_Real ArrowSize,
663                                                               const TopoDS_Shape & FirstShape,
664                                                               const TopoDS_Shape & SecondShape,
665                                                               const Handle( Geom_Plane )& Plane,
666                                                               const Standard_Boolean AutomaticPos,
667                                                               const Standard_Boolean IsSetBndBox,
668                                                               const Bnd_Box & BndBox,
669                                                               gp_Pnt & Position,
670                                                               gp_Pnt & FirstAttach,
671                                                               gp_Pnt & SecondAttach,
672                                                               gp_Pnt& FirstExtreme,
673                                                               gp_Pnt& SecondExtreme, 
674                                                               DsgPrs_ArrowSide & SymbolPrs )
675 {
676   TopoDS_Vertex thevertex;
677   TopoDS_Edge theedge;
678   
679   gp_Pnt ptonedge1,ptonedge2;
680   Handle(Geom_Curve) aCurve;
681   Handle(Geom_Curve) extCurv;
682   Standard_Boolean isInfinite;
683   Standard_Real Val;
684   Standard_Boolean isOnPlanEdge, isOnPlanVertex;
685   Standard_Integer edgenum ;
686
687   if (FirstShape.ShapeType() == TopAbs_VERTEX) {  
688     thevertex = TopoDS::Vertex(FirstShape);
689     theedge   = TopoDS::Edge(SecondShape);
690     edgenum   = 2; //edge is the second shape
691   }
692   else {
693     thevertex = TopoDS::Vertex(SecondShape);
694     theedge   = TopoDS::Edge(FirstShape);
695     edgenum   = 1;//edge is the first shape
696   }
697   if (!PrsDim::ComputeGeometry(theedge,aCurve,ptonedge1,ptonedge2,extCurv,isInfinite,isOnPlanEdge,Plane))
698     return;
699   aPresentation->SetInfiniteState(isInfinite);
700   PrsDim::ComputeGeometry(thevertex, FirstAttach, Plane, isOnPlanVertex);
701
702   if ( aCurve->IsInstance(STANDARD_TYPE(Geom_Line)) ) 
703     {
704   Handle(Geom_Line) geom_lin (Handle(Geom_Line)::DownCast (aCurve));
705   const gp_Lin& l = geom_lin->Lin();
706
707   // computation of Val
708   Val = l.Distance( FirstAttach );
709
710   gp_Dir DirAttach = l.Direction();
711   // size
712   Standard_Real arrsize = ArrowSize;
713   if (Abs(Val) <= Precision::Confusion()) {arrsize = 0.;}
714
715   if (AutomaticPos) {
716     gp_Pnt p = ElCLib::Value(ElCLib::Parameter(l,FirstAttach),l);
717     gp_Pnt curpos((FirstAttach.XYZ()+p.XYZ())*0.5);
718     // make offset 
719     gp_Vec offset(DirAttach);
720     offset = offset*ArrowSize*(-10.);
721     curpos.Translate(offset);
722     Position = curpos;
723   }
724   else { // project point on the plane
725     Position = PrsDim::ProjectPointOnPlane( Position, Plane->Pln() );
726   }
727   
728   if (!isInfinite) {
729     if (Position.Distance(ptonedge1) > Position.Distance(ptonedge2)) SecondAttach = ptonedge2;
730     else SecondAttach = ptonedge1;
731   }
732   else {
733     SecondAttach = ElCLib::Value(ElCLib::Parameter(l,Position),l);
734   }
735   
736   Handle(Prs3d_DimensionAspect) la = aDrawer->DimensionAspect();
737   Handle(Prs3d_ArrowAspect) arr = la->ArrowAspect();  
738   arr->SetLength(arrsize);
739   arr = la->ArrowAspect();
740   arr->SetLength(arrsize);
741
742   if (AutomaticPos && IsSetBndBox)
743     Position = PrsDim::TranslatePointToBound( Position, DirAttach, BndBox );
744   DsgPrs_EqualDistancePresentation::AddInterval(aPresentation,
745                                                 aDrawer,
746                                                 FirstAttach,
747                                                 SecondAttach,
748                                                 DirAttach,
749                                                 Position,
750                                                 SymbolPrs,
751                                                 FirstExtreme,
752                                                 SecondExtreme);
753   
754 }
755   if (aCurve->IsInstance(STANDARD_TYPE(Geom_Circle))){
756     gp_Circ aCirc1 = (Handle(Geom_Circle)::DownCast(aCurve))->Circ();
757     gp_Circ aCirc2(aCirc1); aCirc2.SetRadius(0); //create the second formal circle
758     if(AutomaticPos)
759       {
760         SecondAttach = ptonedge2; //a vertex
761         Position.SetXYZ((SecondAttach.XYZ() + aCirc1.Location().XYZ())*0.5);
762       }
763     else {
764       Standard_Real aPar =  ElCLib::Parameter(aCirc1, Position);
765       SecondAttach =  ElCLib::Value(aPar, aCirc1);
766     }
767
768     Handle(Geom_Circle) aCurve2 = new Geom_Circle(aCirc2);
769     DsgPrs_EqualDistancePresentation::AddIntervalBetweenTwoArcs(aPresentation,
770                                                                 aDrawer,
771                                                                 aCirc1, //circle or arc
772                                                                 aCirc2, //really vertex
773                                                                 ptonedge2, //last point of aCirc1
774                                                                 SecondAttach,
775                                                                 FirstAttach, //vertex really
776                                                                 FirstAttach, 
777                                                                 SymbolPrs);
778     
779     //Assign arc points 
780     if (edgenum == 1){
781       FirstExtreme = SecondAttach; SecondExtreme = FirstAttach;
782       SecondAttach = FirstAttach; FirstAttach = ptonedge2; 
783     } else { //vertex is the first shape, circle is sthe last.
784       FirstExtreme = FirstAttach; SecondExtreme = SecondAttach;
785       SecondAttach = ptonedge2;
786     }
787   }
788  
789   // computation of Val
790   Val = FirstAttach.Distance(SecondAttach);
791   
792   //Display the pieces of attached to the curve if it is not 
793   // in the WP
794   if (!isOnPlanEdge) { // add presentation of projection of the edge in WP
795     PrsDim::ComputeProjEdgePresentation(aPresentation,aDrawer,theedge,aCurve,ptonedge1,ptonedge2);
796       }
797   if (!isOnPlanVertex) { // add presentation of projection of the vertex in WP 
798     PrsDim::ComputeProjVertexPresentation(aPresentation,aDrawer,thevertex,FirstAttach);
799     }
800   
801 }
802 // -- ota -- end