src/Tools/Verima/Scripts/script1.py

   1 import sys
   2 import salome
   3
   4 salome.salome_init()
   5 theStudy = salome.myStudy
   6
   7 import salome_notebook
   8 notebook = salome_notebook.notebook
   9 sys.path.insert( 0, r'/local00/home/A96028/Salome/V7_main/appli_V7_main')
  10
  11 ###
  12 ### GEOM component
  13 ###
  14
  15 import GEOM
  16 from salome.geom import geomBuilder
  17 import math
  18 import SALOMEDS
  19 geompy = geomBuilder.New(theStudy)
  20
  21 ##
  22 ### SMESH component
  23 ###
  24 import  SMESH, SALOMEDS
  25 from salome.smesh import smeshBuilder
  26 smesh = smeshBuilder.New(theStudy)
  27 from salome.NETGENPlugin import NETGENPluginBuilder
  28
  29 # create a box
  30 O = geompy.MakeVertex(0, 0, 0)
  31 OX = geompy.MakeVectorDXDYDZ(1, 0, 0)
  32 OY = geompy.MakeVectorDXDYDZ(0, 1, 0)
  33 OZ = geompy.MakeVectorDXDYDZ(0, 0, 1)
  34 Box_1 = geompy.MakeBoxDXDYDZ(200, 200, 200)
  35 geompy.addToStudy( O, 'O' )
  36 geompy.addToStudy( OX, 'OX' )
  37 geompy.addToStudy( OY, 'OY' )
  38 geompy.addToStudy( OZ, 'OZ' )
  39 geompy.addToStudy( Box_1, 'Box_1' )
  40
  41 Mesh_1 = smesh.Mesh(Box_1)
  42 NETGEN_2D3D = Mesh_1.Tetrahedron(algo=smeshBuilder.NETGEN_1D2D3D)
  43 NETGEN_3D_Simple_Parameters_1 = NETGEN_2D3D.Parameters(smeshBuilder.SIMPLE)
  44 NETGEN_3D_Simple_Parameters_1.SetNumberOfSegments( 15 )
  45 NETGEN_3D_Simple_Parameters_1.LengthFromEdges()
  46 NETGEN_3D_Simple_Parameters_1.LengthFromFaces()
  47 isDone = Mesh_1.Compute()
  48
  49 smesh.SetName(Mesh_1.GetMesh(), 'Mesh_1')
  50 smesh.SetName(NETGEN_2D3D.GetAlgorithm(), 'NETGEN_2D3D')
  51 smesh.SetName(NETGEN_3D_Simple_Parameters_1, 'NETGEN 3D Simple Parameters_1')
  52
  53
  54 fichierMedResult='/tmp/tetra.med'
  55 #Mesh_1.ExportMED(fichierMedResult)
  56
  57 #____________________________________________________________
  58 from Stats.getStats import getStatsMaillage, getStatsGroupes, genHistogram
  59
  60 fichierStatMailles=fichierMedResult.replace('.med','.res')
  61 fichierStatTailles=fichierMedResult.replace('.med','.taille')
  62 fichierStatRatio=fichierMedResult.replace('.med','.ratio')
  63
  64 getStatsMaillage(Mesh_1,fichierStatMailles)
  65 getStatsGroupes(Mesh_1,fichierMedResult)
  66 genHistogram(Mesh_1, SMESH.FT_MaxElementLength3D, 20, False, fichierStatTailles,theStudy)
  67 genHistogram(Mesh_1, SMESH.FT_AspectRatio3D, 20, False, fichierStatRatio,theStudy)
  68
  69
  70