GEOM_TestMeasures.py

GEOM_TestMeasures.py

def TestMeasureOperations (geompy, math):
 
  p0   = geompy.MakeVertex(0 ,  0,  0)
  p137 = geompy.MakeVertex(10, 30, 70)
 
  box = geompy.MakeBoxTwoPnt(p0, p137)
 
  p678 = geompy.MakeVertex(60, 70, 80)
  p789 = geompy.MakeVertex(70, 80, 90)
 
  vz = geompy.MakeVectorDXDYDZ(0, 0, 1)
 
  cube = geompy.MakeBoxTwoPnt(p678, p789)
 
  cylinder = geompy.MakeCylinder(p0, vz, 5, 70)
 
  

 
  Coords = geompy.PointCoordinates(p137)
  if Coords[0] != 10 or Coords[1] != 30 or Coords[2] != 70:
    print("Coordinates of p137 must be (10, 30, 70), but returned (", Coords[0], ", ", Coords[1], ", ", Coords[2], ")")
 
  

 
  (IsValid, err) = geompy.CheckShape(box, 0, 2)
  if IsValid == 0:
    geompy.PrintShapeError(box, err)
    raise RuntimeError("Invalid box created")
  else:
    print("\nBox is valid")
 
  

 
  selfIntersected = geompy.MakeCompound([box, cylinder])
  if geompy.CheckSelfIntersections(selfIntersected):
    raise RuntimeError("Existing self-intersection is not detected")
 
  

 
  if salome_version.getXVersion() > "0x70600":
    if geompy.CheckSelfIntersectionsFast(selfIntersected):
      raise RuntimeError("Existing self-intersection is not detected")
 
  

 
  if not geompy.FastIntersect(box, cylinder)[0]:
    raise RuntimeError("Existing intersection is not detected")
 
  

 
  Descr = geompy.WhatIs(box)
  print("\nBox 10x30x70 description:")
  print(Descr)
 
  
 
  NbSolids = geompy.NbShapes(box, geompy.ShapeType["SOLID"])
  print("\nBox 10x30x70 quantity of solids:", NbSolids)
 
  
 
  BoxInfo = geompy.ShapeInfo(box)
  print("\nBox 10x30x70 shapes:")
  print(BoxInfo)
 
  

 
  Props = geompy.BasicProperties(box)
  print("\nBox 10x30x70 Basic Properties:")
  print(" Wires length: ", Props[0])
  print(" Surface area: ", Props[1])
  print(" Volume      : ", Props[2])
 
  dl = math.sqrt((Props[0] -   880)*(Props[0] -   880))
  da = math.sqrt((Props[1] -  6200)*(Props[1] -  6200))
  dv = math.sqrt((Props[2] - 21000)*(Props[2] - 21000))
  #print "|Props[0] - 880| = ", dl
 
  if dl > 1e-7 or da > 1e-7 or dv > 1e-7:
    print("While must be:")
    print(" Wires length: ", 880)
    print(" Surface area: ", 6200)
    print(" Volume      : ", 21000)
 
  

 
  BB = geompy.BoundingBox(box)
  print("\nBounding Box of box 10x30x70:")
  print(" Xmin = ", BB[0], ", Xmax = ", BB[1])
  print(" Ymin = ", BB[2], ", Ymax = ", BB[3])
  print(" Zmin = ", BB[4], ", Zmax = ", BB[5])
  BB = geompy.MakeBoundingBox(box)
  geompy.addToStudy(BB, "BoundingBox")
 
  

 
  In = geompy.Inertia(box)
  print("\nInertia matrix of box 10x30x70:")
  print(" (", In[0], ", ", In[1], ", ", In[2], ")")
  print(" (", In[3], ", ", In[4], ", ", In[5], ")")
  print(" (", In[6], ", ", In[7], ", ", In[8], ")")
  print("Main moments of inertia of box 10x30x70:")
  print(" Ix = ", In[9], ", Iy = ", In[10], ", Iz = ", In[11])
 
  

 
  Toler = geompy.Tolerance(box)
  print("\nBox 10x30x70 tolerance:")
  print(" Face min. tolerance: ", Toler[0])
  print(" Face max. tolerance: ", Toler[1])
  print(" Edge min. tolerance: ", Toler[2])
  print(" Edge max. tolerance: ", Toler[3])
  print(" Vertex min. tolerance: ", Toler[4])
  print(" Vertex max. tolerance: ", Toler[5])
 
  

 
  pcdg = geompy.MakeCDG(box)
  if pcdg is None:
    raise RuntimeError("MakeCDG(box) failed")
  else:
    print("\nCentre of gravity of box has been successfully obtained:")
    Coords = geompy.PointCoordinates(pcdg)
    print("(", Coords[0], ", ", Coords[1], ", ", Coords[2], ")")
    if Coords[0] != 5 or Coords[1] != 15 or Coords[2] != 35:
      print("But must be (5, 15, 35)")
 
  
 
  faces = geompy.SubShapeAllSortedCentres(box, geompy.ShapeType["FACE"])
  face0 = faces[0]

  vnorm = geompy.GetNormal(face0)
  if vnorm is None:
    raise RuntimeError("GetNormal(face0) failed")
  else:
    geompy.addToStudy(face0, "Face0")
    geompy.addToStudy(vnorm, "Normale to Face0")
    print("\nNormale of face has been successfully obtained:")
    #Coords = geompy.PointCoordinates(pcdg)
    #print "(", Coords[0], ", ", Coords[1], ", ", Coords[2], ")"
    #if Coords[0] != 5 or Coords[1] != 15 or Coords[2] != 35:
    #  print "But must be (5, 15, 35)"
 
  

 
  MinDist = geompy.MinDistance(box, cube)
 
  #print "\nMinimal distance between Box and Cube = ", MinDist[0]
  #print "It is reached at points:"
  #print " On Box  (", MinDist[1], ", ", MinDist[2], ", ", MinDist[3], ")"
  #print " On Cube (", MinDist[4], ", ", MinDist[5], ", ", MinDist[6], ")"
 
  print("\nMinimal distance between Box and Cube = ", MinDist)
 
  MinDistComps = geompy.MinDistanceComponents(box, cube)
  print("\nMinimal distance between Box and Cube = ", MinDistComps[0])
  print("Its components are  (", MinDistComps[1], ", ", MinDistComps[2], ", ", MinDistComps[3], ")")
 
  # Get all closest points
  [nbSols, listCoords] = geompy.ClosestPoints(box, cube)
  for i in range(nbSols):
    v1 = geompy.MakeVertex(listCoords[i*6 + 0], listCoords[i*6 + 1], listCoords[i*6 + 2])
    v2 = geompy.MakeVertex(listCoords[i*6 + 3], listCoords[i*6 + 4], listCoords[i*6 + 5])
 
    geompy.addToStudy(v1, 'MinDist_%d_on_Box'%(i+1))
    geompy.addToStudy(v2, 'MinDist_%d_on_Cube'%(i+1))
    pass
 
  
 
  OX  = geompy.MakeVectorDXDYDZ(10, 0,0)
  OXY = geompy.MakeVectorDXDYDZ(10,10,0)
 
  # in one plane

  Angle = geompy.GetAngle(OX, OXY)
 
  print("\nAngle between OX and OXY = ", Angle)
  if math.fabs(Angle - 45.0) > 1e-05:
    print("  Error: returned angle is", Angle, "while must be 45.0")

 
  Angle = geompy.GetAngleRadians(OX, OXY)
 
  print("\nAngle between OX and OXY in radians = ", Angle)
  if math.fabs(Angle - math.pi/4) > 1e-05:
    print("  Error: returned angle is", Angle, "while must be pi/4")
    pass