]> SALOME platform Git repositories - modules/geom.git/blob - src/GEOM_SWIG/GEOM_TestAll.py
Salome HOME
Update copyright information
[modules/geom.git] / src / GEOM_SWIG / GEOM_TestAll.py
1 #  GEOM GEOM_SWIG : binding of C++ omplementaion with Python
2 #
3 #  Copyright (C) 2003  CEA
4 #
5 #  This library is free software; you can redistribute it and/or
6 #  modify it under the terms of the GNU Lesser General Public
7 #  License as published by the Free Software Foundation; either
8 #  version 2.1 of the License.
9 #
10 #  This library is distributed in the hope that it will be useful,
11 #  but WITHOUT ANY WARRANTY; without even the implied warranty of
12 #  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
13 #  Lesser General Public License for more details.
14 #
15 #  You should have received a copy of the GNU Lesser General Public
16 #  License along with this library; if not, write to the Free Software
17 #  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
18 #
19 #  See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
20 #
21 #
22 #  File   : GEOM_usinggeom.py
23 #  Author : Damien COQUERET, Open CASCADE
24 #  Module : GEOM
25 #  $Header$
26
27 # ! Please, if you edit this example file, update also
28 # ! GEOM_SRC/doc/salome/gui/GEOM/input/tui_test_all.doc
29 # ! as some sequences of symbols from this example are used during
30 # ! documentation generation to identify certain places of this file
31
32 def TestAll (geompy, math):
33
34   #Create base Variables
35   nbtimes1  = 5      #Short
36   nbtimes2  = 5
37   mindeg = 2
38   maxdeg = 5
39   nbiter = 5
40   ShapeTypeFace = geompy.ShapeType["FACE"]
41   ShapeTypeEdge = geompy.ShapeType["EDGE"]
42   WantPlanarFace = 1 #True
43
44   radius  = 10.  #Double
45   radius1 = 100.
46   radius2 = 200.
47   height  = 200.
48   d1      = 10.
49   d2      = 10.
50   step1   = 250.
51   step2   = 250.
52   angle   = 45.
53   angle1  = angle * math.pi / 180
54   angle2  = 2 * angle1
55   factor  = 2.
56   tol3d   = 0.0001
57   tol2d   = 0.0001
58   weight  = 1000000.
59   waterdensity = 1.
60   meshingdeflection = 0.01
61   trimsize  = 1000.
62   precision = 0.00001
63
64   #Create base points
65   p0   = geompy.MakeVertex(0.  , 0.  , 0.  ) #(3 Doubles)->GEOM_Object_ptr
66   px   = geompy.MakeVertex(100., 0.  , 0.  )
67   py   = geompy.MakeVertex(0.  , 100., 0.  )
68   pz   = geompy.MakeVertex(0.  , 0.  , 100.)
69   pxyz = geompy.MakeVertex(100., 100., 100.)
70
71   p200 = geompy.MakeVertexWithRef(pxyz, 100., 100., 100.) #(GEOM_Object_ptr, 3 Doubles)->GEOM_Object_ptr
72
73   #Create base directions
74   vx   = geompy.MakeVector(p0, px) #(GEOM_Object_ptr, GEOM_Object_ptr)->GEOM_Object_ptr
75   vy   = geompy.MakeVector(p0, py)
76   vz   = geompy.MakeVector(p0, pz)
77   vxy  = geompy.MakeVector(px, py)
78
79   vxyz = geompy.MakeVectorDXDYDZ(100., 100., 100.) #(3 Doubles)->GEOM_Object_ptr
80
81   #Create local coordinate systems
82   cs1 = geompy.MakeMarker(50,50,50, 1,0,0, 0,1,0) #(9 Doubles)->GEOM_Object_ptr
83   cs2 = geompy.MakeMarker(70,80,10, 1,0,1, 1,1,0) #(9 Doubles)->GEOM_Object_ptr
84   cs3 = geompy.MakeMarkerPntTwoVec(pz, vxy, vz)   #(3 GEOM_Object_ptr)->GEOM_Object_ptr
85
86   #Create base geometry 2D
87   Line   = geompy.MakeLineTwoPnt(p0, pxyz)                  #(2 GEOM_Object_ptr)->GEOM_Object_ptr
88   Line1  = geompy.MakeLine(pz, vxy)                         #(2 GEOM_Object_ptr)->GEOM_Object_ptr
89   Line2  = geompy.MakeLineTwoPnt(pxyz, pz)                  #(2 GEOM_Object_ptr)->GEOM_Object_ptr
90   Plane  = geompy.MakePlane(pz, vxyz, trimsize)             #(2 GEOM_Object_ptr, Double)->GEOM_Object_ptr
91   Plane1 = geompy.MakePlaneThreePnt(px, pz, p200, trimsize) #(4 Doubles)->GEOM_Object_ptr
92
93   Arc      = geompy.MakeArc(py, pz, px)                   #(3 GEOM_Object_ptr)->GEOM_Object_ptr
94   Arc2     = geompy.MakeArcCenter(py, pz, px,0)           #(3 GEOM_Object_ptr,Boolean)->GEOM_Object_ptr
95   Circle   = geompy.MakeCircle(p0, vz, radius1)           #(2 GEOM_Object_ptr, Double)->GEOM_Object_ptr
96   Circle1  = geompy.MakeCircleThreePnt(p0, pxyz, px)      #(3 GEOM_Object_ptr)->GEOM_Object_ptr
97   Circle2  = geompy.MakeCircleCenter2Pnt(p0, pxyz, py)    #(3 GEOM_Object_ptr)->GEOM_Object_ptr
98   Ellipse  = geompy.MakeEllipse(p0, vy, radius2, radius1) #(2 GEOM_Object_ptr, 2 Doubles)->GEOM_Object_ptr
99   Polyline = geompy.MakePolyline([p0, pz, py, p200])      #(List of GEOM_Object_ptr)->GEOM_Object_ptr
100   Bezier   = geompy.MakeBezier([p0, pz, p200, px])        #(List of GEOM_Object_ptr)->GEOM_Object_ptr
101   Interpol = geompy.MakeInterpol([px, py, p200, pxyz])    #(List of GEOM_Object_ptr)->GEOM_Object_ptr
102   Sketcher = geompy.MakeSketcher("Sketcher:F -100 -100:TT 250 -100:R 0:C 100 150:R 0:L 300:WW",
103                                  [100,0,0, 1,1,1, -1,1,0]) #(String, List of Doubles)->GEOM_Object_ptr
104
105   #Test point on curve creation
106   p_on_arc = geompy.MakeVertexOnCurve(Arc, 0.25) #(GEOM_Object_ptr, Double)->GEOM_Object_ptr
107
108   #Test point on lines intersection
109   p_on_l1l2 = geompy.MakeVertexOnLinesIntersection(Line1, Line2) #(2 GEOM_Object_ptr)->GEOM_Object_ptr
110
111   #Test tangent on curve creation
112   tan_on_arc = geompy.MakeTangentOnCurve(Arc, 0.7) #(GEOM_Object_ptr, Double)->GEOM_Object_ptr
113
114   #Create base geometry 3D
115   Box      = geompy.MakeBoxTwoPnt(p0, p200)                   #(2 GEOM_Object_ptr)->GEOM_Object_ptr
116   Box1     = geompy.MakeBoxDXDYDZ(10, 20, 30)                 #(3 Doubles)->GEOM_Object_ptr
117   Box2     = geompy.MakeBox(10,20,30, 15,25,35)               #(6 Doubles)->GEOM_Object_ptr
118   Cylinder = geompy.MakeCylinder(p0, vz, radius1, height)     #(2 GEOM_Object_ptr, 2 Doubles)->GEOM_Object_ptr
119   Cyl1     = geompy.MakeCylinderRH(radius2, height)           #(2 Doubles)->GEOM_Object_ptr
120   Sphere   = geompy.MakeSpherePntR(p0, radius1)               #(GEOM_Object_ptr, Double)->GEOM_Object_ptr
121   Sphere1  = geompy.MakeSphereR(radius)                       #(Double)->GEOM_Object_ptr
122   Sphere2  = geompy.MakeSphere(50, 70, 30, radius)            #(4 Doubles)->GEOM_Object_ptr
123   Cone     = geompy.MakeCone(p0, vz, radius2, radius, height) #(2 GEOM_Object_ptr, 3 Doubles)->GEOM_Object_ptr
124   Cone1    = geompy.MakeConeR1R2H(radius1, radius, height)    #(3 Doubles)->GEOM_Object_ptr
125   Torus    = geompy.MakeTorus(p0, vz, radius2, radius)        #(2 GEOM_Object_ptr, 2 Doubles)->GEOM_Object_ptr
126   Torus1   = geompy.MakeTorusRR(radius2, radius1)             #(2 Doubles)->GEOM_Object_ptr
127
128   #Boolean (Common, Cut, Fuse, Section)
129   Common  = geompy.MakeBoolean(Box, Sphere, 1) #(2 GEOM_Object_ptr, Short)->GEOM_Object_ptr
130   Cut     = geompy.MakeBoolean(Box, Sphere, 2)
131   Fuse    = geompy.MakeBoolean(Box, Sphere, 3)
132   Section = geompy.MakeBoolean(Box, Sphere, 4)
133
134   #Create base objects
135   Edge     = geompy.MakeEdge(p0, pxyz)               #(2 GEOM_Object_ptr)->GEOM_Object_ptr
136   Wire     = geompy.MakeWire([vxy, Arc])             #(List Of GEOM_Object_ptr)->GEOM_Object_ptr
137   Face     = geompy.MakeFace(Wire, WantPlanarFace)   #(GEOM_Object_ptr, Boolean)->GEOM_Object_ptr
138   Face1    = geompy.MakeFaceWires([Wire, Sketcher],
139                                   WantPlanarFace)    #(List of GEOM_Object_ptr, Boolean)->GEOM_Object_ptr
140   Face2    = geompy.MakeFace(Sketcher, WantPlanarFace)
141   Shell    = geompy.MakeShell([Face, Face1])         #(List of GEOM_Object_ptr)->GEOM_Object_ptr
142
143   Prism1   = geompy.MakePrism(Face2, p0, pxyz)       #(3 GEOM_Object_ptr)->GEOM_Object_ptr
144   prism1_faces = geompy.SubShapeAllSorted(Prism1, ShapeTypeFace)
145   Shell1   = geompy.MakeShell([prism1_faces[0], prism1_faces[1],
146                                prism1_faces[3], prism1_faces[4],
147                                prism1_faces[5], prism1_faces[2]])
148   Solid    = geompy.MakeSolid([Shell1])              #(List of GEOM_Object_ptr)->GEOM_Object_ptr
149
150   ShapeListCompound = []
151   i = 0
152   while i <= 3 :
153         S = geompy.MakeTranslation(Arc, i * 100., i * 100., i * 100.)
154         ShapeListCompound.append(S)
155         i = i + 1
156   Compound = geompy.MakeCompound(ShapeListCompound)  #(List of GEOM_Object_ptr)->GEOM_Object_ptr
157
158   #Test point on surface creation
159   p_on_face = geompy.MakeVertexOnSurface(Face, 0.1, 0.8) #(GEOM_Object_ptr, Double, Double)->GEOM_Object_ptr
160
161   # Test plane from existing face creation
162   Plane2 = geompy.MakePlaneFace(Face, trimsize)      #(GEOM_Object_ptr, Double)->GEOM_Object_ptr
163
164   #ShapeList for Sewing
165   S = geompy.MakeRotation(Face, vxy, angle1)
166
167   #Test Line on Faces Intersection
168   Line3 = geompy.MakeLineTwoFaces(prism1_faces[0], prism1_faces[1]) #(2 GEOM_Object_ptr)->GEOM_Object_ptr
169
170   #Create advanced objects
171   Copy       = geompy.MakeCopy(Box)                      #(GEOM_Object_ptr)->GEOM_Object_ptr
172   Prism      = geompy.MakePrismVecH(Face, vz, 100.0)     #(2 GEOM_Object_ptr, Double)->GEOM_Object_ptr
173   Prism2Ways = geompy.MakePrismVecH2Ways(Face, vz, 10.0) #(2 GEOM_Object_ptr, Double)->GEOM_Object_ptr
174   Revolution = geompy.MakeRevolution(Face, vz, angle2)   #
175   Revolution2Ways = geompy.MakeRevolution(Face, vz, angle1)   #
176   Filling    = geompy.MakeFilling(Compound, mindeg, maxdeg,
177                                   tol2d, tol3d, nbiter)  #(GEOM_Object_ptr, 4 Doubles, Short)->GEOM_Object_ptr
178   Pipe       = geompy.MakePipe(Wire, Edge)               #(2 GEOM_Object_ptr)->GEOM_Object_ptr
179   Sewing     = geompy.MakeSewing([Face, S], precision)   #(List Of GEOM_Object_ptr, Double)->GEOM_Object_ptr
180
181   #Transform objects
182   Translation = geompy.MakeTranslationTwoPoints(Box, px, pz)    #(3 GEOM_Object_ptr)->GEOM_Object_ptr
183   TranslVect  = geompy.MakeTranslationVector(Box, vxyz)         #(2 GEOM_Object_ptr)->GEOM_Object_ptr
184   TranslVectD = geompy.MakeTranslationVectorDistance(Box, vxyz, 50.0)   #(2 GEOM_Object_ptr)->GEOM_Object_ptr  
185   Rotation    = geompy.MakeRotation(Box, vz, angle1)            #(2 GEOM_Object_ptr, Double)->GEOM_Object_ptr
186   RotatPnt    = geompy.MakeRotationThreePoints(Box, px, py, pz) #(4 GEOM_Object_ptr)->GEOM_Object_ptr
187
188   #Scale by factor relatively given point
189   Scale1 = geompy.MakeScaleTransform(Box, pxyz, factor)      #(2 GEOM_Object_ptr, Double)->GEOM_Object_ptr
190   #Scale by factor relatively the origin of global CS
191   Scale2 = geompy.MakeScaleTransform(Box, None, factor)      #
192   #Scale along axes of global CS by different factors. Scale relatively given point
193   Scale3 = geompy.MakeScaleAlongAxes(Box, pxyz, 1.5, 0.5, 3) #(2 GEOM_Object_ptr, 3 Doubles)->GEOM_Object_ptr
194   #Scale along axes of global CS by different factors. Scale relatively the origin of global CS
195   Scale4 = geompy.MakeScaleAlongAxes(Box, None, 1.5, 0.5, 3) #
196
197   Mirror      = geompy.MakeMirrorByPlane(Box, Plane) #(2 GEOM_Object_ptr)->GEOM_Object_ptr
198   MirrorAxis  = geompy.MakeMirrorByAxis(Box, Line1)  #
199   MirrorPnt   = geompy.MakeMirrorByPoint(Box, p200)  #
200   Position    = geompy.MakePosition(Box, cs1, cs2)   #(3 GEOM_Object_ptr)->GEOM_Object_ptr
201   Offset      = geompy.MakeOffset(Box, 10.)          #(GEOM_Object_ptr, Double)->GEOM_Object_ptr
202   Orientation = geompy.ChangeOrientation(Box)
203
204   #IDList for Fillet/Chamfer
205   prism_edges = geompy.SubShapeAllSorted(Prism, ShapeTypeEdge)
206
207   for anEdge in prism_edges:
208     eid = geompy.GetSubShapeID(Prism, anEdge)
209     sse = geompy.GetSubShape(Prism, [eid])
210
211     sse_id = geompy.GetSubShapeID(Prism, sse)
212     if sse_id != eid:
213       print "Error: GetSubShape() or GetSubShapeID() has failed!"
214
215   IDlist_e = []
216   IDlist_e.append(geompy.GetSubShapeID(Prism, prism_edges[0]))
217   IDlist_e.append(geompy.GetSubShapeID(Prism, prism_edges[1]))
218   IDlist_e.append(geompy.GetSubShapeID(Prism, prism_edges[2]))
219
220   prism_faces = geompy.SubShapeAllSorted(Prism, ShapeTypeFace)
221
222   f_ind_1 = geompy.GetSubShapeID(Prism, prism_faces[0])
223   f_ind_2 = geompy.GetSubShapeID(Prism, prism_faces[1])
224
225   IDlist_f = [f_ind_1, f_ind_2]
226   
227   #Local operations
228   Fillet   = geompy.MakeFillet (Prism, radius, ShapeTypeEdge,
229                                 IDlist_e) #(GEOM_Object_ptr, Double, Short, ListOfLong)->GEOM_Object_ptr
230   Fillet2  = geompy.MakeFilletR1R2 (Prism, 7., 13., ShapeTypeEdge,
231                                     IDlist_e) #(GEOM_Object_ptr, Double, Double, Short, ListOfLong)->GEOM_Object_ptr
232   Chamfer  = geompy.MakeChamferEdge(Prism, d1, d2,
233                                     f_ind_1, f_ind_2) #(GEOM_Object_ptr, 2 Doubles, 2 Long)->GEOM_Object_ptr
234   Chamfer2 = geompy.MakeChamferFaces(Prism, d1, d2,
235                                      IDlist_f) #(GEOM_Object_ptr, 2 Doubles, ListOfLong)->GEOM_Object_ptr
236   Chamfer3 = geompy.MakeChamferEdges(Prism, d1, d2,
237                                      IDlist_e) #(GEOM_Object_ptr, 2 Doubles, ListOfLong)->GEOM_Object_ptr
238   Chamfer4 = geompy.MakeChamferFacesAD(Prism, d1, 20. * math.pi / 180.,
239                                        IDlist_f) #(GEOM_Object_ptr, 2 Doubles, ListOfLong)->GEOM_Object_ptr
240   #End of Local operations
241
242   #Create Patterns
243   MultiTrans1D = geompy.MakeMultiTranslation1D(Fillet, vz, step1, nbtimes1)
244   MultiTrans2D = geompy.MakeMultiTranslation2D(Fillet, vz, step1, nbtimes1, vy, step2, nbtimes2)
245   #!!!!Angle In Degree!!!!
246   MultiRot1D   = geompy.MultiRotate1D(Chamfer, vx, nbtimes1)
247   MultiRot2D   = geompy.MultiRotate2D(Chamfer, vx, angle, nbtimes1, step1, nbtimes2)
248
249   #Create Informations objects
250   CDG        = geompy.MakeCDG(Prism)               #(GEOM_Object_ptr)->GEOM_Object_ptr
251   Archimede  = geompy.Archimede(Box, weight, waterdensity,
252                                 meshingdeflection) #(GEOM_Object_ptr, 3 Doubles)->GEOM_Object_ptr
253   mindist = geompy.MinDistanceComponents(TranslVect, Mirror) #(2 GEOM_Object_ptr)->4 Doubles
254   print "Minumal distance between TranslVect and Mirror is", mindist[0],
255   print "by components:", mindist[1], ",", mindist[2], ",", mindist[3]
256   CheckShape = geompy.CheckShape(Prism)            #(GEOM_Object_ptr)->Boolean
257   print "CheckShape(Prism) = ", CheckShape
258
259   #Partition objects
260   Partition  = geompy.MakePartition([Box], [Plane]) #(2 Lists Of GEOM_Object_ptr)->GEOM_Object_ptr
261   Partition1 = geompy.MakeHalfPartition(Box, Plane) #(2 GEOM_Object_ptr)->GEOM_Object_ptr
262
263   #Add In Study
264   id_p0   = geompy.addToStudy(p0,   "Vertex 0")
265   id_px   = geompy.addToStudy(px,   "Vertex X")
266   id_py   = geompy.addToStudy(py,   "Vertex Y")
267   id_pz   = geompy.addToStudy(pz,   "Vertex Z")
268   id_pxyz = geompy.addToStudy(pxyz, "Vertex XYZ")
269   id_p200 = geompy.addToStudy(p200, "Vertex 200")
270
271   id_vx   = geompy.addToStudy(vx,   "Vector X")
272   id_vy   = geompy.addToStudy(vy,   "Vector Y")
273   id_vz   = geompy.addToStudy(vz,   "Vector Z")
274   id_vxy  = geompy.addToStudy(vxy,  "Vector XY")
275   id_vxyz = geompy.addToStudy(vxyz, "Vector XYZ")
276
277   id_cs1 = geompy.addToStudy(cs1, "CS 50,50,50, 1,0,0, 0,1,0")
278   id_cs2 = geompy.addToStudy(cs2, "CS 70,80,10, 1,0,1, 1,1,0")
279   id_cs3 = geompy.addToStudy(cs3, "CS: pz, vxy, vz")
280
281   id_Line   = geompy.addToStudy(Line,   "Line")
282   id_Line1  = geompy.addToStudy(Line1,  "Line by point and vector")
283   id_Line3  = geompy.addToStudy(Line3,  "Line on Two Faces Intersection")
284   id_Plane  = geompy.addToStudy(Plane,  "Plane")
285   id_Plane1 = geompy.addToStudy(Plane1,  "Plane by 3 points")
286
287   id_Arc      = geompy.addToStudy(Arc,      "Arc")
288   id_Circle   = geompy.addToStudy(Circle,   "Circle")
289   id_Circle1  = geompy.addToStudy(Circle1,  "Circle by 3 points")
290   id_Circle2  = geompy.addToStudy(Circle2,  "Circle by center and 2 points")
291   id_Ellipse  = geompy.addToStudy(Ellipse,  "Ellipse")
292   id_Polyline = geompy.addToStudy(Polyline, "Polyline")
293   id_Bezier   = geompy.addToStudy(Bezier,   "Bezier")
294   id_Interpol = geompy.addToStudy(Interpol, "Interpol")
295   id_Sketcher = geompy.addToStudy(Sketcher, "Sketcher")
296
297   id_p_on_arc  = geompy.addToStudy(p_on_arc,  "Vertex on Arc (0.25)")
298   id_p_on_l1l2 = geompy.addToStudy(p_on_l1l2, "Vertex on Lines Intersection")
299
300   id_tan_on_arc = geompy.addToStudy(tan_on_arc, "Tangent on Arc (0.7)")
301
302   id_Box      = geompy.addToStudy(Box,      "Box")
303   id_Box1     = geompy.addToStudy(Box1,     "Box 10x20x30")
304   id_Box2     = geompy.addToStudy(Box2,     "Box (10,20,30)-(15,25,35)")
305   id_Cylinder = geompy.addToStudy(Cylinder, "Cylinder")
306   id_Cyl1     = geompy.addToStudy(Cyl1,     "Cylinder RH")
307   id_Sphere   = geompy.addToStudy(Sphere,   "Sphere Pnt R")
308   id_Sphere1  = geompy.addToStudy(Sphere1,  "Sphere R")
309   id_Sphere2  = geompy.addToStudy(Sphere2,  "Sphere")
310   id_Cone     = geompy.addToStudy(Cone,     "Cone")
311   id_Cone1    = geompy.addToStudy(Cone1,    "Cone R1R2H")
312   id_Torus    = geompy.addToStudy(Torus,    "Torus")
313   id_Torus1   = geompy.addToStudy(Torus1,   "Torus RR")
314
315   id_Common  = geompy.addToStudy(Common,  "Common")
316   id_Cut     = geompy.addToStudy(Cut,     "Cut")
317   id_Fuse    = geompy.addToStudy(Fuse,    "Fuse")
318   id_Section = geompy.addToStudy(Section, "Section")
319
320   id_Edge     = geompy.addToStudy(Edge,     "Edge")
321   id_Wire     = geompy.addToStudy(Wire,     "Wire")
322   id_Face     = geompy.addToStudy(Face,     "Face")
323   id_Face1    = geompy.addToStudy(Face1,    "Face from two wires")
324   id_Face2    = geompy.addToStudy(Face2,    "Face from Sketcher")
325   id_Shell    = geompy.addToStudy(Shell,    "Shell")
326
327   id_p_on_face = geompy.addToStudy(p_on_face, "Vertex on Face (0.1, 0.8)")
328
329   id_Prism1   = geompy.addToStudy(Prism1,     "Prism by Two Pnt")
330   id_Shell1   = geompy.addToStudy(Shell1,   "Shell from Prism1 faces")
331   id_Solid    = geompy.addToStudy(Solid,    "Solid")
332   id_Compound = geompy.addToStudy(Compound, "Compound")
333
334   id_Plane2   = geompy.addToStudy(Plane2,   "Plane on Face")
335
336   id_Copy       = geompy.addToStudy(Copy,       "Copy")
337   id_Prism      = geompy.addToStudy(Prism,      "Prism")
338   id_Prism2Ways = geompy.addToStudy(Prism2Ways, "Prism2Ways")
339   id_Revolution = geompy.addToStudy(Revolution, "Revolution")
340   id_Revolution2Ways = geompy.addToStudy(Revolution2Ways, "Revolution2Ways")
341   id_Filling    = geompy.addToStudy(Filling,    "Filling")
342   id_Pipe       = geompy.addToStudy(Pipe,       "Pipe")
343   id_Sewing     = geompy.addToStudy(Sewing,     "Sewing")
344
345   id_Translation = geompy.addToStudy(Translation, "Translation")
346   id_TranslVect  = geompy.addToStudy(TranslVect , "Translation along vector")
347   id_TranslVectD = geompy.addToStudy(TranslVectD, "Translation along vector with defined distance")
348   id_Rotation    = geompy.addToStudy(Rotation,    "Rotation")
349   id_RotatPnt    = geompy.addToStudy(RotatPnt,    "Rotation by three points")
350   id_Scale1      = geompy.addToStudy(Scale1,      "Scale1")
351   id_Scale2      = geompy.addToStudy(Scale2,      "Scale2")
352   id_Scale3      = geompy.addToStudy(Scale3,      "Scale3")
353   id_Scale4      = geompy.addToStudy(Scale4,      "Scale4")
354   id_Mirror      = geompy.addToStudy(Mirror,      "Mirror by Plane")
355   id_MirrorAxis  = geompy.addToStudy(MirrorAxis,  "Mirror by Axis")
356   id_MirrorPnt   = geompy.addToStudy(MirrorPnt,   "Mirror by Point")
357   id_Position    = geompy.addToStudy(Position,    "Positioned box")
358   id_Offset      = geompy.addToStudy(Offset,      "Offset")
359   id_Orientation = geompy.addToStudy(Orientation, "Orientation")
360
361   id_Fillet   = geompy.addToStudy(Fillet,   "Fillet")
362   id_Fillet2  = geompy.addToStudy(Fillet2,  "Fillet2")
363
364   id_Chamfer  = geompy.addToStudy(Chamfer,  "Chamfer on Edge")
365   id_Chamfer2 = geompy.addToStudy(Chamfer2, "Chamfer on Faces")
366   id_Chamfer3 = geompy.addToStudy(Chamfer3, "Chamfer on Edges")
367   id_Chamfer4 = geompy.addToStudy(Chamfer4, "Chamfer on Faces with params D Angle")
368
369   id_MultiTrans1D = geompy.addToStudy(MultiTrans1D, "MultiTrans1D")
370   id_MultiTrans2D = geompy.addToStudy(MultiTrans2D, "MultiTrans2D")
371   id_MultiRot1D   = geompy.addToStudy(MultiRot1D,   "MultiRot1D")
372   id_MultiRot2D   = geompy.addToStudy(MultiRot2D,   "MultiRot2D")
373
374   id_CDG       = geompy.addToStudy(CDG,       "CDG")
375   id_Archimede = geompy.addToStudy(Archimede, "Archimede")
376
377   id_Partition  = geompy.addToStudy(Partition, "Partition")
378   id_Partition1 = geompy.addToStudy(Partition1, "Half Partition")
379
380   #Decompose objects
381
382   # SubShape
383   SubFace    = geompy.SubShape(Box, geompy.ShapeType["FACE"], [2])
384   name       = geompy.SubShapeName(SubFace, Box)
385   id_SubFace = geompy.addToStudyInFather(Box, SubFace, name)
386
387   # SubShapeSorted
388   SubFaceS   = geompy.SubShapeSorted(Box, geompy.ShapeType["FACE"], [5])
389   nameS      = geompy.SubShapeName(SubFaceS, Box)
390   id_SubFace = geompy.addToStudyInFather(Box, SubFaceS, nameS)
391
392   # SubShapeAll
393   SubEdgeList = geompy.SubShapeAll(SubFace, geompy.ShapeType["EDGE"])
394   i=0
395   for SubEdge in SubEdgeList :
396     name = geompy.SubShapeName(SubEdge, SubFace)
397     id_SubEdge = geompy.addToStudyInFather(SubFace, SubEdge, name)
398
399   # SubShapeAllIDs
400   SubEdgeIDsList = geompy.SubShapeAllIDs(SubFace, geompy.ShapeType["EDGE"])
401   print "IDs of edges of SubFace:", SubEdgeIDsList, "(unsorted)"
402   group = geompy.CreateGroup(SubFace, geompy.ShapeType["EDGE"])
403   geompy.UnionIDs(group, SubEdgeIDsList)
404   geompy.addToStudyInFather(SubFace, group, "Group of all edges")
405
406   # SubShapeAllSortedIDs
407   SubEdgeIDsList = geompy.SubShapeAllSortedIDs(SubFace, geompy.ShapeType["EDGE"])
408   print "IDs of edges of SubFace:", SubEdgeIDsList, "(sorted)"
409
410   # GetSubShape and GetSubShapeID
411   for ind in SubEdgeIDsList:
412     edge = geompy.GetSubShape(SubFace, [ind])
413     ind_e = geompy.GetSubShapeID(SubFace, edge)
414     if ind_e != ind:
415       print "Error in GetSubShape or GetSubShapeID"
416
417   print "DONE"