Salome HOME
Copyright update 2021
[modules/shaper.git] / src / FeaturesPlugin / Test / TestMeasurementLength.py
1 # Copyright (C) 2014-2021  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 email : webmaster.salome@opencascade.com
18 #
19
20 from salome.shaper import model
21
22 import math
23
24 model.begin()
25 partSet = model.moduleDocument()
26 Part_1 = model.addPart(partSet)
27 Part_1_doc = Part_1.document()
28 ParamR = model.addParameter(Part_1_doc, "R", "50")
29 ParamL = model.addParameter(Part_1_doc, "L", "40")
30 ParamH = model.addParameter(Part_1_doc, "H", "20")
31 Point_2 = model.addPoint(Part_1_doc, 0, 100, 100)
32 Plane_4 = model.addPlane(Part_1_doc, model.selection("EDGE", "PartSet/OX"), model.selection("VERTEX", "Point_1"), False)
33 Sketch_1 = model.addSketch(Part_1_doc, model.defaultPlane("XOY"))
34 SketchCircle_1 = Sketch_1.addCircle(50, 50, 50)
35 SketchConstraintRadius_1 = Sketch_1.setRadius(SketchCircle_1.results()[1], "R")
36 SketchProjection_1 = Sketch_1.addProjection(model.selection("EDGE", "PartSet/OX"), False)
37 SketchLine_1 = SketchProjection_1.createdFeature()
38 SketchConstraintTangent_1 = Sketch_1.setTangent(SketchLine_1.result(), SketchCircle_1.results()[1])
39 SketchProjection_2 = Sketch_1.addProjection(model.selection("EDGE", "PartSet/OY"), False)
40 SketchLine_2 = SketchProjection_2.createdFeature()
41 SketchConstraintTangent_2 = Sketch_1.setTangent(SketchLine_2.result(), SketchCircle_1.results()[1])
42 SketchLine_3 = Sketch_1.addLine(70, 40, 30, 40)
43 SketchLine_4 = Sketch_1.addLine(30, 40, 30, 60)
44 SketchLine_5 = Sketch_1.addLine(30, 60, 70, 60)
45 SketchLine_6 = Sketch_1.addLine(70, 60, 70, 40)
46 SketchConstraintCoincidence_1 = Sketch_1.setCoincident(SketchLine_6.endPoint(), SketchLine_3.startPoint())
47 SketchConstraintCoincidence_2 = Sketch_1.setCoincident(SketchLine_3.endPoint(), SketchLine_4.startPoint())
48 SketchConstraintCoincidence_3 = Sketch_1.setCoincident(SketchLine_4.endPoint(), SketchLine_5.startPoint())
49 SketchConstraintCoincidence_4 = Sketch_1.setCoincident(SketchLine_5.endPoint(), SketchLine_6.startPoint())
50 SketchConstraintHorizontal_1 = Sketch_1.setHorizontal(SketchLine_3.result())
51 SketchConstraintVertical_1 = Sketch_1.setVertical(SketchLine_4.result())
52 SketchConstraintHorizontal_2 = Sketch_1.setHorizontal(SketchLine_5.result())
53 SketchConstraintVertical_2 = Sketch_1.setVertical(SketchLine_6.result())
54 SketchConstraintLength_1 = Sketch_1.setLength(SketchLine_3.result(), "L")
55 SketchConstraintLength_2 = Sketch_1.setLength(SketchLine_4.result(), "H")
56 SketchConstraintDistanceHorizontal_1 = Sketch_1.setHorizontalDistance(SketchLine_4.startPoint(), SketchCircle_1.center(), "L/2")
57 SketchConstraintDistanceVertical_1 = Sketch_1.setVerticalDistance(SketchLine_4.startPoint(), SketchCircle_1.center(), "H/2")
58 model.do()
59 Extrusion_1 = model.addExtrusion(Part_1_doc, [model.selection("FACE", "Sketch_1/Face-SketchCircle_1_2r-SketchLine_3f-SketchLine_4f-SketchLine_5f-SketchLine_6f")], model.selection(), model.selection("FACE", "Plane_1"), 0, model.selection(), 10)
60 model.do()
61
62 TOLERANCE = 1.e-6
63
64 # reference data
65 REF_DATA = [("[Extrusion_1_1/From_Face][Extrusion_1_1/Generated_Face&Sketch_1/SketchLine_6]", ParamH.value()),
66             ("[Extrusion_1_1/From_Face][Extrusion_1_1/Generated_Face&Sketch_1/SketchLine_3]", ParamL.value()),
67             ("[Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/From_Face]", ParamR.value() * 2.0 * math.pi),
68             ("[Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face]", 381.981436664),
69             ("([Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/From_Face])([Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face])([Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/From_Face])([Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2][Extrusion_1_1/To_Face])", 60),
70             ("Sketch_1/SketchLine_3", ParamL.value()),
71             ("Sketch_1/SketchLine_4", ParamH.value()),
72             ("Sketch_1/SketchCircle_1_2", ParamR.value() * 2.0 * math.pi)
73            ]
74
75 for ref in REF_DATA:
76     length = model.measureLength(Part_1_doc, model.selection("EDGE", ref[0]))
77     assert(math.fabs(length - ref[1]) < TOLERANCE), "Length {} differs from expected value {}".format(length, ref[1])
78
79 # select incorrect data
80 length = model.measureLength(Part_1_doc, model.selection("FACE", "Extrusion_1_1/Generated_Face&Sketch_1/SketchCircle_1_2"))
81 assert(length == -1)
82
83 model.end()