src/SMESH_SWIG/ex30_tepal.py

   1 #  -*- coding: iso-8859-1 -*-
   2 # Copyright (C) 2007-2013  CEA/DEN, EDF R&D, OPEN CASCADE
   3 #
   4 # This library is free software; you can redistribute it and/or
   5 # modify it under the terms of the GNU Lesser General Public
   6 # License as published by the Free Software Foundation; either
   7 # version 2.1 of the License.
   8 #
   9 # This library is distributed in the hope that it will be useful,
  10 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 # Lesser General Public License for more details.
  13 #
  14 # You should have received a copy of the GNU Lesser General Public
  15 # License along with this library; if not, write to the Free Software
  16 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20
  21 # ====================================================================
  22 #
  23 import os
  24
  25 import geompy
  26 import smesh
  27
  28 # Parameters
  29 # ----------
  30
  31 results = "/tmp/ZZ"
  32
  33 radius =  50
  34 height = 200
  35
  36 # Build a cylinder
  37 # ----------------
  38
  39 base = geompy.MakeVertex(0, 0, 0)
  40 direction = geompy.MakeVectorDXDYDZ(0, 0, 1)
  41
  42 cylinder = geompy.MakeCylinder(base, direction, radius, height)
  43
  44 geompy.addToStudy(cylinder, "Cylinder")
  45
  46 # Define a mesh on a geometry
  47 # ---------------------------
  48
  49 smesh.SetCurrentStudy(salome.myStudy)
  50
  51 m = smesh.Mesh(cylinder)
  52
  53 # 2D mesh with BLSURF
  54 # -------------------
  55
  56 algo2d = m.Triangle(smesh.BLSURF)
  57
  58 algo2d.SetPhysicalMesh(1)
  59 algo2d.SetPhySize(5)
  60
  61 algo2d.SetGeometricMesh(0)
  62
  63 # 3D mesh with tepal
  64 # ------------------
  65
  66 algo3d = m.Tetrahedron(smesh.GHS3DPRL)
  67
  68 algo3d.SetMEDName(results)
  69 algo3d.SetNbPart(4)
  70 algo3d.SetBackground(False)
  71 algo3d.SetKeepFiles(False)
  72 algo3d.SetToMeshHoles(True)
  73
  74 # Launch meshers
  75 # --------------
  76
  77 status = m.Compute()
  78
  79 # Test if ok
  80 # ----------
  81
  82 if os.access(results+".xml", os.F_OK):
  83     print "Ok: tepal"
  84 else:
  85     print "KO: tepal"