src/DsgPrs/DsgPrs_AnglePresentation.cxx

   1 // Created on: 1995-02-07
   2 // Copyright (c) 1995-1999 Matra Datavision
   3 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   4 //
   5 // This file is part of Open CASCADE Technology software library.
   6 //
   7 // This library is free software; you can redistribute it and/or modify it under
   8 // the terms of the GNU Lesser General Public License version 2.1 as published
   9 // by the Free Software Foundation, with special exception defined in the file
  10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  11 // distribution for complete text of the license and disclaimer of any warranty.
  12 //
  13 // Alternatively, this file may be used under the terms of Open CASCADE
  14 // commercial license or contractual agreement.
  15
  16 #include <DsgPrs_AnglePresentation.hxx>
  17
  18 #include <DsgPrs.hxx>
  19 #include <ElCLib.hxx>
  20 #include <GC_MakeCircle.hxx>
  21 #include <gce_MakePln.hxx>
  22 #include <Geom_Circle.hxx>
  23 #include <Geom_Line.hxx>
  24 #include <GeomAPI_ExtremaCurveCurve.hxx>
  25 #include <gp_Ax1.hxx>
  26 #include <gp_Circ.hxx>
  27 #include <gp_Dir.hxx>
  28 #include <gp_Lin.hxx>
  29 #include <gp_Pnt.hxx>
  30 #include <Graphic3d_ArrayOfPolylines.hxx>
  31 #include <Graphic3d_ArrayOfSegments.hxx>
  32 #include <Graphic3d_AspectLine3d.hxx>
  33 #include <Graphic3d_AspectMarker3d.hxx>
  34 #include <Graphic3d_Group.hxx>
  35 #include <Precision.hxx>
  36 #include <Prs3d_Arrow.hxx>
  37 #include <Prs3d_ArrowAspect.hxx>
  38 #include <Prs3d_DimensionAspect.hxx>
  39 #include <Prs3d_LineAspect.hxx>
  40 #include <Prs3d_Presentation.hxx>
  41 #include <Prs3d_Text.hxx>
  42 #include <TCollection_AsciiString.hxx>
  43 #include <TCollection_ExtendedString.hxx>
  44 #include <UnitsAPI.hxx>
  45
  46 #include <stdio.h>
  47 //------------------------------------------------------------------------------------------------------------------
  48 // Returns 1 if C is above of CMin; 0 if C is bitween CMin and CMax; -1 if C is Below CMax
  49 //-----------------------------------------------------------------------------------------------------------------
  50 static Standard_Integer AboveInBelowCone(const gp_Circ &CMax, const gp_Circ &CMin, const gp_Circ &C)
  51 {
  52   const Standard_Real D = CMax.Location().Distance( CMin.Location() );
  53   const Standard_Real D1 = CMax.Location().Distance( C.Location() );
  54   const Standard_Real D2 = CMin.Location().Distance( C.Location() );
  55
  56   if ( D >= D1 && D >= D2 ) return 0;
  57   if ( D < D2 && D1 < D2 ) return -1;
  58   if ( D < D1 && D2 < D1 ) return 1;
  59
  60   return 0;
  61 }
  62
  63
  64 //==========================================================================
  65 // function : DsgPrs_AnglePresentation::Add
  66 // purpose  : draws the presentation of the cone's angle;
  67 //==========================================================================
  68 void DsgPrs_AnglePresentation::Add (const Handle(Prs3d_Presentation)& aPresentation,
  69                                     const Handle(Prs3d_Drawer)& aDrawer,
  70                                     const Standard_Real /*aVal*/,
  71                                     const TCollection_ExtendedString& aText,
  72                                     const gp_Circ& aCircle,
  73                                     const gp_Pnt& aPosition,
  74                                     const gp_Pnt& Apex,
  75                                     const gp_Circ& VminCircle,
  76                                     const gp_Circ& VmaxCircle,
  77                                     const Standard_Real aArrowSize)
  78 {
  79   Handle(Prs3d_DimensionAspect) aDimensionAspect = aDrawer->DimensionAspect();
  80
  81   TCollection_ExtendedString txt(aText);
  82
  83   const Standard_Real myArrowSize = ( aArrowSize == 0.0 )? (0.1 * aCircle.Radius()) : aArrowSize;
  84
  85   aDimensionAspect->ArrowAspect()->SetLength(myArrowSize);
  86   aDrawer->ArrowAspect()->SetLength(myArrowSize);
  87
  88   Standard_Boolean IsConeTrimmed = Standard_False;
  89   gp_Circ myCircle = aCircle;
  90   if( VminCircle.Radius() > 0.01 ) {
  91     IsConeTrimmed = Standard_True;
  92     if( AboveInBelowCone( VmaxCircle, VminCircle, myCircle ) == 1 ) myCircle = VminCircle;
  93   }
  94
  95   gp_Pnt P1 = ElCLib::Value(0., myCircle);
  96   gp_Pnt P2 = ElCLib::Value(M_PI, myCircle);
  97
  98   gce_MakePln mkPln(P1, P2, Apex); // create a plane whitch defines plane for projection aPosition on it
  99
 100   gp_Vec aVector( mkPln.Value().Location(), aPosition );     //project aPosition on a plane
 101   gp_Vec Normal = mkPln.Value().Axis().Direction();
 102   Normal = (aVector * Normal) * Normal;
 103
 104   gp_Pnt aPnt = aPosition;
 105   aPnt = aPnt.Translated( -Normal );
 106
 107   gp_Pnt tmpPnt = aPnt;
 108
 109   gp_Pnt AttachmentPnt, OppositePnt;
 110   if( aPnt.Distance(P1) <  aPnt.Distance(P2) ) {
 111     AttachmentPnt = P1;
 112     OppositePnt = P2;
 113   }
 114   else {
 115     AttachmentPnt = P2;
 116     OppositePnt = P1;
 117   }
 118
 119   aPnt = AttachmentPnt ;                          // Creating of circle whitch defines a plane for a dimension arc
 120   gp_Vec Vec(AttachmentPnt, Apex);                // Dimension arc is a part of the circle
 121   Vec.Scale(2.);
 122   aPnt.Translate(Vec);
 123   GC_MakeCircle mkCirc(AttachmentPnt, OppositePnt, aPnt);
 124   gp_Circ aCircle2 = mkCirc.Value()->Circ();
 125
 126   Standard_Integer i;
 127   Standard_Real AttParam = ElCLib::Parameter(aCircle2, AttachmentPnt);  //must be equal to zero (look circle construction)
 128   Standard_Real OppParam = ElCLib::Parameter(aCircle2, OppositePnt);
 129
 130   while ( AttParam >= 2. * M_PI ) AttParam -= 2. * M_PI;
 131   while ( OppParam >= 2. * M_PI ) OppParam -= 2. * M_PI;
 132
 133   //-------------------------- Compute angle ------------------------
 134   if( txt.Length() == 0 ) {
 135     Standard_Real angle = UnitsAPI::CurrentFromLS( Abs( OppParam ),"PLANE ANGLE");
 136     char res[80];
 137     sprintf(res, "%g", angle );
 138     txt = TCollection_ExtendedString(res);
 139   }
 140   //-----------------------------------------------------------------
 141
 142   Standard_Boolean IsArrowOut = Standard_True;    //Is arrows inside or outside of the cone
 143   if( ElCLib::Parameter(aCircle2, tmpPnt) < OppParam )
 144     if( 2. * myCircle.Radius() > 4. * myArrowSize ) IsArrowOut = Standard_False;  //four times more than an arrow size
 145
 146   Standard_Real angle = OppParam - AttParam;
 147   Standard_Real param = AttParam;
 148
 149   gp_Dir aDir, aDir2;
 150   if(IsArrowOut) {
 151     aDir = gp_Dir( ( gp_Vec( ElCLib::Value( AttParam - M_PI / 12., aCircle2 ), AttachmentPnt) ) );
 152     aDir2 = gp_Dir( ( gp_Vec( ElCLib::Value( OppParam + M_PI / 12., aCircle2 ), OppositePnt) ) );
 153   }
 154   else {
 155     aDir = gp_Dir( ( gp_Vec( ElCLib::Value( AttParam + M_PI / 12., aCircle2 ), AttachmentPnt ) ) );
 156     aDir2 = gp_Dir( ( gp_Vec( ElCLib::Value( OppParam - M_PI / 12., aCircle2 ), OppositePnt ) ) );
 157   }
 158
 159   while ( angle > 2. * M_PI ) angle -= 2. * M_PI;
 160
 161   Handle(Graphic3d_ArrayOfPrimitives) aPrims = new Graphic3d_ArrayOfPolylines(12);
 162   for( i = 0; i <= 11; i++ )
 163     aPrims->AddVertex(ElCLib::Value(param + angle/11 * i, aCircle2));
 164   Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 165
 166   DsgPrs::ComputeSymbol(aPresentation, aDimensionAspect, AttachmentPnt,
 167                         AttachmentPnt, aDir, aDir, DsgPrs_AS_LASTAR);
 168   DsgPrs::ComputeSymbol(aPresentation, aDimensionAspect, OppositePnt,
 169                         OppositePnt, aDir2, aDir2, DsgPrs_AS_LASTAR);
 170
 171   param = ElCLib::Parameter(aCircle2, tmpPnt);
 172   tmpPnt = ElCLib::Value(param, aCircle2);
 173   tmpPnt = tmpPnt.Translated(gp_Vec(0, 0, -2));
 174   Prs3d_Text::Draw(aPresentation, aDimensionAspect->TextAspect(), txt, tmpPnt);   //add the TCollection_ExtendedString
 175
 176   angle = 2. * M_PI - param ;
 177   if( param > OppParam )
 178   {
 179     while ( angle > 2. * M_PI ) angle -= 2. * M_PI;
 180     aPrims = new Graphic3d_ArrayOfPolylines(12);
 181     for( i = 11; i >= 0; i-- )
 182       aPrims->AddVertex(ElCLib::Value(-angle/11 * i, aCircle2));
 183     Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 184   }
 185
 186   if( AboveInBelowCone( VmaxCircle, VminCircle, myCircle ) == 1 && !IsConeTrimmed )         //above
 187   {
 188     aPrims = new Graphic3d_ArrayOfPolylines(3);
 189     aPrims->AddVertex(AttachmentPnt);
 190     aPrims->AddVertex(Apex);
 191     aPrims->AddVertex(OppositePnt);
 192   }
 193   else
 194   {
 195     aPnt = OppositePnt ;
 196     if ( AboveInBelowCone( VmaxCircle, VminCircle, myCircle ) == 0 ) return;
 197
 198     gp_Pnt P11 = ElCLib::Value( 0., VmaxCircle );
 199     gp_Pnt P12 = ElCLib::Value( M_PI, VmaxCircle );
 200
 201     aPrims = new Graphic3d_ArrayOfSegments(4);
 202     aPrims->AddVertex(AttachmentPnt);
 203     aPrims->AddVertex(( aPnt.Distance(P1) < aPnt.Distance(P2) )? P12 : P11);
 204     aPrims->AddVertex(OppositePnt);
 205     aPrims->AddVertex(( aPnt.Distance(P1) < aPnt.Distance(P2) )? P11 : P12);
 206   }
 207   Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 208 }
 209
 210
 211 //==========================================================================
 212 // function : DsgPrs_AnglePresentation::Add
 213 // purpose  :
 214 //
 215 //==========================================================================
 216
 217 void DsgPrs_AnglePresentation::Add (const Handle(Prs3d_Presentation)& aPresentation,
 218                                     const Handle(Prs3d_Drawer)& aDrawer,
 219                                     const Standard_Real theval,
 220                                     const TCollection_ExtendedString& aText,
 221                                     const gp_Pnt& CenterPoint,
 222                                     const gp_Pnt& AttachmentPoint1,
 223                                     const gp_Pnt& AttachmentPoint2,
 224                                     const gp_Dir& dir1,
 225                                     const gp_Dir& dir2,
 226                                     const gp_Dir& axisdir,
 227                                     const gp_Pnt& OffsetPoint)
 228 {
 229   char valcar[80];
 230   sprintf(valcar,"%5.2f",theval);
 231
 232   Handle(Prs3d_DimensionAspect) LA = aDrawer->DimensionAspect();
 233   Prs3d_Root::CurrentGroup(aPresentation)->SetPrimitivesAspect(LA->LineAspect()->Aspect());
 234
 235   gp_Ax2 ax(CenterPoint,axisdir,dir1);
 236   gp_Circ cer(ax,CenterPoint.Distance(OffsetPoint));
 237   gp_Vec vec1(dir1);
 238   vec1 *= cer.Radius();
 239   gp_Vec vec2(dir2);
 240   vec2 *= cer.Radius();
 241   gp_Pnt p2 = CenterPoint.Translated(vec2);
 242
 243   Standard_Real uc1 = 0.;
 244   Standard_Real uc2 = ElCLib::Parameter(cer,p2);
 245   Standard_Real uco = ElCLib::Parameter(cer,OffsetPoint);
 246
 247   Standard_Real udeb = uc1;
 248   Standard_Real ufin = uc2;
 249
 250   if (uco > ufin) {
 251     if (Abs(theval)<M_PI) {
 252       // test if uco is in the opposite sector
 253       if (uco > udeb+M_PI && uco < ufin+M_PI) {
 254         udeb += M_PI;
 255         ufin += M_PI;
 256         uc1 = udeb;
 257         uc2 = ufin;
 258       }
 259     }
 260   }
 261
 262   if (uco > ufin) {
 263     if ((uco-uc2) < (uc1-uco+(2.*M_PI))) {
 264       ufin = uco;
 265     }
 266     else {
 267       udeb = uco - 2.*M_PI;
 268     }
 269   }
 270
 271   const Standard_Real alpha = Abs(ufin-udeb);
 272   const Standard_Integer nbp = Max (4 , Standard_Integer (50. * alpha / M_PI));
 273   const Standard_Real dteta = alpha/(nbp-1);
 274
 275   Handle(Graphic3d_ArrayOfPolylines) aPrims = new Graphic3d_ArrayOfPolylines(nbp+4,3);
 276   aPrims->AddBound(nbp);
 277   for (Standard_Integer i = 1; i<=nbp; i++)
 278     aPrims->AddVertex(ElCLib::Value(udeb+ dteta*(i-1),cer));
 279
 280   Prs3d_Text::Draw(aPresentation,LA->TextAspect(),aText,OffsetPoint);
 281
 282   Standard_Real length = LA->ArrowAspect()->Length();
 283   if (length <  Precision::Confusion()) length = 1.e-04;
 284
 285   gp_Vec vecarr;
 286   gp_Pnt ptarr;
 287   ElCLib::D1(uc1,cer,ptarr,vecarr);
 288
 289   gp_Ax1 ax1(ptarr, axisdir);
 290   gp_Dir dirarr(-vecarr);
 291
 292   //calculate angle of rotation
 293   gp_Pnt ptarr2(ptarr.XYZ() + length*dirarr.XYZ());
 294   const Standard_Real parcir = ElCLib::Parameter(cer, ptarr2);
 295   gp_Pnt ptarr3 = ElCLib::Value(parcir, cer);
 296   gp_Vec v1(ptarr,ptarr2);
 297   gp_Vec v2(ptarr,ptarr3);
 298   const Standard_Real beta = v1.Angle(v2);
 299   dirarr.Rotate(ax1, beta);
 300   Prs3d_Arrow::Draw(aPresentation,ptarr,dirarr,LA->ArrowAspect()->Angle(),length);
 301
 302   aPrims->AddBound(2);
 303   aPrims->AddVertex(AttachmentPoint1);
 304   aPrims->AddVertex(ptarr);
 305
 306   ElCLib::D1(uc2,cer,ptarr,vecarr);
 307
 308   ax1.SetLocation(ptarr);
 309   gp_Dir dirarr2(vecarr);
 310   dirarr2.Rotate(ax1,-beta);
 311   Prs3d_Arrow::Draw(aPresentation,ptarr,dirarr2,LA->ArrowAspect()->Angle(),length);
 312
 313   aPrims->AddBound(2);
 314   aPrims->AddVertex(AttachmentPoint2);
 315   aPrims->AddVertex(ptarr);
 316
 317   Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 318 }
 319
 320
 321 //==========================================================================
 322 // function : DsgPrs_AnglePresentation::Add
 323 // purpose  : Adds prezentation of angle between two faces
 324 //==========================================================================
 325
 326 void DsgPrs_AnglePresentation::Add (const Handle(Prs3d_Presentation)& aPresentation,
 327                                     const Handle(Prs3d_Drawer)& aDrawer,
 328                                     const Standard_Real theval,
 329                                     const TCollection_ExtendedString& aText,
 330                                     const gp_Pnt& CenterPoint,
 331                                     const gp_Pnt& AttachmentPoint1,
 332                                     const gp_Pnt& AttachmentPoint2,
 333                                     const gp_Dir& dir1,
 334                                     const gp_Dir& dir2,
 335                                     const gp_Dir& axisdir,
 336                                     const Standard_Boolean isPlane,
 337                                     const gp_Ax1& AxisOfSurf,
 338                                     const gp_Pnt& OffsetPoint,
 339                                     const DsgPrs_ArrowSide ArrowPrs )
 340 {
 341   char valcar[80];
 342   sprintf(valcar,"%5.2f",theval);
 343
 344   Handle( Prs3d_DimensionAspect ) LA = aDrawer->DimensionAspect();
 345   Prs3d_Root::CurrentGroup( aPresentation )->SetPrimitivesAspect( LA->LineAspect()->Aspect() );
 346
 347   gp_Circ AngleCirc, AttachCirc;
 348   Standard_Real FirstParAngleCirc, LastParAngleCirc, FirstParAttachCirc, LastParAttachCirc;
 349   gp_Pnt EndOfArrow1, EndOfArrow2, ProjAttachPoint2;
 350   gp_Dir DirOfArrow1, DirOfArrow2;
 351   DsgPrs::ComputeFacesAnglePresentation( LA->ArrowAspect()->Length(),
 352                                          theval,
 353                                          CenterPoint,
 354                                          AttachmentPoint1,
 355                                          AttachmentPoint2,
 356                                          dir1,
 357                                          dir2,
 358                                          axisdir,
 359                                          isPlane,
 360                                          AxisOfSurf,
 361                                          OffsetPoint,
 362                                          AngleCirc,
 363                                          FirstParAngleCirc,
 364                                          LastParAngleCirc,
 365                                          EndOfArrow1,
 366                                          EndOfArrow2,
 367                                          DirOfArrow1,
 368                                          DirOfArrow2,
 369                                          ProjAttachPoint2,
 370                                          AttachCirc,
 371                                          FirstParAttachCirc,
 372                                          LastParAttachCirc );
 373
 374   // Creating the angle's arc or line if null angle
 375   Handle(Graphic3d_ArrayOfPrimitives) aPrims;
 376   if (theval > Precision::Angular() && Abs( M_PI-theval ) > Precision::Angular())
 377   {
 378     const Standard_Real Alpha  = Abs( LastParAngleCirc - FirstParAngleCirc );
 379     const Standard_Integer NodeNumber = Max (4 , Standard_Integer (50. * Alpha / M_PI));
 380     const Standard_Real delta = Alpha / (Standard_Real)( NodeNumber - 1 );
 381
 382     aPrims = new Graphic3d_ArrayOfPolylines(NodeNumber+4,3);
 383     aPrims->AddBound(NodeNumber);
 384     for (Standard_Integer i = 0; i < NodeNumber; i++, FirstParAngleCirc += delta)
 385       aPrims->AddVertex(ElCLib::Value( FirstParAngleCirc, AngleCirc ));
 386
 387     Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 388     aPrims = new Graphic3d_ArrayOfSegments(4);
 389   }
 390   else // null angle
 391   {
 392     aPrims = new Graphic3d_ArrayOfSegments(6);
 393     aPrims->AddVertex(OffsetPoint);
 394     aPrims->AddVertex(EndOfArrow1);
 395   }
 396
 397   // Add presentation of arrows
 398   DsgPrs::ComputeSymbol( aPresentation, LA, EndOfArrow1, EndOfArrow2, DirOfArrow1, DirOfArrow2, ArrowPrs );
 399
 400   // Drawing the text
 401   Prs3d_Text::Draw( aPresentation, LA->TextAspect(), aText, OffsetPoint );
 402
 403   // Line from AttachmentPoint1 to end of Arrow1
 404   aPrims->AddVertex(AttachmentPoint1);
 405   aPrims->AddVertex(EndOfArrow1);
 406   // Line from "projection" of AttachmentPoint2 to end of Arrow2
 407   aPrims->AddVertex(ProjAttachPoint2);
 408   aPrims->AddVertex(EndOfArrow2);
 409
 410   Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 411
 412   // Line or arc from AttachmentPoint2 to its "projection"
 413   if (AttachmentPoint2.Distance( ProjAttachPoint2 ) > Precision::Confusion())
 414   {
 415     if (isPlane)
 416         {
 417           // Creating the line from AttachmentPoint2 to its projection
 418       aPrims = new Graphic3d_ArrayOfSegments(2);
 419       aPrims->AddVertex(AttachmentPoint2);
 420       aPrims->AddVertex(ProjAttachPoint2);
 421         }
 422     else
 423         {
 424           // Creating the arc from AttachmentPoint2 to its projection
 425       const Standard_Real Alpha = Abs( LastParAttachCirc - FirstParAttachCirc );
 426           const Standard_Integer NodeNumber = Max (4 , Standard_Integer (50. * Alpha / M_PI));
 427       const Standard_Real delta = Alpha / (Standard_Real)( NodeNumber - 1 );
 428
 429       aPrims = new Graphic3d_ArrayOfPolylines(NodeNumber);
 430           for (Standard_Integer i = 0; i < NodeNumber; i++, FirstParAttachCirc += delta)
 431         aPrims->AddVertex(ElCLib::Value( FirstParAttachCirc, AttachCirc ));
 432         }
 433     Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 434   }
 435 }
 436
 437
 438 //==========================================================================
 439 // function : DsgPrs_AnglePresentation::Add
 440 // purpose  :
 441 //
 442 //==========================================================================
 443
 444 void DsgPrs_AnglePresentation::Add (const Handle(Prs3d_Presentation)& aPresentation,
 445                                     const Handle(Prs3d_Drawer)& aDrawer,
 446                                     const Standard_Real theval,
 447                                     const TCollection_ExtendedString& aText,
 448                                     const gp_Pnt& CenterPoint,
 449                                     const gp_Pnt& AttachmentPoint1,
 450                                     const gp_Pnt& AttachmentPoint2,
 451                                     const gp_Dir& dir1,
 452                                     const gp_Dir& dir2,
 453                                     const gp_Pnt& OffsetPoint)
 454 {
 455   char valcar[80];
 456   sprintf(valcar,"%5.2f",theval);
 457
 458   Handle(Prs3d_DimensionAspect) LA = aDrawer->DimensionAspect();
 459   Prs3d_Root::CurrentGroup(aPresentation)->SetPrimitivesAspect(LA->LineAspect()->Aspect());
 460
 461   gp_Dir Norm;
 462   if (!dir1.IsParallel(dir2, Precision::Angular())) {
 463     Norm = dir1.Crossed(dir2);
 464   }
 465   else {
 466     gp_Dir dir2B = gp_Dir(gp_Vec(CenterPoint, OffsetPoint));
 467     Norm = dir1.Crossed(dir2B);
 468   }
 469
 470   if (Abs(theval) > M_PI) Norm.Reverse();
 471
 472   gp_Ax2 ax(CenterPoint,Norm,dir1);
 473   gp_Circ cer(ax,CenterPoint.Distance(OffsetPoint));
 474   gp_Vec vec1(dir1);
 475   vec1 *= cer.Radius();
 476   gp_Vec vec2(dir2);
 477   vec2 *= cer.Radius();
 478   gp_Pnt p2 = CenterPoint.Translated(vec2);
 479
 480   Standard_Real uc1 = 0.;
 481   Standard_Real uc2 = ElCLib::Parameter(cer,p2);
 482   Standard_Real uco = ElCLib::Parameter(cer,OffsetPoint);
 483
 484   Standard_Real udeb = uc1;
 485   Standard_Real ufin = uc2;
 486
 487   if (uco > ufin) {
 488     if (Abs(theval)<M_PI) {
 489       // test if uco is in the opposite sector
 490       if (uco > udeb+M_PI && uco < ufin+M_PI) {
 491         udeb += M_PI;
 492         ufin += M_PI;
 493         uc1 = udeb;
 494         uc2 = ufin;
 495       }
 496     }
 497   }
 498
 499   if (uco > ufin) {
 500     if ((uco-uc2) < (uc1-uco+(2.*M_PI))) {
 501       ufin = uco;
 502     }
 503     else {
 504       udeb = uco - 2.*M_PI;
 505     }
 506   }
 507
 508   const Standard_Real alpha = Abs(ufin-udeb);
 509   const Standard_Integer nbp = Max (4 , Standard_Integer (50. * alpha / M_PI));
 510   const Standard_Real dteta = alpha/(nbp-1);
 511
 512   Handle(Graphic3d_ArrayOfPolylines) aPrims = new Graphic3d_ArrayOfPolylines(nbp+4,3);
 513   aPrims->AddBound(nbp);
 514   for (Standard_Integer i = 1; i<=nbp; i++)
 515     aPrims->AddVertex(ElCLib::Value(udeb+ dteta*(i-1),cer));
 516
 517   Prs3d_Text::Draw(aPresentation,LA->TextAspect(),aText,OffsetPoint);
 518
 519   Standard_Real length = LA->ArrowAspect()->Length();
 520   if (length <  Precision::Confusion()) length = 1.e-04;
 521
 522   gp_Vec vecarr;
 523   gp_Pnt ptarr;
 524   ElCLib::D1(uc1,cer,ptarr,vecarr);
 525
 526   gp_Ax1 ax1(ptarr, Norm);
 527   gp_Dir dirarr(-vecarr);
 528   //calculate the angle of rotation
 529   gp_Pnt ptarr2(ptarr.XYZ() + length*dirarr.XYZ());
 530   const Standard_Real parcir = ElCLib::Parameter(cer, ptarr2);
 531   gp_Pnt ptarr3 = ElCLib::Value(parcir, cer);
 532   gp_Vec v1(ptarr,ptarr2);
 533   gp_Vec v2(ptarr,ptarr3);
 534   const Standard_Real beta = v1.Angle(v2);
 535   dirarr.Rotate(ax1, beta);
 536   Prs3d_Arrow::Draw(aPresentation,ptarr,dirarr,LA->ArrowAspect()->Angle(),length);
 537
 538   aPrims->AddBound(2);
 539   aPrims->AddVertex(AttachmentPoint1);
 540   aPrims->AddVertex(ptarr);
 541
 542   ElCLib::D1(uc2,cer,ptarr,vecarr);
 543
 544   ax1.SetLocation(ptarr);
 545   gp_Dir dirarr2(vecarr);
 546   dirarr2.Rotate(ax1,  - beta);
 547   Prs3d_Arrow::Draw(aPresentation,ptarr,dirarr2,LA->ArrowAspect()->Angle(),length);
 548
 549   aPrims->AddBound(2);
 550   aPrims->AddVertex(AttachmentPoint2);
 551   aPrims->AddVertex(ptarr);
 552
 553   Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 554 }
 555
 556
 557 //==========================================================================
 558 // function : DsgPrs_AnglePresentation::Add
 559 // purpose  : It is possible to choose the symbol of extremities of the face (arrow, point...)
 560 //==========================================================================
 561
 562 void DsgPrs_AnglePresentation::Add (const Handle(Prs3d_Presentation)& aPresentation,
 563                                     const Handle(Prs3d_Drawer)& aDrawer,
 564                                     const Standard_Real theval,
 565                                     const TCollection_ExtendedString& aText,
 566                                     const gp_Pnt& CenterPoint,
 567                                     const gp_Pnt& AttachmentPoint1,
 568                                     const gp_Pnt& AttachmentPoint2,
 569                                     const gp_Dir& dir1,
 570                                     const gp_Dir& dir2,
 571                                     const gp_Pnt& OffsetPoint,
 572                                     const DsgPrs_ArrowSide ArrowPrs)
 573 {
 574   char valcar[80];
 575   sprintf(valcar,"%5.2f",theval);
 576
 577   Handle(Prs3d_DimensionAspect) LA = aDrawer->DimensionAspect();
 578   Prs3d_Root::CurrentGroup(aPresentation)->SetPrimitivesAspect(LA->LineAspect()->Aspect());
 579
 580   gp_Dir Norm = dir1.Crossed(dir2);
 581   if (Abs(theval) > M_PI) Norm.Reverse();
 582
 583   gp_Ax2 ax(CenterPoint,Norm,dir1);
 584   gp_Circ cer(ax,CenterPoint.Distance(OffsetPoint));
 585   gp_Vec vec1(dir1);
 586   vec1 *= cer.Radius();
 587   gp_Vec vec2(dir2);
 588   vec2 *= cer.Radius();
 589   gp_Pnt p2 = CenterPoint.Translated(vec2);
 590
 591   Standard_Real uc1 = 0.;
 592   Standard_Real uc2 = ElCLib::Parameter(cer,p2);
 593   Standard_Real uco = ElCLib::Parameter(cer,OffsetPoint);
 594
 595   Standard_Real udeb = uc1;
 596   Standard_Real ufin = uc2;
 597
 598   if (uco > ufin) {
 599     if (Abs(theval)<M_PI) {
 600       // test if uco is in the opposite sector
 601       if (uco > udeb+M_PI && uco < ufin+M_PI) {
 602         udeb += M_PI;
 603         ufin += M_PI;
 604         uc1 = udeb;
 605         uc2 = ufin;
 606       }
 607     }
 608   }
 609
 610   if (uco > ufin) {
 611     if ((uco-uc2) < (uc1-uco+(2.*M_PI))) {
 612       ufin = uco;
 613     }
 614     else {
 615       udeb = uco - 2.*M_PI;
 616     }
 617   }
 618
 619   const Standard_Real alpha = Abs(ufin-udeb);
 620   const Standard_Integer nbp = Max (4 , Standard_Integer (50. * alpha / M_PI));
 621   const Standard_Real dteta = alpha/(nbp-1);
 622
 623   Handle(Graphic3d_ArrayOfPolylines) aPrims = new Graphic3d_ArrayOfPolylines(nbp+4,3);
 624   aPrims->AddBound(nbp);
 625   for (Standard_Integer i = 1; i<=nbp; i++)
 626     aPrims->AddVertex(ElCLib::Value(udeb+ dteta*(i-1),cer));
 627
 628   Prs3d_Text::Draw(aPresentation,LA->TextAspect(),aText,OffsetPoint);
 629
 630   Standard_Real length = LA->ArrowAspect()->Length();
 631   if (length <  Precision::Confusion()) length = 1.e-04;
 632
 633   // Lines of recall
 634   gp_Vec vecarr;
 635   gp_Pnt ptarr;
 636   ElCLib::D1(uc1,cer,ptarr,vecarr);
 637
 638   gp_Ax1 ax1(ptarr, Norm);
 639   gp_Dir dirarr(-vecarr);
 640   //calculate angle of rotation
 641   gp_Pnt ptarr2(ptarr.XYZ() + length*dirarr.XYZ());
 642   const Standard_Real parcir = ElCLib::Parameter(cer, ptarr2);
 643   gp_Pnt ptarr3 = ElCLib::Value(parcir, cer);
 644   gp_Vec v1(ptarr,ptarr2 );
 645   gp_Vec v2(ptarr, ptarr3);
 646   const Standard_Real beta = v1.Angle(v2);
 647   dirarr.Rotate(ax1, beta);
 648
 649   aPrims->AddBound(2);
 650   aPrims->AddVertex(AttachmentPoint1);
 651   aPrims->AddVertex(ptarr);
 652
 653   gp_Vec vecarr1;
 654   gp_Pnt ptarr1;
 655   ElCLib::D1(uc2,cer,ptarr1,vecarr1);
 656   ax1.SetLocation(ptarr1);
 657   gp_Dir dirarr2(vecarr1);
 658   dirarr2.Rotate(ax1,  - beta);
 659
 660   aPrims->AddBound(2);
 661   aPrims->AddVertex(AttachmentPoint2);
 662   aPrims->AddVertex(ptarr1);
 663
 664   Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 665
 666   // One traces the arrows
 667   DsgPrs::ComputeSymbol(aPresentation,LA,ptarr,ptarr1,dirarr,dirarr2,ArrowPrs);
 668 }
 669
 670
 671 //==========================================================================
 672 // function : DsgPrs_AnglePresentation::Add
 673 // purpose  :
 674 //
 675 //==========================================================================
 676
 677 void DsgPrs_AnglePresentation::Add (const Handle(Prs3d_Presentation)& aPresentation,
 678                                     const Handle(Prs3d_Drawer)& aDrawer,
 679                                     const Standard_Real theval,
 680                                     const gp_Pnt& CenterPoint,
 681                                     const gp_Pnt& AttachmentPoint1,
 682                                     const gp_Pnt& AttachmentPoint2,
 683                                     const gp_Dir& dir1,
 684                                     const gp_Dir& dir2,
 685                                     const gp_Pnt& OffsetPoint)
 686 {
 687   char valcar[80];
 688   sprintf(valcar,"%5.2f",theval);
 689
 690   TCollection_AsciiString valas(valcar);
 691   TCollection_ExtendedString aText(valas);
 692
 693   Handle(Prs3d_DimensionAspect) LA = aDrawer->DimensionAspect();
 694   Prs3d_Root::CurrentGroup(aPresentation)->SetPrimitivesAspect(LA->LineAspect()->Aspect());
 695
 696   gp_Dir Norm = dir1.Crossed(dir2);
 697   if (Abs(theval) > M_PI) Norm.Reverse();
 698
 699   gp_Ax2 ax(CenterPoint,Norm,dir1);
 700   gp_Circ cer(ax,CenterPoint.Distance(OffsetPoint));
 701   gp_Vec vec1(dir1);
 702   vec1 *= cer.Radius();
 703   gp_Vec vec2(dir2);
 704   vec2 *= cer.Radius();
 705   gp_Pnt p2 = CenterPoint.Translated(vec2);
 706
 707   Standard_Real uc1 = 0.;
 708   Standard_Real uc2 = ElCLib::Parameter(cer,p2);
 709   Standard_Real uco = ElCLib::Parameter(cer,OffsetPoint);
 710
 711   Standard_Real udeb = uc1;
 712   Standard_Real ufin = uc2;
 713
 714   if (uco > ufin) {
 715     if (Abs(theval)<M_PI) {
 716       // test if uco is in the opposite sector
 717       if (uco > udeb+M_PI && uco < ufin+M_PI) {
 718         udeb += M_PI;
 719         ufin += M_PI;
 720         uc1 = udeb;
 721         uc2 = ufin;
 722       }
 723     }
 724   }
 725
 726   if (uco > ufin) {
 727     if ((uco-uc2) < (uc1-uco+(2.*M_PI))) {
 728       ufin = uco;
 729     }
 730     else {
 731       udeb = uco - 2.*M_PI;
 732     }
 733   }
 734
 735   const Standard_Real alpha = Abs(ufin-udeb);
 736   const Standard_Integer nbp = Max (4 , Standard_Integer (50. * alpha / M_PI));
 737   const Standard_Real dteta = alpha/(nbp-1);
 738
 739   Handle(Graphic3d_ArrayOfPolylines) aPrims = new Graphic3d_ArrayOfPolylines(nbp+4,3);
 740   aPrims->AddBound(nbp);
 741   for (Standard_Integer i = 1; i<=nbp; i++)
 742     aPrims->AddVertex(ElCLib::Value(udeb+ dteta*(i-1),cer));
 743
 744   Prs3d_Text::Draw(aPresentation,LA->TextAspect(),aText,OffsetPoint);
 745
 746   Standard_Real length = LA->ArrowAspect()->Length();
 747   if (length <  Precision::Confusion()) length = 1.e-04;
 748
 749   gp_Vec vecarr;
 750   gp_Pnt ptarr;
 751   ElCLib::D1(uc1,cer,ptarr,vecarr);
 752
 753   gp_Ax1 ax1(ptarr, Norm);
 754   gp_Dir dirarr(-vecarr);
 755   //calculate the angle of rotation
 756   gp_Pnt ptarr2(ptarr.XYZ() + length*dirarr.XYZ());
 757   const Standard_Real parcir = ElCLib::Parameter(cer, ptarr2);
 758   gp_Pnt ptarr3 = ElCLib::Value(parcir, cer);
 759   gp_Vec v1(ptarr,ptarr2 );
 760   gp_Vec v2(ptarr, ptarr3);
 761   const Standard_Real beta = v1.Angle(v2);
 762   dirarr.Rotate(ax1, beta);
 763
 764   Prs3d_Arrow::Draw(aPresentation,ptarr,dirarr,LA->ArrowAspect()->Angle(),length);
 765
 766   aPrims->AddBound(2);
 767   aPrims->AddVertex(AttachmentPoint1);
 768   aPrims->AddVertex(ptarr);
 769
 770   ElCLib::D1(uc2,cer,ptarr,vecarr);
 771   ax1.SetLocation(ptarr);
 772   gp_Dir dirarr2(vecarr);
 773   dirarr2.Rotate(ax1, -beta);
 774
 775   Prs3d_Arrow::Draw(aPresentation,ptarr,dirarr2,LA->ArrowAspect()->Angle(),length);
 776
 777   aPrims->AddBound(2);
 778   aPrims->AddVertex(AttachmentPoint2);
 779   aPrims->AddVertex(ptarr);
 780
 781   Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 782 }
 783
 784 void DsgPrs_AnglePresentation::Add (const Handle(Prs3d_Presentation)& aPresentation,
 785                                     const Handle(Prs3d_Drawer)& aDrawer,
 786                                     const Standard_Real theval,
 787                                     const gp_Pnt& CenterPoint,
 788                                     const gp_Pnt& AttachmentPoint1,
 789                                     const gp_Ax1& theAxe,
 790                                     const DsgPrs_ArrowSide ArrowSide)
 791 {
 792   Handle(Prs3d_DimensionAspect) LA = aDrawer->DimensionAspect();
 793   Prs3d_Root::CurrentGroup(aPresentation)->SetPrimitivesAspect(LA->LineAspect()->Aspect());
 794
 795   gp_Dir dir1(gp_Vec(CenterPoint, AttachmentPoint1));
 796   gp_Ax2 ax(CenterPoint,theAxe.Direction(),dir1);
 797   gp_Circ cer(ax,CenterPoint.Distance(AttachmentPoint1));
 798
 799   const Standard_Integer nbp = Max (4 , Standard_Integer (50. * theval / M_PI));
 800   const Standard_Real dteta = theval/(nbp-1);
 801
 802   Handle(Graphic3d_ArrayOfPolylines) aPrims = new Graphic3d_ArrayOfPolylines(nbp);
 803   for (Standard_Integer i = 1; i<=nbp; i++)
 804     aPrims->AddVertex(ElCLib::Value(dteta*(i-1),cer));
 805   Prs3d_Root::CurrentGroup(aPresentation)->AddPrimitiveArray(aPrims);
 806
 807   Standard_Real uc1 = 0.;
 808   Standard_Real uc2 = ElCLib::Parameter(cer,AttachmentPoint1.Rotated(theAxe,theval));
 809
 810   Standard_Real length = LA->ArrowAspect()->Length();
 811   if (length < Precision::Confusion()) length = 1.e-04;
 812
 813   gp_Vec vecarr;
 814   gp_Pnt ptarr;
 815   switch(ArrowSide)
 816   {
 817     case DsgPrs_AS_FIRSTAR:
 818     {
 819       ElCLib::D1(uc1,cer,ptarr,vecarr);
 820       Prs3d_Arrow::Draw(aPresentation,ptarr,gp_Dir(-vecarr),LA->ArrowAspect()->Angle(),length);
 821       break;
 822     }
 823     case DsgPrs_AS_LASTAR:
 824     {
 825       ElCLib::D1(uc2,cer,ptarr,vecarr);
 826       Prs3d_Arrow::Draw(aPresentation,ptarr,gp_Dir(vecarr),LA->ArrowAspect()->Angle(),length);
 827       break;
 828     }
 829     case DsgPrs_AS_BOTHAR:
 830     {
 831       ElCLib::D1(uc1,cer,ptarr,vecarr);
 832       Prs3d_Arrow::Draw(aPresentation,ptarr,gp_Dir(-vecarr),LA->ArrowAspect()->Angle(),length);
 833       ElCLib::D1(uc2,cer,ptarr,vecarr);
 834       Prs3d_Arrow::Draw(aPresentation,ptarr,gp_Dir(vecarr),LA->ArrowAspect()->Angle(),length);
 835       break;
 836     }
 837     default: break;
 838   }
 839 }