# -*- coding: iso-8859-1 -*-
-# Copyright (C) 2007-2014 CEA/DEN, EDF R&D, OPEN CASCADE
+# Copyright (C) 2007-2020 CEA/DEN, EDF R&D, OPEN CASCADE
#
# Copyright (C) 2003-2007 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
# CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
# See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
#
+import salome_version
+
def TestMeasureOperations (geompy, math):
p0 = geompy.MakeVertex(0 , 0, 0)
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)
+
####### PointCoordinates #######
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], ")"
+ print("Coordinates of p137 must be (10, 30, 70), but returned (", Coords[0], ", ", Coords[1], ", ", Coords[2], ")")
####### CheckShape #######
- IsValid = geompy.CheckShape(box)
+ (IsValid, err) = geompy.CheckShape(box, 0, 2)
if IsValid == 0:
- raise RuntimeError, "Invalid box created"
+ geompy.PrintShapeError(box, err)
+ raise RuntimeError("Invalid box created")
else:
- print "\nBox is valid"
+ print("\nBox is valid")
####### Detect Self-intersections #######
- [Face_1,Face_2] = geompy.SubShapes(box, [33, 23])
- Translation_1 = geompy.MakeTranslation(Face_1, 5, -15, -40)
- Compound_1 = geompy.MakeCompound([Face_2, Translation_1])
- if geompy.CheckSelfIntersections(Compound_1) == True:
- raise RuntimeError, "Existing self-intersection is not detected"
+ selfIntersected = geompy.MakeCompound([box, cylinder])
+ if geompy.CheckSelfIntersections(selfIntersected):
+ raise RuntimeError("Existing self-intersection is not detected")
+
+ ####### Detect Self-intersections fast #######
+
+ if salome_version.getXVersion() > "0x70600":
+ if geompy.CheckSelfIntersectionsFast(selfIntersected):
+ raise RuntimeError("Existing self-intersection is not detected")
+
+ ####### Fast intersection #######
+
+ if not geompy.FastIntersect(box, cylinder)[0]:
+ raise RuntimeError("Existing intersection is not detected")
####### WhatIs #######
Descr = geompy.WhatIs(box)
- print "\nBox 10x30x70 description:"
- print Descr
+ print("\nBox 10x30x70 description:")
+ print(Descr)
####### NbShapes #######
NbSolids = geompy.NbShapes(box, geompy.ShapeType["SOLID"])
- print "\nBox 10x30x70 quantity of solids:", NbSolids
+ print("\nBox 10x30x70 quantity of solids:", NbSolids)
####### ShapeInfo #######
BoxInfo = geompy.ShapeInfo(box)
- print "\nBox 10x30x70 shapes:"
- print BoxInfo
+ print("\nBox 10x30x70 shapes:")
+ print(BoxInfo)
####### BasicProperties #######
Props = geompy.BasicProperties(box)
- print "\nBox 10x30x70 Basic Properties:"
- print " Wires length: ", Props[0]
- print " Surface area: ", Props[1]
- print " Volume : ", Props[2]
+ 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))
#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
+ print("While must be:")
+ print(" Wires length: ", 880)
+ print(" Surface area: ", 6200)
+ print(" Volume : ", 21000)
####### BoundingBox #######
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]
+ 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")
####### Inertia #######
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]
+ 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])
####### Tolerance #######
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]
+ 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])
####### MakeCDG #######
pcdg = geompy.MakeCDG(box)
if pcdg is None:
- raise RuntimeError, "MakeCDG(box) failed"
+ raise RuntimeError("MakeCDG(box) failed")
else:
- print "\nCentre of gravity of box has been successfully obtained:"
+ print("\nCentre of gravity of box has been successfully obtained:")
Coords = geompy.PointCoordinates(pcdg)
- print "(", Coords[0], ", ", Coords[1], ", ", Coords[2], ")"
+ 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)"
+ print("But must be (5, 15, 35)")
####### GetNormal #######
face0 = faces[0]
vnorm = geompy.GetNormal(face0)
if vnorm is None:
- raise RuntimeError, "GetNormal(face0) failed"
+ raise RuntimeError("GetNormal(face0) failed")
else:
geompy.addToStudy(face0, "Face0")
geompy.addToStudy(vnorm, "Normale to Face0")
- print "\nNormale of face has been successfully obtained:"
+ 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 " 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
+ 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], ")"
+ 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)
# in one plane
Angle = geompy.GetAngle(OX, OXY)
- print "\nAngle between OX and OXY = ", Angle
+ 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"
+ 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
+ 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"
+ print(" Error: returned angle is", Angle, "while must be pi/4")
pass
# not in one plane
OXY_shift = geompy.MakeTranslation(OXY,10,-10,20)
Angle = geompy.GetAngle(OX, OXY_shift)
- print "Angle between OX and OXY_shift = ", Angle
+ print("Angle between OX and OXY_shift = ", Angle)
if math.fabs(Angle - 45.0) > 1e-05:
- print " Error: returned angle is", Angle, "while must be 45.0"
+ print(" Error: returned angle is", Angle, "while must be 45.0")
####### Position (LCS) #######
Pos = geompy.GetPosition(box)
- print "\nPosition(LCS) of box 10x30x70:"
- print "Origin: (", Pos[0], ", ", Pos[1], ", ", Pos[2], ")"
- print "Z axis: (", Pos[3], ", ", Pos[4], ", ", Pos[5], ")"
- print "X axis: (", Pos[6], ", ", Pos[7], ", ", Pos[8], ")"
+ print("\nPosition(LCS) of box 10x30x70:")
+ print("Origin: (", Pos[0], ", ", Pos[1], ", ", Pos[2], ")")
+ print("Z axis: (", Pos[3], ", ", Pos[4], ", ", Pos[5], ")")
+ print("X axis: (", Pos[6], ", ", Pos[7], ", ", Pos[8], ")")
####### KindOfShape #######
Kind = geompy.KindOfShape(box)
- print "\nKindOfShape(box 10x30x70):", Kind
+ print("\nKindOfShape(box 10x30x70):", Kind)
#if Kind[0] != geompy.kind.BOX:
# print "Error: returned type is", Kind[0], "while must be", geompy.kind.BOX
Kind = geompy.KindOfShape(p137)
- print "\nKindOfShape(p137):", Kind
+ print("\nKindOfShape(p137):", Kind)
if Kind[0] != geompy.kind.VERTEX:
- print " Error: returned type is", Kind[0], "while must be", geompy.kind.VERTEX
+ print(" Error: returned type is", Kind[0], "while must be", geompy.kind.VERTEX)
else:
dx = math.fabs(Kind[1] - 10)
dy = math.fabs(Kind[2] - 30)
dz = math.fabs(Kind[3] - 70)
if (dx + dy + dz) > 1e-5:
- print " Error: coordinates are (", Kind[1], ",", Kind[2], ",", Kind[3], ") while must be (10, 20, 30)"
+ print(" Error: coordinates are (", Kind[1], ",", Kind[2], ",", Kind[3], ") while must be (10, 20, 30)")
pass
