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