1 # Copyright (C) 2018-2019 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 email : webmaster.salome@opencascade.com
21 from SketchAPI import *
23 from salome.shaper import model
29 def assertCircle(theEdge, theCenter, theRadius):
30 assert(theEdge.isCircle())
31 aCircle = theEdge.circle()
32 assert(aCircle.center().distance(theCenter) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(aCircle.center().x(), aCircle.center().y(), aCircle.center().z(), theCenter.x(), theCenter.y(), theCenter.z())
33 assert(math.fabs(aCircle.radius() - theRadius) < TOLERANCE), "Radius {} != {}".format(aCircle.radius(), theRadius)
35 def checkCircleEdge(theDocument, theEdgeName, theCenter, theRadius):
36 anEdge = model.addEdge(theDocument, [model.selection("EDGE", theEdgeName)])
37 aShape = anEdge.result().resultSubShapePair()[0].shape()
38 assert(aShape.isEdge())
39 assert(aShape.edge().isClosed())
40 assertCircle(aShape.edge(), theCenter, theRadius)
41 theDocument.removeFeature(anEdge.feature())
43 def checkCircleFace(theDocument, theFaceName, theCenter, theRadius):
44 aFaceFeature = model.addFace(theDocument, [model.selection("FACE", theFaceName)])
45 aShape = aFaceFeature.result().resultSubShapePair()[0].shape()
46 assert(aShape.isFace())
47 aFace = aShape.face();
48 aSubs = aFace.subShapes(GeomAPI.GeomAPI_Shape.EDGE)
49 assert(aSubs.size() == 1)
50 assertCircle(aSubs[0].edge(), theCenter, theRadius)
51 theDocument.removeFeature(aFaceFeature.feature())
53 def assertEllipse(theEdge, theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius, theNormal = None):
54 assert(theEdge.isEllipse())
55 anEllipse = theEdge.ellipse()
56 assert(anEllipse.firstFocus().distance(theFirstFocus) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(anEllipse.firstFocus().x(), anEllipse.firstFocus().y(), anEllipse.firstFocus().z(), theFirstFocus.x(), theFirstFocus.y(), theFirstFocus.z())
57 assert(anEllipse.secondFocus().distance(theSecondFocus) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(anEllipse.secondFocus().x(), anEllipse.secondFocus().y(), anEllipse.secondFocus().z(), theSecondFocus.x(), theSecondFocus.y(), theSecondFocus.z())
58 assert(math.fabs(anEllipse.majorRadius() - theMajorRadius) < TOLERANCE), "Major radius {} != {}".format(anEllipse.majorRadius(), theMajorRadius)
59 assert(math.fabs(anEllipse.minorRadius() - theMinorRadius) < TOLERANCE), "Minor radius {} != {}".format(anEllipse.minorRadius(), theMinorRadius)
61 center = GeomAPI.GeomAPI_Pnt((theFirstFocus.x() + theSecondFocus.x()) * 0.5, (theFirstFocus.y() + theSecondFocus.y()) * 0.5, (theFirstFocus.z() + theSecondFocus.z()) * 0.5)
62 assert(anEllipse.center().distance(center) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(anEllipse.center().x(), anEllipse.center().y(), anEllipse.center().z(), center.x(), center.y(), center.z())
64 if theNormal is not None:
65 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())
67 def checkEllipseEdge(theDocument, theEdgeName, theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius, theNormal):
68 anEdge = model.addEdge(theDocument, [model.selection("EDGE", theEdgeName)])
69 aShape = anEdge.result().resultSubShapePair()[0].shape()
70 assert(aShape.isEdge())
71 assertEllipse(aShape.edge(), theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius, theNormal)
72 theDocument.removeFeature(anEdge.feature())
74 def checkEllipseFace(theDocument, theFaceName, theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius):
75 aFaceFeature = model.addFace(theDocument, [model.selection("FACE", theFaceName)])
76 aShape = aFaceFeature.result().resultSubShapePair()[0].shape()
77 assert(aShape.isFace())
78 aFace = aShape.face();
79 aSubs = aFace.subShapes(GeomAPI.GeomAPI_Shape.EDGE)
80 assert(aSubs.size() == 1)
81 assertEllipse(aSubs[0].edge(), theFirstFocus, theSecondFocus, theMajorRadius, theMinorRadius)
82 theDocument.removeFeature(aFaceFeature.feature())
84 def assertCone(theCone, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
85 assert(theCone is not None)
86 assert(theCone.isSemiInfinite() == False)
87 assert(theCone.isInfinite() == False)
88 anApex = theCone.apex()
89 anAxis = theCone.axis()
90 assert(anApex.distance(theApex) < TOLERANCE), "({}, {}, {}) != expected ({}, {}, {})".format(anApex.x(), anApex.y(), anApex.z(), theApex.x(), theApex.y(), theApex.z())
91 assert(anAxis.isParallel(theAxis, TOLERANCE)), "dir({}, {}, {}) is not parallel to dir({}, {}, {})".format(anAxis.x(), anAxis.y(), anAxis.z(), theAxis.x(), theAxis.y(), theAxis.z())
92 assert(math.fabs(theCone.semiAngle() - theSemiAngle) < TOLERANCE), "SemiAngle {} != {}".format(theCone.semiAngle(), theSemiAngle)
93 assert(math.fabs(theCone.radius1() - theRadius1) < TOLERANCE), "Radius1 {} != {}".format(theCone.radius1(), theRadius1)
94 assert(math.fabs(theCone.radius2() - theRadius2) < TOLERANCE), "Radius2 {} != {}".format(theCone.radius2(), theRadius2)
95 assert(math.fabs(theCone.height() - theHeight) < TOLERANCE), "Height {} != {}".format(theCone.height(), theHeight)
97 def checkConeFace(theDocument, theFaceName, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
99 aFace = model.addFace(theDocument, [model.selection("FACE", theFaceName)])
100 aShape = aFace.result().resultSubShapePair()[0].shape()
101 assert(aShape.isFace())
102 assertCone(aShape.face().getCone(), theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
103 theDocument.removeFeature(aFace.feature())
105 def checkConeShell(theDocument, theFaceNames, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
106 # check conical shell
108 for name in theFaceNames:
109 aSelection.append(model.selection("FACE", name))
110 aShell = model.addShell(theDocument, aSelection)
111 aShape = aShell.result().resultSubShapePair()[0].shape()
112 assert(aShape.isShell())
113 assertCone(aShape.shell().getCone(), theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
114 theDocument.removeFeature(aShell.feature())
116 def checkConeSolid(theDocument, theFaceNames, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
119 for name in theFaceNames:
120 aSelection.append(model.selection("FACE", name))
121 aSolid = model.addSolid(theDocument, aSelection)
122 aShape = aSolid.result().resultSubShapePair()[0].shape()
123 assert(aShape.isSolid())
124 assertCone(aShape.solid().getCone(), theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
125 theDocument.removeFeature(aSolid.feature())
127 def checkConeAll(theDocument, theFeature, theFaceName, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight):
129 aShape = theFeature.result().resultSubShapePair()[0].shape()
130 assert(aShape.isSolid())
131 assertCone(aShape.solid().getCone(), theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
133 checkConeShell(theDocument, [theFaceName], theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
134 checkConeFace(theDocument, theFaceName, theApex, theAxis, theSemiAngle, theRadius1, theRadius2, theHeight)
136 def checkNonCone(theFeature):
137 aShape = theFeature.result().resultSubShapePair()[0].shape()
138 assert(aShape.isSolid())
139 assert(aShape.solid().getCone() is None)
141 def checkNonConeShell(theFeature):
142 aShape = theFeature.result().resultSubShapePair()[0].shape()
143 assert(aShape.isShell())
144 assert(aShape.shell().getCone() is None)
146 def checkSegment(theDocument, theEdgeName, theStartPoint, theEndPoint):
147 anEdgeFeature = model.addEdge(theDocument, [model.selection("EDGE", theEdgeName)])
148 aShape = anEdgeFeature.result().resultSubShapePair()[0].shape()
149 assert(aShape.isEdge())
150 anEdge = aShape.edge()
151 assert(anEdge.isLine())
152 assert(anEdge.firstPoint().distance(theStartPoint) < TOLERANCE)
153 assert(anEdge.lastPoint().distance(theEndPoint) < TOLERANCE)
154 theDocument.removeFeature(anEdgeFeature.feature())
156 def checkVertex(theDocument, theVertexName, theCoordinates):
157 aVertex = model.addVertex(theDocument, [model.selection("VERTEX", theVertexName)])
158 aShape = aVertex.result().resultSubShapePair()[0].shape()
159 assert(aShape.isVertex())
160 assert(aShape.vertex().point().distance(theCoordinates) < TOLERANCE)
161 theDocument.removeFeature(aVertex.feature())
163 def semiAngle(theRadius1, theRadius2, theHeight):
164 return math.atan(math.fabs(theRadius1 - theRadius2) / theHeight)
168 partSet = model.moduleDocument()
169 Part_1 = model.addPart(partSet)
170 Part_1_doc = Part_1.document()
171 ParamR1 = model.addParameter(Part_1_doc, "R1", "50")
172 ParamR2 = model.addParameter(Part_1_doc, "R2", "5")
173 ParamH = model.addParameter(Part_1_doc, "H", "70")
174 ParamShift = model.addParameter(Part_1_doc, "Shift", "5")
175 ParamAngle = model.addParameter(Part_1_doc, "Angle", "60")
176 Cone_1 = model.addCone(Part_1_doc, model.selection("VERTEX", "PartSet/Origin"), model.selection("EDGE", "PartSet/OZ"), "R1", "R2", "H")
180 aSemiAngle = semiAngle(ParamR1.value(), ParamR2.value(), ParamH.value())
181 anApex = GeomAPI.GeomAPI_Pnt(0, 0, ParamR1.value() / math.tan(aSemiAngle))
182 anAxis = GeomAPI.GeomAPI_Dir(0, 0, 1)
183 checkConeAll(Part_1_doc, Cone_1, "Cone_1_1/Face_1", anApex, anAxis, aSemiAngle, ParamR2.value(), ParamR1.value(), ParamH.value())
184 checkCircleFace(Part_1_doc, "Cone_1_1/Face_2", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
185 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_2]", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
186 checkCircleFace(Part_1_doc, "Cone_1_1/Face_3", GeomAPI.GeomAPI_Pnt(0, 0, 0), ParamR1.value())
187 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_3]", GeomAPI.GeomAPI_Pnt(0, 0, 0), ParamR1.value())
188 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()))
190 # Test 2. Update cone radii
193 aSemiAngle = semiAngle(ParamR1.value(), ParamR2.value(), ParamH.value())
195 checkConeAll(Part_1_doc, Cone_1, "Cone_1_1/Face_1", anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
196 checkCircleFace(Part_1_doc, "Cone_1_1/Face_2", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
197 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_2]", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
198 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()))
199 checkVertex(Part_1_doc, "Cone_1_1/Vertex_2", GeomAPI.GeomAPI_Pnt(0, 0, 0))
204 aSemiAngle = semiAngle(ParamR1.value(), ParamR2.value(), ParamH.value())
205 anApex.setZ(-ParamR1.value() / math.tan(aSemiAngle))
206 checkConeAll(Part_1_doc, Cone_1, "Cone_1_1/Face_1", anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
207 checkCircleFace(Part_1_doc, "Cone_1_1/Face_2", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
208 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_2]", GeomAPI.GeomAPI_Pnt(0, 0, ParamH.value()), ParamR2.value())
209 checkCircleFace(Part_1_doc, "Cone_1_1/Face_3", GeomAPI.GeomAPI_Pnt(0, 0, 0), ParamR1.value())
210 checkCircleEdge(Part_1_doc, "[Cone_1_1/Face_1][Cone_1_1/Face_3]", GeomAPI.GeomAPI_Pnt(0, 0, 0), ParamR1.value())
211 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()))
213 # Test 3. Translate cone
214 Translation_1 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Cone_1_1")], model.selection("EDGE", "PartSet/OX"), "Shift")
215 anApex.setX(anApex.x() + ParamShift.value())
216 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())
217 checkCircleFace(Part_1_doc, "Translation_1_1/MF:Translated&Cone_1_1/Face_2", GeomAPI.GeomAPI_Pnt(ParamShift.value(), 0, ParamH.value()), ParamR2.value())
218 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())
219 checkCircleFace(Part_1_doc, "Translation_1_1/MF:Translated&Cone_1_1/Face_3", GeomAPI.GeomAPI_Pnt(ParamShift.value(), 0, 0), ParamR1.value())
220 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())
221 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()))
223 # Test 4. Rotate cone
224 Rotation_1 = model.addRotation(Part_1_doc, [model.selection("SOLID", "Translation_1_1")], model.selection("EDGE", "PartSet/OY"), "Angle")
225 anAngle = ParamAngle.value() * math.pi / 180.0
226 anAxis = GeomAPI.GeomAPI_Dir(math.sin(anAngle), 0, math.cos(anAngle))
227 x, z = anApex.x(), anApex.z()
228 anApex.setX(x * math.cos(anAngle) + z * math.sin(anAngle))
229 anApex.setZ(-x * math.sin(anAngle) + z * math.cos(anAngle))
230 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())
231 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))
232 checkCircleFace(Part_1_doc, "Rotation_1_1/MF:Rotated&Cone_1_1/Face_2", aCenter, ParamR2.value())
233 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())
234 aCenter = GeomAPI.GeomAPI_Pnt(ParamShift.value() * math.cos(anAngle), 0, -ParamShift.value() * math.sin(anAngle))
235 checkCircleFace(Part_1_doc, "Rotation_1_1/MF:Rotated&Cone_1_1/Face_3", aCenter, ParamR1.value())
236 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())
238 # Test 5. Split cone by plane and check conical shell and elliptic face
239 Plane_4 = model.addPlane(Part_1_doc, model.selection("FACE", "PartSet/YOZ"), 20, False)
240 Partition_1 = model.addPartition(Part_1_doc, [model.selection("SOLID", "Rotation_1_1"), model.selection("FACE", "Plane_1")])
241 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())
243 aFirstFocus = GeomAPI.GeomAPI_Pnt(20, 0, 31.062397266842858)
244 aSecondFocus = GeomAPI.GeomAPI_Pnt(20, 0, -1.0935246846933797)
245 aMajorRadius = 27.91915871311068
246 aMinorRadius = 22.824955511666207
247 aNormal = GeomAPI.GeomAPI_Dir(1, 0, 0)
248 checkEllipseFace(Part_1_doc, "_weak_name_1_Partition_1_1_2", aFirstFocus, aSecondFocus, aMajorRadius, aMinorRadius)
249 checkEllipseEdge(Part_1_doc, "Partition_1_1_1/Generated_Edge&Cone_1_1/Face_1", aFirstFocus, aSecondFocus, aMajorRadius, aMinorRadius, aNormal)
251 # Test 6. Compose a conical solid
252 Solid_1_objects = ["Rotation_1_1/MF:Rotated&Cone_1_1/Face_3",
253 "Partition_1_1_1/Modified_Face&Cone_1_1/Face_1",
254 "Partition_1_1_2/Modified_Face&Cone_1_1/Face_1",
255 "Rotation_1_1/MF:Rotated&Cone_1_1/Face_2"]
256 checkConeSolid(Part_1_doc, Solid_1_objects, anApex, anAxis, aSemiAngle, ParamR1.value(), ParamR2.value(), ParamH.value())
258 # Test 7. Check non-cone
259 Cone_2 = model.addCone(Part_1_doc, model.selection("VERTEX", "PartSet/Origin"), model.selection("EDGE", "PartSet/OZ"), 10, 5, 10)
260 Cone_3 = model.addCone(Part_1_doc, model.selection("VERTEX", "PartSet/Origin"), model.selection("EDGE", "PartSet/OZ"), 5, 10, 20)
261 Fuse_1 = model.addFuse(Part_1_doc, [model.selection("SOLID", "Cone_2_1"), model.selection("SOLID", "Cone_3_1")], True)
264 Solid_2_objects = [model.selection("FACE", "Fuse_1_1/Modified_Face&Cone_2_1/Face_3&Cone_3_1/Face_3"),
265 model.selection("FACE", "Fuse_1_1/Modified_Face&Cone_2_1/Face_1"),
266 model.selection("FACE", "Fuse_1_1/Modified_Face&Cone_3_1/Face_1"),
267 model.selection("FACE", "Cone_3_1/Face_2")]
268 Solid_2 = model.addSolid(Part_1_doc, Solid_2_objects)
269 checkNonCone(Solid_2)
272 # in order to use study objects once again, undo Test7 actions
275 # Test 8. Check non-conical shell
277 Shell_1_objects = [model.selection("FACE", "Rotation_1_1/MF:Rotated&Cone_1_1/Face_3"),
278 model.selection("FACE", "Partition_1_1_1/Modified_Face&Cone_1_1/Face_1"),
279 model.selection("FACE", "Partition_1_1_2/Modified_Face&Cone_1_1/Face_1"),
280 model.selection("FACE", "Rotation_1_1/MF:Rotated&Cone_1_1/Face_2")]
281 Shell_1 = model.addShell(Part_1_doc, Shell_1_objects)
282 checkNonConeShell(Shell_1)
284 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")])
285 checkNonConeShell(Shell_2)
288 # in order to use study objects once again, undo Test8 actions
291 # Test 9. Check error on conversion to wrong type of curve
293 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]")])
294 aShape = anEdge.result().resultSubShapePair()[0].shape()
295 assert(aShape.isEdge())
296 assert(aShape.edge().ellipse() is None)
297 assert(aShape.edge().line() is None)
304 anEdge = model.addEdge(Part_1_doc, [model.selection("EDGE", "[Partition_1_1_2/Modified_Face&Cone_1_1/Face_1][weak_name_2]")])
305 aShape = anEdge.result().resultSubShapePair()[0].shape()
306 assert(aShape.isEdge())
307 assert(aShape.edge().circle() is None)
308 assert(aShape.isSelfIntersected() == False)
310 # Test 10. Test cone constructors
311 apex = GeomAPI_Pnt(10, 0, 0)
312 dir = GeomAPI_Dir(1, 0, 0)
313 semiAngle = math.pi / 4
314 cone = GeomAPI_Cone(apex, dir, semiAngle)
315 assert(cone.location().distance(apex) < TOLERANCE)
318 cone = GeomAPI_Cone(apex, dir, semiAngle, radius)
319 assert(cone.location().distance(apex) < TOLERANCE)