--- /dev/null
+# CEA/LGLS 2004, Francis KLOSS (OCC)
+# ==================================
+
+# Import
+# ------
+
+from geompy import *
+from meshpy import *
+
+# Piece
+# -----
+
+# grid compound of 3 x 3 elements
+# an element is compound of 3 cylinders concentriques
+# an element is centered in a square of the grid
+# the smaller cylinder is a hole
+
+# prisme the grid
+
+# Geometry
+# --------
+
+g_x = 0
+g_y = 0
+g_z = 0
+
+g_arete = 50
+g_hauteur = 30
+
+g_rayon1 = 20
+g_rayon2 = 30
+g_rayon3 = 40
+
+g_grid = 3
+
+g_trim = 1000
+
+# Element
+# -------
+
+e_boite = MakeBox(g_x-g_arete, g_y-g_hauteur, g_z-g_arete, g_x+g_arete, g_y+g_hauteur, g_z+g_arete)
+
+e_pi4 = 3.141592653/4
+e_hauteur = 2*g_hauteur
+e_centre = MakeVertex(g_x, g_y-g_hauteur, g_z)
+e_dir = MakeVectorDXDYDZ(0, 1, 0)
+
+e_cyl0 = MakeCylinder(e_centre, e_dir, g_rayon3, e_hauteur)
+e_cyl1 = MakeRotation(e_cyl0, e_dir, e_pi4)
+
+e_blo1 = MakeCut(e_boite, e_cyl1)
+
+e_cyl0 = MakeCylinder(e_centre, e_dir, g_rayon2, e_hauteur)
+e_cyl2 = MakeRotation(e_cyl0, e_dir, e_pi4)
+
+e_blo2 = MakeCut(e_cyl1, e_cyl2)
+
+e_cyl0 = MakeCylinder(e_centre, e_dir, g_rayon1, e_hauteur)
+e_cyl3 = MakeRotation(e_cyl0, e_dir, e_pi4)
+
+e_blo3 = MakeCut(e_cyl2, e_cyl3)
+
+# Partition
+# ---------
+
+p_tools = []
+p_tools.append(MakePlane(e_centre, MakeVectorDXDYDZ( 1, 0, 1), g_trim))
+p_tools.append(MakePlane(e_centre, MakeVectorDXDYDZ(-1, 0, 1), g_trim))
+
+p_element = MakePartition([e_blo1, e_blo2, e_blo3], p_tools, [], [], ShapeType["SOLID"])
+
+# Grid
+# ----
+
+piece = MakeMultiTranslation2D(p_element, MakeVectorDXDYDZ(1, 0, 0), 2*g_arete, g_grid,
+ MakeVectorDXDYDZ(0, 0, 1), 2*g_arete, g_grid)
+
+# Add in study
+# ------------
+
+piece_id = addToStudy(piece, "Grid3partition")
+
+# Meshing
+# =======
+
+# Create hexahedrical mesh on piece
+# ---------------------------------
+
+m_hexa=MeshHexa(piece, 4, "Grid3partitionHexa")
+
+# Compute
+# -------
+
+m_hexa.Compute()
