1 ## Copyright (C) 2018-20xx CEA/DEN, EDF R&D
3 ## This library is free software; you can redistribute it and/or
4 ## modify it under the terms of the GNU Lesser General Public
5 ## License as published by the Free Software Foundation; either
6 ## version 2.1 of the License, or (at your option) any later version.
8 ## This library is distributed in the hope that it will be useful,
9 ## but WITHOUT ANY WARRANTY; without even the implied warranty of
10 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 ## Lesser General Public License for more details.
13 ## You should have received a copy of the GNU Lesser General Public
14 ## License along with this library; if not, write to the Free Software
15 ## Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 ## See http:##www.salome-platform.org/ or
18 ## email : webmaster.salome@opencascade.com<mailto:webmaster.salome@opencascade.com>
22 from SketchAPI import *
24 from salome.shaper import model
30 def assertCircle(theEdge, theCenter, theRadius):
31 assert(theEdge.isCircle())
32 aCircle = theEdge.circle()
33 assert(aCircle.center().distance(theCenter) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(aCircle.center().x(), aCircle.center().y(), aCircle.center().z(), theCenter.x(), theCenter.y(), theCenter.z())
34 assert(math.fabs(aCircle.radius() - theRadius) < TOLERANCE), "Radius {} != {}".format(aCircle.radius(), theRadius)
36 def checkCircleEdge(theDocument, theEdgeName, theCenter, theRadius):
37 anEdge = model.addEdge(theDocument, [model.selection("EDGE", theEdgeName)])
38 aShape = anEdge.result().resultSubShapePair()[0].shape()
39 assert(aShape.isEdge())
40 assert(aShape.edge().isClosed())
41 assertCircle(aShape.edge(), theCenter, theRadius)
42 theDocument.removeFeature(anEdge.feature())
44 def checkCircleFace(theDocument, theFaceName, theCenter, theRadius):
45 aFaceFeature = model.addFace(theDocument, [model.selection("FACE", theFaceName)])
46 aShape = aFaceFeature.result().resultSubShapePair()[0].shape()
47 assert(aShape.isFace())
48 aFace = aShape.face();
49 aSubs = aFace.subShapes(GeomAPI.GeomAPI_Shape.EDGE)
50 assert(aSubs.size() == 1)
51 assertCircle(aSubs[0].edge(), theCenter, theRadius)
52 theDocument.removeFeature(aFaceFeature.feature())
54 def assertEllipse(theEdge, theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius, theNormal = None):
55 assert(theEdge.isEllipse())
56 anEllipse = theEdge.ellipse()
57 assert(anEllipse.firstFocus().distance(theFirstFocus) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(anEllipse.firstFocus().x(), anEllipse.firstFocus().y(), anEllipse.firstFocus().z(), theFirstFocus.x(), theFirstFocus.y(), theFirstFocus.z())
58 assert(anEllipse.secondFocus().distance(theSecondFocus) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(anEllipse.secondFocus().x(), anEllipse.secondFocus().y(), anEllipse.secondFocus().z(), theSecondFocus.x(), theSecondFocus.y(), theSecondFocus.z())
59 assert(math.fabs(anEllipse.majorRadius() - theMajorRadius) < TOLERANCE), "Major radius {} != {}".format(anEllipse.majorRadius(), theMajorRadius)
60 assert(math.fabs(anEllipse.minorRadius() - theMinorRadius) < TOLERANCE), "Minor radius {} != {}".format(anEllipse.minorRadius(), theMinorRadius)
62 center = GeomAPI.GeomAPI_Pnt((theFirstFocus.x() + theSecondFocus.x()) * 0.5, (theFirstFocus.y() + theSecondFocus.y()) * 0.5, (theFirstFocus.z() + theSecondFocus.z()) * 0.5)
63 assert(anEllipse.center().distance(center) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(anEllipse.center().x(), anEllipse.center().y(), anEllipse.center().z(), center.x(), center.y(), center.z())
65 if theNormal is not None:
66 assert(math.fabs(anEllipse.normal().dot(theNormal) - 1.) < TOLERANCE), "Normal ({}, {}, {}) != ({}, {}, {})".format(anEllipse.normal().x(), anEllipse.normal().y(), anEllipse.normal().z(), theNormal.x(), theNormal.y(), theNormal.z())
68 def checkEllipseEdge(theDocument, theEdgeName, theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius, theNormal):
69 anEdge = model.addEdge(theDocument, [model.selection("EDGE", theEdgeName)])
70 aShape = anEdge.result().resultSubShapePair()[0].shape()
71 assert(aShape.isEdge())
72 assertEllipse(aShape.edge(), theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius, theNormal)
73 theDocument.removeFeature(anEdge.feature())
75 def checkEllipseFace(theDocument, theFaceName, theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius):
76 aFaceFeature = model.addFace(theDocument, [model.selection("FACE", theFaceName)])
77 aShape = aFaceFeature.result().resultSubShapePair()[0].shape()
78 assert(aShape.isFace())
79 aFace = aShape.face();
80 aSubs = aFace.subShapes(GeomAPI.GeomAPI_Shape.EDGE)
81 assert(aSubs.size() == 1)
82 assertEllipse(aSubs[0].edge(), theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius)
83 theDocument.removeFeature(aFaceFeature.feature())
85 def assertCone(theCone, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
86 assert(theCone is not None)
87 assert(theCone.isSemiInfinite() == False)
88 assert(theCone.isInfinite() == False)
89 anApex = theCone.apex()
90 anAxis = theCone.axis()
91 assert(anApex.distance(theApex) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(anApex.x(), anApex.y(), anApex.z(), theApex.x(), theApex.y(), theApex.z())
92 assert(anAxis.isParallel(theAxis, TOLERANCE)), "dir({}, {}, {}) is not parallel to dir({}, {}, {})".format(anAxis.x(), anAxis.y(), anAxis.z(), theAxis.x(), theAxis.y(), theAxis.z())
93 assert(math.fabs(theCone.semiAngle() - theSemiAngle) < TOLERANCE), "SemiAngle {} != {}".format(theCone.semiAngle(), theSemiAngle)
94 assert(math.fabs(theCone.radius1() - theRadius1) < TOLERANCE), "Radius1 {} != {}".format(theCone.radius1(), theRadius1)
95 assert(math.fabs(theCone.radius2() - theRadius2) < TOLERANCE), "Radius2 {} != {}".format(theCone.radius2(), theRadius2)
96 assert(math.fabs(theCone.height() - theHeight) < TOLERANCE), "Height {} != {}".format(theCone.height(), theHeight)
98 def checkConeFace(theDocument, theFaceName, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
100 aFace = model.addFace(theDocument, [model.selection("FACE", theFaceName)])
101 aShape = aFace.result().resultSubShapePair()[0].shape()
102 assert(aShape.isFace())
103 assertCone(aShape.face().getCone(), theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
104 theDocument.removeFeature(aFace.feature())
106 def checkConeShell(theDocument, theFaceNames, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
107 # check conical shell
109 for name in theFaceNames:
110 aSelection.append(model.selection("FACE", name))
111 aShell = model.addShell(theDocument, aSelection)
112 aShape = aShell.result().resultSubShapePair()[0].shape()
113 assert(aShape.isShell())
114 assertCone(aShape.shell().getCone(), theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
115 theDocument.removeFeature(aShell.feature())
117 def checkConeSolid(theDocument, theFaceNames, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
120 for name in theFaceNames:
121 aSelection.append(model.selection("FACE", name))
122 aSolid = model.addSolid(theDocument, aSelection)
123 aShape = aSolid.result().resultSubShapePair()[0].shape()
124 assert(aShape.isSolid())
125 assertCone(aShape.solid().getCone(), theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
126 theDocument.removeFeature(aSolid.feature())
128 def checkConeAll(theDocument, theFeature, theFaceName, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
130 aShape = theFeature.result().resultSubShapePair()[0].shape()
131 assert(aShape.isSolid())
132 assertCone(aShape.solid().getCone(), theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
134 checkConeShell(theDocument, [theFaceName], theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
135 checkConeFace(theDocument, theFaceName, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
137 def checkNonCone(theFeature):
138 aShape = theFeature.result().resultSubShapePair()[0].shape()
139 assert(aShape.isSolid())
140 assert(aShape.solid().getCone() is None)
142 def checkNonConeShell(theFeature):
143 aShape = theFeature.result().resultSubShapePair()[0].shape()
144 assert(aShape.isShell())
145 assert(aShape.shell().getCone() is None)
147 def checkSegment(theDocument, theEdgeName, theStartPoint, theEndPoint):
148 anEdgeFeature = model.addEdge(theDocument, [model.selection("EDGE", theEdgeName)])
149 aShape = anEdgeFeature.result().resultSubShapePair()[0].shape()
150 assert(aShape.isEdge())
151 anEdge = aShape.edge()
152 assert(anEdge.isLine())
153 assert(anEdge.firstPoint().distance(theStartPoint) < TOLERANCE)
154 assert(anEdge.lastPoint().distance(theEndPoint) < TOLERANCE)
155 theDocument.removeFeature(anEdgeFeature.feature())
157 def checkVertex(theDocument, theVertexName, theCoordinates):
158 aVertex = model.addVertex(theDocument, [model.selection("VERTEX", theVertexName)])
159 aShape = aVertex.result().resultSubShapePair()[0].shape()
160 assert(aShape.isVertex())
161 assert(aShape.vertex().point().distance(theCoordinates) < TOLERANCE)
162 theDocument.removeFeature(aVertex.feature())
164 def semiAngle(theRadius1, theRadius2, theHeight):
165 return math.atan(math.fabs(theRadius1 - theRadius2) / theHeight)
169 partSet = model.moduleDocument()
170 Part_1 = model.addPart(partSet)
171 Part_1_doc = Part_1.document()
172 ParamR1 = model.addParameter(Part_1_doc, "R1", "50")
173 ParamR2 = model.addParameter(Part_1_doc, "R2", "5")
174 ParamH = model.addParameter(Part_1_doc, "H", "70")
175 ParamShift = model.addParameter(Part_1_doc, "Shift", "5")
176 ParamAngle = model.addParameter(Part_1_doc, "Angle", "60")
177 Cone_1 = model.addCone(Part_1_doc, model.selection("VERTEX", "PartSet/Origin"), model.selection("EDGE", "PartSet/OZ"), "R1", "R2", "H")
181 aSemiAngle = semiAngle(ParamR1.value(), ParamR2.value(), ParamH.value())
182 anApex = GeomAPI.GeomAPI_Pnt(0, 0, ParamR1.value() / math.tan(aSemiAngle))
183 anAxis = GeomAPI.GeomAPI_Dir(0, 0, 1)
184 checkConeAll(Part_1_doc, Cone_1, "Cone_1_1/Face_1", anApex, anAxis, aSemiAngle, ParamR2.value(), ParamR1.value(), ParamH.value())
185 checkCircleFace(Part_1_doc, "Cone_1_1/Face_2", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
186 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_2]", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
187 checkCircleFace(Part_1_doc, "Cone_1_1/Face_3", GeomAPI.GeomAPI_Pnt(0, 0, 0), ParamR1.value())
188 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_3]", GeomAPI.GeomAPI_Pnt(0, 0, 0), ParamR1.value())
189 checkSegment(Part_1_doc, "[Cone_1_1/Face_1][weak_name_3]", GeomAPI.GeomAPI_Pnt(ParamR1.value(), 0, 0), GeomAPI.GeomAPI_Pnt(ParamR2.value(), 0, ParamH.value()))
191 # Test 2. Update cone radii
194 aSemiAngle = semiAngle(ParamR1.value(), ParamR2.value(), ParamH.value())
196 checkConeAll(Part_1_doc, Cone_1, "Cone_1_1/Face_1", anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
197 checkCircleFace(Part_1_doc, "Cone_1_1/Face_2", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
198 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_2]", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
199 checkSegment(Part_1_doc, "[Cone_1_1/Face_1][weak_name_3]", GeomAPI.GeomAPI_Pnt(0, 0, 0), GeomAPI.GeomAPI_Pnt(ParamR2.value(), 0, ParamH.value()))
200 checkVertex(Part_1_doc, "Cone_1_1/Vertex_2", GeomAPI.GeomAPI_Pnt(0, 0, 0))
205 aSemiAngle = semiAngle(ParamR1.value(), ParamR2.value(), ParamH.value())
206 anApex.setZ(-ParamR1.value() / math.tan(aSemiAngle))
207 checkConeAll(Part_1_doc, Cone_1, "Cone_1_1/Face_1", anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
208 checkCircleFace(Part_1_doc, "Cone_1_1/Face_2", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
209 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_2]", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
210 checkCircleFace(Part_1_doc, "Cone_1_1/Face_3", GeomAPI.GeomAPI_Pnt(0, 0, 0), ParamR1.value())
211 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_3]", GeomAPI.GeomAPI_Pnt(0, 0, 0), ParamR1.value())
212 checkSegment(Part_1_doc, "[Cone_1_1/Face_1][weak_name_3]", GeomAPI.GeomAPI_Pnt(ParamR1.value(), 0, 0), GeomAPI.GeomAPI_Pnt(ParamR2.value(), 0, ParamH.value()))
214 # Test 3. Translate cone
215 Translation_1 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Cone_1_1")], model.selection("EDGE", "PartSet/OX"), "Shift")
216 anApex.setX(anApex.x() + ParamShift.value())
217 checkConeAll(Part_1_doc, Translation_1, "Translation_1_1/MF:Translated&Cone_1_1/Face_1", anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
218 checkCircleFace(Part_1_doc, "Translation_1_1/MF:Translated&Cone_1_1/Face_2", GeomAPI.GeomAPI_Pnt(ParamShift.value(), 0, ParamH.value()), ParamR2.value())
219 checkCircleEdge(Part_1_doc, "[Translation_1_1/MF:Translated&Cone_1_1/Face_1][Translation_1_1/MF:Translated&Cone_1_1/Face_2]", GeomAPI.GeomAPI_Pnt(ParamShift.value(), 0, ParamH.value()), ParamR2.value())
220 checkCircleFace(Part_1_doc, "Translation_1_1/MF:Translated&Cone_1_1/Face_3", GeomAPI.GeomAPI_Pnt(ParamShift.value(), 0, 0), ParamR1.value())
221 checkCircleEdge(Part_1_doc, "[Translation_1_1/MF:Translated&Cone_1_1/Face_1][Translation_1_1/MF:Translated&Cone_1_1/Face_3]", GeomAPI.GeomAPI_Pnt(ParamShift.value(), 0, 0), ParamR1.value())
222 checkSegment(Part_1_doc, "[Translation_1_1/MF:Translated&Cone_1_1/Face_1][weak_name_3]", GeomAPI.GeomAPI_Pnt(ParamR1.value() + ParamShift.value(), 0, 0), GeomAPI.GeomAPI_Pnt(ParamR2.value() + ParamShift.value(), 0, ParamH.value()))
224 # Test 4. Rotate cone
225 Rotation_1 = model.addRotation(Part_1_doc, [model.selection("SOLID", "Translation_1_1")], model.selection("EDGE", "PartSet/OY"), "Angle")
226 anAngle = ParamAngle.value() * math.pi / 180.0
227 anAxis = GeomAPI.GeomAPI_Dir(math.sin(anAngle), 0, math.cos(anAngle))
228 x, z = anApex.x(), anApex.z()
229 anApex.setX(x * math.cos(anAngle) + z * math.sin(anAngle))
230 anApex.setZ(-x * math.sin(anAngle) + z * math.cos(anAngle))
231 checkConeAll(Part_1_doc, Rotation_1, "Rotation_1_1/MF:Rotated&Cone_1_1/Face_1", anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
232 aCenter = GeomAPI.GeomAPI_Pnt(ParamShift.value() * math.cos(anAngle) + ParamH.value() * math.sin(anAngle), 0, -ParamShift.value() * math.sin(anAngle) + ParamH.value() * math.cos(anAngle))
233 checkCircleFace(Part_1_doc, "Rotation_1_1/MF:Rotated&Cone_1_1/Face_2", aCenter, ParamR2.value())
234 checkCircleEdge(Part_1_doc, "[Rotation_1_1/MF:Rotated&Cone_1_1/Face_1][Rotation_1_1/MF:Rotated&Cone_1_1/Face_2]", aCenter, ParamR2.value())
235 aCenter = GeomAPI.GeomAPI_Pnt(ParamShift.value() * math.cos(anAngle), 0, -ParamShift.value() * math.sin(anAngle))
236 checkCircleFace(Part_1_doc, "Rotation_1_1/MF:Rotated&Cone_1_1/Face_3", aCenter, ParamR1.value())
237 checkCircleEdge(Part_1_doc, "[Rotation_1_1/MF:Rotated&Cone_1_1/Face_1][Rotation_1_1/MF:Rotated&Cone_1_1/Face_3]", aCenter, ParamR1.value())
239 # Test 5. Split cone by plane and check conical shell and elliptic face
240 Plane_4 = model.addPlane(Part_1_doc, model.selection("FACE", "PartSet/YOZ"), 20, False)
241 Partition_1 = model.addPartition(Part_1_doc, [model.selection("SOLID", "Rotation_1_1"), model.selection("FACE", "Plane_1")])
242 checkConeShell(Part_1_doc, ["Partition_1_1_1/Modified_Face&Cone_1_1/Face_1", "Partition_1_1_2/Modified_Face&Cone_1_1/Face_1"], anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
244 aFirstFocus = GeomAPI.GeomAPI_Pnt(20, 0, 31.062397266842858)
245 aSecondFocus = GeomAPI.GeomAPI_Pnt(20, 0, -1.0935246846933797)
246 aMajorRadius = 27.91915871311068
247 aMinorRadius = 22.824955511666207
248 aNormal = GeomAPI.GeomAPI_Dir(1, 0, 0)
249 checkEllipseFace(Part_1_doc, "_weak_name_1_Partition_1_1_2", aFirstFocus, aSecondFocus, aMajorRadius, aMinorRadius)
250 checkEllipseEdge(Part_1_doc, "Partition_1_1_1/Generated_Edge&Cone_1_1/Face_1", aFirstFocus, aSecondFocus, aMajorRadius, aMinorRadius, aNormal)
252 # Test 6. Compose a conical solid
253 Solid_1_objects = ["Rotation_1_1/MF:Rotated&Cone_1_1/Face_3",
254 "Partition_1_1_1/Modified_Face&Cone_1_1/Face_1",
255 "Partition_1_1_2/Modified_Face&Cone_1_1/Face_1",
256 "Rotation_1_1/MF:Rotated&Cone_1_1/Face_2"]
257 checkConeSolid(Part_1_doc, Solid_1_objects, anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
259 # Test 7. Check non-cone
260 Cone_2 = model.addCone(Part_1_doc, model.selection("VERTEX", "PartSet/Origin"), model.selection("EDGE", "PartSet/OZ"), 10, 5, 10)
261 Cone_3 = model.addCone(Part_1_doc, model.selection("VERTEX", "PartSet/Origin"), model.selection("EDGE", "PartSet/OZ"), 5, 10, 20)
262 Fuse_1 = model.addFuse(Part_1_doc, [model.selection("SOLID", "Cone_2_1"), model.selection("SOLID", "Cone_3_1")], True)
263 Solid_2_objects = [model.selection("FACE", "Fuse_1_1/Modified_Face&Cone_2_1/Face_3&Cone_3_1/Face_3"),
264 model.selection("FACE", "Fuse_1_1/Modified_Face&Cone_2_1/Face_1"),
265 model.selection("FACE", "Fuse_1_1/Modified_Face&Cone_3_1/Face_1"),
266 model.selection("FACE", "Cone_3_1/Face_2")]
267 Solid_2 = model.addSolid(Part_1_doc, Solid_2_objects)
268 checkNonCone(Solid_2)
270 # Test 8. Check non-conical shell
271 Shell_1_objects = [model.selection("FACE", "Rotation_1_1/MF:Rotated&Cone_1_1/Face_3"),
272 model.selection("FACE", "Partition_1_1_1/Modified_Face&Cone_1_1/Face_1"),
273 model.selection("FACE", "Partition_1_1_2/Modified_Face&Cone_1_1/Face_1"),
274 model.selection("FACE", "Rotation_1_1/MF:Rotated&Cone_1_1/Face_2")]
275 Shell_1 = model.addShell(Part_1_doc, Shell_1_objects)
276 checkNonConeShell(Shell_1)
278 Shell_2 = model.addShell(Part_1_doc, [model.selection("FACE", "Fuse_1_1/Modified_Face&Cone_2_1/Face_1"), model.selection("FACE", "Fuse_1_1/Modified_Face&Cone_3_1/Face_1")])
279 checkNonConeShell(Shell_2)
281 # Test 9. Check error on conversion to wrong type of curve
282 anEdge = model.addEdge(Part_1_doc, [model.selection("EDGE", "[Partition_1_1_2/Modified_Face&Cone_1_1/Face_1][Rotation_1_1/MF:Rotated&Cone_1_1/Face_2]")])
283 aShape = anEdge.result().resultSubShapePair()[0].shape()
284 assert(aShape.isEdge())
285 assert(aShape.edge().ellipse() is None)
286 assert(aShape.edge().line() is None)
288 anEdge = model.addEdge(Part_1_doc, [model.selection("EDGE", "[Partition_1_1_2/Modified_Face&Cone_1_1/Face_1][weak_name_2]")])
289 aShape = anEdge.result().resultSubShapePair()[0].shape()
290 assert(aShape.isEdge())
291 assert(aShape.edge().circle() is None)
292 assert(aShape.isSelfIntersected() == False)
294 # Test 10. Test cone constructors
295 apex = GeomAPI_Pnt(10, 0, 0)
296 dir = GeomAPI_Dir(1, 0, 0)
297 semiAngle = math.pi / 4
298 cone = GeomAPI_Cone(apex, dir, semiAngle)
299 assert(cone.location().distance(apex) < TOLERANCE)
302 cone = GeomAPI_Cone(apex, dir, semiAngle, radius)
303 assert(cone.location().distance(apex) < TOLERANCE)