doc/salome/gui/GEOM/input/tui_angle.doc

   1 /*!
   2
   3 \page tui_angle_page Angle
   4
   5 \code
   6 import salome
   7 salome.salome_init()
   8
   9 import math
  10 import geompy
  11 geompy.init_geom(salome.myStudy)
  12
  13 OX  = geompy.MakeVectorDXDYDZ(10, 0,0)
  14 OXY = geompy.MakeVectorDXDYDZ(10,10,0)
  15
  16 # in one plane
  17 Angle = geompy.GetAngle(OX, OXY)
  18
  19 print "\nAngle between OX and OXY = ", Angle
  20 if math.fabs(Angle - 45.0) > 1e-05:
  21     print "  Error: returned angle is", Angle, "while must be 45.0"
  22     pass
  23
  24 Angle = geompy.GetAngleRadians(OX, OXY)
  25
  26 print "\nAngle between OX and OXY in radians = ", Angle
  27 if math.fabs(Angle - math.pi/4) > 1e-05:
  28     print "  Error: returned angle is", Angle, "while must be pi/4"
  29     pass
  30
  31 # not in one plane
  32 OXY_shift = geompy.MakeTranslation(OXY,10,-10,20)
  33 Angle = geompy.GetAngle(OX, OXY_shift)
  34
  35 print "Angle between OX and OXY_shift = ", Angle
  36 if math.fabs(Angle - 45.0) > 1e-05:
  37     print "  Error: returned angle is", Angle, "while must be 45.0"
  38     pass
  39
  40 # not linear
  41 pnt1 = geompy.MakeVertex(0, 0, 0)
  42 pnt2 = geompy.MakeVertex(10, 0, 0)
  43 pnt3 = geompy.MakeVertex(20, 10, 0)
  44 arc  = geompy.MakeArc(pnt1, pnt2, pnt3)
  45 Angle = geompy.GetAngle(OX, arc)
  46
  47 if (math.fabs(Angle + 1.0) > 1e-6 or geompy.MeasuOp.IsDone()):
  48     print "Error. Angle must not be computed on curvilinear edges"
  49     pass
  50
  51 \endcode
  52
  53 */