Salome HOME
New big test 'sprocket' came from issue #2380
[modules/shaper.git] / src / SketchPlugin / Test / TestMultiRotationWithParameter.py
1 ## Copyright (C) 2018-20xx  CEA/DEN, EDF R&D
2 ##
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.
7 ##
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.
12 ##
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
16 ##
17 ## See http:##www.salome-platform.org/ or
18 ## email : webmaster.salome@opencascade.com<mailto:webmaster.salome@opencascade.com>
19 ##
20
21 from GeomAPI import *
22 from GeomAlgoAPI import GeomAlgoAPI_ShapeTools
23 from SketchAPI import *
24 from math import *
25
26 from salome.shaper import model
27
28 model.begin()
29 partSet = model.moduleDocument()
30 Part_1 = model.addPart(partSet)
31 Part_1_doc = Part_1.document()
32 Param = model.addParameter(Part_1_doc, "N", "3")
33 Sketch_1 = model.addSketch(Part_1_doc, model.defaultPlane("XOY"))
34 SketchLine_1 = Sketch_1.addLine(30, 20, 10, 20)
35 SketchLine_2 = Sketch_1.addLine(10, 20, 10, 40)
36 SketchLine_3 = Sketch_1.addLine(10, 40, 30, 40)
37 SketchLine_4 = Sketch_1.addLine(30, 40, 30, 20)
38 SketchConstraintCoincidence_1 = Sketch_1.setCoincident(SketchLine_4.endPoint(), SketchLine_1.startPoint())
39 SketchConstraintCoincidence_2 = Sketch_1.setCoincident(SketchLine_1.endPoint(), SketchLine_2.startPoint())
40 SketchConstraintCoincidence_3 = Sketch_1.setCoincident(SketchLine_2.endPoint(), SketchLine_3.startPoint())
41 SketchConstraintCoincidence_4 = Sketch_1.setCoincident(SketchLine_3.endPoint(), SketchLine_4.startPoint())
42 SketchConstraintHorizontal_1 = Sketch_1.setHorizontal(SketchLine_1.result())
43 SketchConstraintVertical_1 = Sketch_1.setVertical(SketchLine_2.result())
44 SketchConstraintHorizontal_2 = Sketch_1.setHorizontal(SketchLine_3.result())
45 SketchConstraintVertical_2 = Sketch_1.setVertical(SketchLine_4.result())
46 SketchProjection_1 = Sketch_1.addProjection(model.selection("VERTEX", "PartSet/Origin"), False)
47 SketchPoint_1 = SketchProjection_1.createdFeature()
48 model.do()
49 SketchMultiRotation_1 = Sketch_1.addRotation([SketchLine_2.result()], SketchAPI_Point(SketchPoint_1).coordinates(), "360/N", "N")
50 Extrusion_1 = model.addExtrusion(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchLine_1r-SketchLine_2r-SketchLine_3r-SketchLine_4r")], model.selection(), 10, 0)
51 model.do()
52
53 # change parameters
54 N =[6, 5]
55 for n in N:
56     Param.setValue(n)
57     model.do()
58     Part_1_doc.setCurrentFeature(Extrusion_1.feature(), True)
59     model.do()
60
61 # check coordinates of center of mass of each line
62 REF_DATA = [[-25.4415255451, 18.7810749942],
63             [-25.7237275125, -18.3926573083],
64             [9.54338762502, -30.1483623542],
65             [31.6218654326, -0.240055331703]]
66 TOLERANCE = 1.e-7
67 aLines = SketchMultiRotation_1.rotated()
68 ind = 0
69 for lin in aLines:
70     curShape = lin.feature().lastResult().shape()
71     massCenter = GeomAlgoAPI_ShapeTools.centreOfMass(curShape)
72     assert(fabs(massCenter.x() - REF_DATA[ind][0]) < TOLERANCE and fabs(massCenter.y() - REF_DATA[ind][1]) < TOLERANCE and massCenter.z() == 0), "({}, {}, {}) != ({}, {}, {})".format(massCenter.x(), massCenter.y(), massCenter.z(), REF_DATA[ind][0], REF_DATA[ind][1], 0.0)
73     ind += 1
74
75 model.end()
76
77 assert(model.checkPythonDump())