From: nadir <nadir>
Date: Thu, 16 Oct 2003 13:14:47 +0000 (+0000)
Subject: first version of this file needed for the Netgen integration in the SMESH
X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=bf3481306b2595978e7d51d7269bb276aec7af42;p=modules%2Fsmesh.git

first version of this file needed for the Netgen integration in the SMESH
Engine.
---

diff --git a/src/SMESH_SWIG/SMESH_box2_tetra.py b/src/SMESH_SWIG/SMESH_box2_tetra.py
new file mode 100644
index 000000000..1baf0be33
--- /dev/null
+++ b/src/SMESH_SWIG/SMESH_box2_tetra.py
@@ -0,0 +1,191 @@
+#
+# Tetrahedrization of the geometry union of 2 boxes having a face in common
+# Hypothesis and algorithms for the mesh generation are global
+#
+
+import salome
+from salome import sg
+
+import geompy
+
+import SMESH
+import smeshpy
+
+##geom = salome.lcc.FindOrLoadComponent("FactoryServer", "Geometry")
+##myBuilder = salome.myStudy.NewBuilder()
+
+geom = geompy.geom
+myBuilder = geompy.myBuilder
+
+ShapeTypeShell     = 3
+ShapeTypeFace      = 4
+ShapeTypeEdge      = 6
+
+# ---- define 2 boxes box1 and box2
+
+box1 = geompy.MakeBox(0., 0., 0., 100., 200., 300.)
+
+idbox1 = geompy.addToStudy(box1,"box1")
+
+print "Analysis of the geometry box1 :"
+subShellList=geompy.SubShapeAll(box1,ShapeTypeShell)
+subFaceList=geompy.SubShapeAll(box1,ShapeTypeFace)
+subEdgeList=geompy.SubShapeAll(box1,ShapeTypeEdge)
+
+print "number of Shells in box1 : ",len(subShellList)
+print "number of Faces in box1 : ",len(subFaceList)
+print "number of Edges in box1 : ",len(subEdgeList)
+
+box2 = geompy.MakeBox(100., 0., 0., 200., 200., 300.)
+
+idbox2 = geompy.addToStudy(box2,"box2")
+
+print "Analysis of the geometry box2 :"
+subShellList=geompy.SubShapeAll(box2,ShapeTypeShell)
+subFaceList=geompy.SubShapeAll(box2,ShapeTypeFace)
+subEdgeList=geompy.SubShapeAll(box2,ShapeTypeEdge)
+
+print "number of Shells in box2 : ",len(subShellList)
+print "number of Faces in box2 : ",len(subFaceList)
+print "number of Edges in box2 : ",len(subEdgeList)
+
+blocs = []
+blocs.append(box1._get_Name())
+blocs.append(box2._get_Name())
+
+tol3d = 1.e-5
+
+shell = geom.MakeSewing(blocs,tol3d)
+idshell = geompy.addToStudy(shell,"shell")
+
+print "Analysis of the geometry shell (union of box1 and box2) :"
+subShellList=geompy.SubShapeAll(shell,ShapeTypeShell)
+subFaceList=geompy.SubShapeAll(shell,ShapeTypeFace)
+subEdgeList=geompy.SubShapeAll(shell,ShapeTypeEdge)
+
+print "number of Shells in shell : ",len(subShellList)
+print "number of Faces in shell : ",len(subFaceList)
+print "number of Edges in shell : ",len(subEdgeList)
+
+# ---- launch SMESH
+
+smeshgui = salome.ImportComponentGUI("SMESH")
+smeshgui.Init(salome.myStudyId)
+
+gen=smeshpy.smeshpy()
+
+# ---- create Hypothesis
+
+print "-------------------------- create Hypothesis"
+
+print "-------------------------- NumberOfSegments"
+
+numberOfSegments = 10
+
+hyp1=gen.CreateHypothesis("NumberOfSegments")
+hypNbSeg=hyp1._narrow(SMESH.SMESH_NumberOfSegments)
+hypNbSeg.SetNumberOfSegments(numberOfSegments)
+hypNbSegID = hypNbSeg.GetId()
+print hypNbSeg.GetName()
+print hypNbSegID
+print hypNbSeg.GetNumberOfSegments()
+
+idseg = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypNbSeg) )
+smeshgui.SetName(idseg, "NumberOfSegments")
+
+print "-------------------------- MaxElementArea"
+
+maxElementArea = 500
+
+hyp2=gen.CreateHypothesis("MaxElementArea")
+hypArea=hyp2._narrow(SMESH.SMESH_MaxElementArea)
+hypArea.SetMaxElementArea(maxElementArea)
+print hypArea.GetName()
+print hypArea.GetId()
+print hypArea.GetMaxElementArea()
+
+idarea = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypArea) )
+smeshgui.SetName(idarea, "MaxElementArea")
+
+print "-------------------------- MaxElementVolume"
+
+maxElementVolume = 500
+
+hyp3=gen.CreateHypothesis("MaxElementVolume")
+hypVolume=hyp3._narrow(SMESH.SMESH_MaxElementVolume)
+hypVolume.SetMaxElementVolume(maxElementVolume)
+print hypVolume.GetName()
+print hypVolume.GetId()
+print hypVolume.GetMaxElementVolume()
+
+idvolume = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypVolume) )
+smeshgui.SetName(idvolume, "MaxElementVolume")
+
+# ---- create Algorithms
+
+print "-------------------------- create Algorithms"
+
+print "-------------------------- Regular_1D"
+
+hypothesis=gen.CreateHypothesis("Regular_1D")
+regular1D = hypothesis._narrow(SMESH.SMESH_Regular_1D)
+regularID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(regular1D) )
+smeshgui.SetName(regularID, "Wire Discretisation")
+
+print "-------------------------- MEFISTO_2D"
+
+hypothesis=gen.CreateHypothesis("MEFISTO_2D")
+mefisto2D = hypothesis._narrow(SMESH.SMESH_MEFISTO_2D)
+mefistoID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(mefisto2D) )
+smeshgui.SetName(mefistoID, "MEFISTO_2D")
+
+print "-------------------------- Tetra_3D"
+
+hypothesis=gen.CreateHypothesis("Tetra_3D")
+tetra3D = hypothesis._narrow(SMESH.SMESH_Tetra_3D)
+tetraID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(tetra3D) )
+smeshgui.SetName(tetraID, "Tetra_3D")
+
+# ---- init a Mesh with the shell
+
+mesh=gen.Init(idshell)
+idmesh = smeshgui.AddNewMesh( salome.orb.object_to_string(mesh) )
+smeshgui.SetName(idmesh, "MeshBox2")
+smeshgui.SetShape(idshell, idmesh)
+
+# ---- add hypothesis to shell
+
+print "-------------------------- add hypothesis to shell"
+
+ret=mesh.AddHypothesis(shell,regular1D)
+print ret
+ret=mesh.AddHypothesis(shell,hypNbSeg)
+print ret
+ret=mesh.AddHypothesis(shell,mefisto2D)
+print ret
+ret=mesh.AddHypothesis(shell,hypArea)
+print ret
+ret=mesh.AddHypothesis(shell,tetra3D)
+print ret
+ret=mesh.AddHypothesis(shell,hypVolume)
+print ret
+
+smeshgui.SetAlgorithms( idmesh, regularID)
+smeshgui.SetHypothesis( idmesh, idseg )
+smeshgui.SetAlgorithms( idmesh, mefistoID )
+smeshgui.SetHypothesis( idmesh, idarea )
+smeshgui.SetAlgorithms( idmesh, tetraID )
+smeshgui.SetHypothesis( idmesh, idvolume )
+
+sg.updateObjBrowser(1)
+
+
+print "-------------------------- compute shell"
+ret=gen.Compute(mesh,idshell)
+print ret
+log=mesh.GetLog(0) # no erase trace
+for linelog in log:
+    print linelog
+
+
+sg.updateObjBrowser(1)
diff --git a/src/SMESH_SWIG/SMESH_box_tetra.py b/src/SMESH_SWIG/SMESH_box_tetra.py
new file mode 100644
index 000000000..3a49053d7
--- /dev/null
+++ b/src/SMESH_SWIG/SMESH_box_tetra.py
@@ -0,0 +1,162 @@
+#
+# Tetrahedrization of a simple box. Hypothesis and algorithms for
+# the mesh generation are global
+#
+
+import salome
+from salome import sg
+
+import geompy
+
+import SMESH
+import smeshpy
+
+# -----------------------------------------------------------------------------
+
+##geom = salome.lcc.FindOrLoadComponent("FactoryServer", "Geometry")
+##myBuilder = salome.myStudy.NewBuilder()
+
+geom = geompy.geom
+myBuilder = geompy.myBuilder
+
+ShapeTypeShell     = 3
+ShapeTypeFace      = 4
+ShapeTypeEdge      = 6
+
+# ---- define a boxe
+
+box = geompy.MakeBox(0., 0., 0., 100., 200., 300.)
+
+idbox = geompy.addToStudy(box,"box")
+
+print "Analysis of the geometry box :"
+subShellList=geompy.SubShapeAll(box,ShapeTypeShell)
+subFaceList=geompy.SubShapeAll(box,ShapeTypeFace)
+subEdgeList=geompy.SubShapeAll(box,ShapeTypeEdge)
+
+print "number of Shells in box : ",len(subShellList)
+print "number of Faces in box : ",len(subFaceList)
+print "number of Edges in box : ",len(subEdgeList)
+
+# ---- launch SMESH
+
+smeshgui = salome.ImportComponentGUI("SMESH")
+smeshgui.Init(salome.myStudyId)
+
+gen=smeshpy.smeshpy()
+
+# ---- create Hypothesis
+
+print "-------------------------- create Hypothesis"
+
+print "-------------------------- NumberOfSegments"
+
+numberOfSegments = 10
+
+hyp1=gen.CreateHypothesis("NumberOfSegments")
+hypNbSeg=hyp1._narrow(SMESH.SMESH_NumberOfSegments)
+hypNbSeg.SetNumberOfSegments(numberOfSegments)
+hypNbSegID = hypNbSeg.GetId()
+print hypNbSeg.GetName()
+print hypNbSegID
+print hypNbSeg.GetNumberOfSegments()
+
+idseg = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypNbSeg) )
+smeshgui.SetName(idseg, "NumberOfSegments")
+
+print "-------------------------- MaxElementArea"
+
+maxElementArea = 500
+
+hyp2=gen.CreateHypothesis("MaxElementArea")
+hypArea=hyp2._narrow(SMESH.SMESH_MaxElementArea)
+hypArea.SetMaxElementArea(maxElementArea)
+print hypArea.GetName()
+print hypArea.GetId()
+print hypArea.GetMaxElementArea()
+
+idarea = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypArea) )
+smeshgui.SetName(idarea, "MaxElementArea")
+
+print "-------------------------- MaxElementVolume"
+
+maxElementVolume = 500
+
+hyp3=gen.CreateHypothesis("MaxElementVolume")
+hypVolume=hyp3._narrow(SMESH.SMESH_MaxElementVolume)
+hypVolume.SetMaxElementVolume(maxElementVolume)
+print hypVolume.GetName()
+print hypVolume.GetId()
+print hypVolume.GetMaxElementVolume()
+
+idvolume = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypVolume) )
+smeshgui.SetName(idvolume, "MaxElementVolume")
+
+# ---- create Algorithms
+
+print "-------------------------- create Algorithms"
+
+print "-------------------------- Regular_1D"
+
+hypothesis=gen.CreateHypothesis("Regular_1D")
+regular1D = hypothesis._narrow(SMESH.SMESH_Regular_1D)
+regularID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(regular1D) )
+smeshgui.SetName(regularID, "Wire Discretisation")
+
+print "-------------------------- MEFISTO_2D"
+
+hypothesis=gen.CreateHypothesis("MEFISTO_2D")
+mefisto2D = hypothesis._narrow(SMESH.SMESH_MEFISTO_2D)
+mefistoID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(mefisto2D) )
+smeshgui.SetName(mefistoID, "MEFISTO_2D")
+
+print "-------------------------- Tetra_3D"
+
+hypothesis=gen.CreateHypothesis("Tetra_3D")
+tetra3D = hypothesis._narrow(SMESH.SMESH_Tetra_3D)
+tetraID = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(tetra3D) )
+smeshgui.SetName(tetraID, "Tetra_3D")
+
+# ---- init a Mesh with the boxe
+
+mesh=gen.Init(idbox)
+idmesh = smeshgui.AddNewMesh( salome.orb.object_to_string(mesh) )
+smeshgui.SetName(idmesh, "MeshBox")
+smeshgui.SetShape(idbox, idmesh)
+
+# ---- add hypothesis to the boxe
+
+print "-------------------------- add hypothesis to the boxe"
+
+ret=mesh.AddHypothesis(box,regular1D)
+print ret
+ret=mesh.AddHypothesis(box,hypNbSeg)
+print ret
+ret=mesh.AddHypothesis(box,mefisto2D)
+print ret
+ret=mesh.AddHypothesis(box,hypArea)
+print ret
+ret=mesh.AddHypothesis(box,tetra3D)
+print ret
+ret=mesh.AddHypothesis(box,hypVolume)
+print ret
+
+smeshgui.SetAlgorithms( idmesh, regularID)
+smeshgui.SetHypothesis( idmesh, idseg )
+smeshgui.SetAlgorithms( idmesh, mefistoID )
+smeshgui.SetHypothesis( idmesh, idarea )
+smeshgui.SetAlgorithms( idmesh, tetraID )
+smeshgui.SetHypothesis( idmesh, idvolume )
+
+sg.updateObjBrowser(1)
+
+
+print "-------------------------- compute the mesh of the boxe"
+ret=gen.Compute(mesh,idbox)
+print ret
+log=mesh.GetLog(0) # no erase trace
+for linelog in log:
+    print linelog
+
+
+sg.updateObjBrowser(1)
diff --git a/src/SMESH_SWIG/SMESH_mechanic_tetra.py b/src/SMESH_SWIG/SMESH_mechanic_tetra.py
new file mode 100644
index 000000000..b83c66472
--- /dev/null
+++ b/src/SMESH_SWIG/SMESH_mechanic_tetra.py
@@ -0,0 +1,269 @@
+#  Copyright (C) 2003  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
+#  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS 
+# 
+#  This library is free software; you can redistribute it and/or 
+#  modify it under the terms of the GNU Lesser General Public 
+#  License as published by the Free Software Foundation; either 
+#  version 2.1 of the License. 
+# 
+#  This library is distributed in the hope that it will be useful, 
+#  but WITHOUT ANY WARRANTY; without even the implied warranty of 
+#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU 
+#  Lesser General Public License for more details. 
+# 
+#  You should have received a copy of the GNU Lesser General Public 
+#  License along with this library; if not, write to the Free Software 
+#  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA 
+# 
+#  See http://www.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org 
+#
+#
+#
+#  File   : SMESH_withHole.py
+#  Author : Lucien PIGNOLONI
+#  Module : SMESH
+#  $Header$
+
+import SMESH
+import smeshpy
+import salome
+from salome import sg
+import math
+#import SMESH_BasicHypothesis_idl
+
+import geompy
+
+
+
+# ---------------------------- GEOM --------------------------------------
+geom = salome.lcc.FindOrLoadComponent("FactoryServer", "GEOM")
+myBuilder = salome.myStudy.NewBuilder()
+#from geompy import gg
+
+smeshgui = salome.ImportComponentGUI("SMESH")
+smeshgui.Init(salome.myStudyId)
+
+ShapeTypeCompSolid = 1
+ShapeTypeSolid     = 2
+ShapeTypeShell     = 3
+ShapeTypeFace      = 4
+ShapeTypeWire      = 5
+ShapeTypeEdge      = 6
+ShapeTypeVertex    = 7
+
+
+# ---- define contigous arcs and segment to define a closed wire
+
+p1   = geom.MakePointStruct( 100.0,   0.0,  0.0 )
+p2   = geom.MakePointStruct(  50.0,  50.0,  0.0 )
+p3   = geom.MakePointStruct( 100.0, 100.0,  0.0 ) 
+arc1 = geom.MakeArc( p1, p2, p3 )
+
+p4   = geom.MakePointStruct( 170.0, 100.0, 0.0 )
+seg1 = geom.MakeVector( p3, p4 )
+
+p5   = geom.MakePointStruct( 200.0, 70.0, 0.0 )
+p6   = geom.MakePointStruct( 170.0, 40.0, 0.0 )
+arc2 = geom.MakeArc( p4, p5, p6 )
+
+p7   = geom.MakePointStruct( 120.0, 30.0, 0.0 )
+arc3 = geom.MakeArc( p6, p7, p1 )
+
+
+# ---- define a closed wire with arcs and segment
+
+List1 = []
+List1.append( arc1 )
+List1.append( seg1 )
+List1.append( arc2 )
+List1.append( arc3 )
+
+ListIOR1 = []
+for S in List1 :
+    ListIOR1.append( S._get_Name() )
+wire1 = geom.MakeWire( ListIOR1 )
+
+Id_wire1 = geompy.addToStudy( wire1, "wire1")
+
+
+# ---- define a planar face with wire
+WantPlanarFace = 1 #True
+face1 = geom.MakeFace( wire1, WantPlanarFace )
+Id_face1 = geompy.addToStudy( face1, "face1")
+
+
+# ---- create a shape by extrusion
+pO = geom.MakePointStruct( 0.0, 0.0,   0.0 )
+pz = geom.MakePointStruct( 0.0, 0.0, 100.0 )
+
+prism1    = geom.MakePrism( face1, pO, pz )
+Id_prism1 = geompy.addToStudy( prism1, "prism1")
+
+# ---- create two cylinders
+
+pc1 = geom.MakePointStruct(  90.0, 50.0, -40.0 )
+pc2 = geom.MakePointStruct( 170.0, 70.0, -40.0 )
+vz  = geom.MakeDirection( pz )
+radius = 20.0
+height = 180.0
+cyl1  = geom.MakeCylinder( pc1, vz, radius, height )
+cyl2  = geom.MakeCylinder( pc2, vz, radius, height )
+
+Id_Cyl1 = geompy.addToStudy( cyl1, "cyl1" )
+Id_Cyl2 = geompy.addToStudy( cyl2, "cyl2" )
+
+# ---- cut with cyl1 
+shape  = geom.MakeBoolean( prism1, cyl1, 2 )
+
+# ---- fuse with cyl2 
+shape1 =  geom.MakeBoolean( shape, cyl2, 3 )
+
+Id_shape1 = geompy.addToStudy( shape1, "shape1")
+
+# ---- add a face sub shape in study to be meshed differently
+
+IdSubFaceList = []
+IdSubFaceList.append(10)
+sub_face = geompy.SubShapeSorted( shape1, ShapeTypeFace, IdSubFaceList )
+name     = geompy.SubShapeName( sub_face._get_Name(), shape1._get_Name() )
+
+Id_SubFace = geompy.addToStudyInFather( shape1, sub_face, name )
+
+# ---- add a face sub shape in study to be meshed differently
+
+IdSubFaceL = []
+IdSubFaceL.append(7)
+sub_face2 = geompy.SubShapeSorted( shape1, ShapeTypeFace, IdSubFaceL )
+name      = geompy.SubShapeName( sub_face2._get_Name(), shape1._get_Name() )
+
+Id_SubFace2 = geompy.addToStudyInFather( shape1, sub_face2, name )
+
+### ---------------------------- SMESH --------------------------------------
+
+# ---- launch SMESH, init a Mesh with shape 'shape1'
+gen  = smeshpy.smeshpy()
+mesh = gen.Init( Id_shape1 )
+
+idmesh = smeshgui.AddNewMesh( salome.orb.object_to_string(mesh) )
+smeshgui.SetName( idmesh, "Mesh_meca" )
+smeshgui.SetShape( Id_shape1, idmesh )
+
+print "-------------------------- NumberOfSegments"
+
+numberOfSegment = 10
+
+hypNumberOfSegment = gen.CreateHypothesis( "NumberOfSegments" )
+hypNbSeg = hypNumberOfSegment._narrow( SMESH.SMESH_NumberOfSegments )
+hypNbSeg.SetNumberOfSegments(numberOfSegment)
+print hypNbSeg.GetName()
+print hypNbSeg.GetId()
+print hypNbSeg.GetNumberOfSegments()
+
+idSeg = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypNbSeg) )
+smeshgui.SetName(idSeg, "NumberOfSegments")
+
+print "-------------------------- MaxElementArea"
+
+maxElementArea = 20
+
+hypMaxElementArea = gen.CreateHypothesis( "MaxElementArea" )
+hypArea = hypMaxElementArea._narrow( SMESH.SMESH_MaxElementArea )
+hypArea.SetMaxElementArea(maxElementArea)
+print hypArea.GetName()
+print hypArea.GetId()
+print hypArea.GetMaxElementArea()
+
+idArea = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypArea) )
+smeshgui.SetName(idArea, "MaxElementArea")
+
+print "-------------------------- MaxElementVolume"
+
+maxElementVolume = 20
+
+hypMaxElementVolume = gen.CreateHypothesis( "MaxElementVolume" )
+hypVolume = hypMaxElementVolume._narrow( SMESH.SMESH_MaxElementVolume )
+hypVolume.SetMaxElementVolume(maxElementVolume)
+print hypVolume.GetName()
+print hypVolume.GetId()
+print hypVolume.GetMaxElementVolume()
+
+idVolume = smeshgui.AddNewHypothesis( salome.orb.object_to_string(hypVolume) )
+smeshgui.SetName(idVolume, "MaxElementArea")
+
+print "-------------------------- Regular_1D"
+
+alg1D = gen.CreateHypothesis( "Regular_1D" )
+algo1D   = alg1D._narrow( SMESH.SMESH_Algo )
+listHyp =algo1D.GetCompatibleHypothesis()
+for hyp in listHyp:
+    print hyp
+algoReg1D = alg1D._narrow( SMESH.SMESH_Regular_1D )
+print algoReg1D.GetName()
+print algoReg1D.GetId()
+
+idReg1D = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(algoReg1D) )
+smeshgui.SetName( idReg1D, "Regular_1D" )
+
+print "-------------------------- MEFISTO_2D"
+
+alg2D = gen.CreateHypothesis( "MEFISTO_2D" )
+algo2D = alg2D._narrow( SMESH.SMESH_Algo )
+listHyp = algo2D.GetCompatibleHypothesis()
+for hyp in listHyp:
+    print hyp
+algoMef = alg2D._narrow( SMESH.SMESH_MEFISTO_2D )
+print algoMef.GetName()
+print algoMef.GetId()
+
+idMef = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(algoMef) )
+smeshgui.SetName( idMef, "MEFISTO_2D" )
+
+print "-------------------------- Tetra_3D"
+
+alg3D = gen.CreateHypothesis( "Tetra_3D" )
+algo3D = alg3D._narrow( SMESH.SMESH_Algo )
+listHyp = algo3D.GetCompatibleHypothesis()
+for hyp in listHyp:
+    print hyp
+algoNg = alg3D._narrow( SMESH.SMESH_Tetra_3D )
+print algoNg.GetName()
+print algoNg.GetId()
+
+idNg = smeshgui.AddNewAlgorithms( salome.orb.object_to_string(algoNg) )
+smeshgui.SetName( idNg, "Tetra_2D" )
+
+print "-------------------------- add hypothesis to main shape1"
+
+shape_mesh = salome.IDToObject( Id_shape1  )
+submesh    = mesh.GetElementsOnShape( shape_mesh )
+
+ret = mesh.AddHypothesis( shape_mesh, algoReg1D )   # Regular 1D/wire discretisation
+print ret
+ret = mesh.AddHypothesis( shape_mesh, algoMef )     # MEFISTO 2D
+print ret
+ret = mesh.AddHypothesis( shape_mesh, algoNg )     # Tetra 3D
+print ret
+ret = mesh.AddHypothesis( shape_mesh, hypNbSeg )   # nb segments
+print ret
+ret = mesh.AddHypothesis( shape_mesh, hypArea )    # max area
+print ret
+ret = mesh.AddHypothesis( shape_mesh, hypVolume )    # max volume
+print ret
+
+smeshgui.SetAlgorithms( idmesh, idReg1D );  # Regular 1D/wire discretisation
+smeshgui.SetAlgorithms( idmesh, idMef );    # MEFISTO 2D
+smeshgui.SetAlgorithms( idmesh, idNg );    # Tetra 3D
+smeshgui.SetHypothesis( idmesh, idSeg );    # nb segments
+smeshgui.SetHypothesis( idmesh, idArea );  # max area
+smeshgui.SetHypothesis( idmesh, idVolume );  # max volume
+
+sg.updateObjBrowser(1);
+
+print "-------------------------- compute the mesh of the boxe"
+ret=gen.Compute(mesh,Id_shape1)
+print ret
+log=mesh.GetLog(0) # no erase trace
+for linelog in log:
+    print linelog
+
+sg.updateObjBrowser(1)