X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FSMESH_SWIG%2Fbatchmode_mefisto.py;h=d3e5ced2da7534d4b0f160eb41edbb0b110f20fd;hp=6fa30717e72f5d491af3f8741d887ecea981ff2b;hb=1821a9c35f90be93d3d55b1e1db9114902eaf6f8;hpb=c3bf92bd87b770fd81631a3853f7f5bb1ac6a4e8 diff --git a/src/SMESH_SWIG/batchmode_mefisto.py b/src/SMESH_SWIG/batchmode_mefisto.py index 6fa30717e..d3e5ced2d 100644 --- a/src/SMESH_SWIG/batchmode_mefisto.py +++ b/src/SMESH_SWIG/batchmode_mefisto.py @@ -1,20 +1,37 @@ +# -*- coding: iso-8859-1 -*- +# Copyright (C) 2007-2014 CEA/DEN, EDF R&D, OPEN CASCADE +# +# Copyright (C) 2003-2007 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, or (at your option) any later version. +# +# 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.salome-platform.org/ or email : webmaster.salome@opencascade.com +# + import os import re import batchmode_salome import batchmode_geompy import batchmode_smesh +from salome.StdMeshers import StdMeshersBuilder -geom = batchmode_geompy.geom smesh = batchmode_smesh.smesh - -geom.GetCurrentStudy(batchmode_salome.myStudyId) smesh.SetCurrentStudy(batchmode_salome.myStudy) -ShapeType = batchmode_smesh.ShapeType - -import StdMeshers - def CreateMesh (theFileName, area, len = None, nbseg = None): if not(os.path.isfile(theFileName)) or re.search("\.brep$", theFileName) is None : @@ -31,120 +48,80 @@ def CreateMesh (theFileName, area, len = None, nbseg = None): # ---- Import shape from BREP file and add it to the study - shape_mesh = geom.ImportBREP(theFileName) - Id_shape = batchmode_geompy.addToStudy( shape_mesh, "shape_mesh") - + shape_mesh = batchmode_geompy.Import(theFileName, "BREP") + Id_shape = batchmode_geompy.addToStudy(shape_mesh, "shape_mesh") # ---- SMESH + print "-------------------------- create mesh" + mesh = smesh.Mesh(shape_mesh) - # ---- create Hypothesis - print "-------------------------- create Hypothesis" if (len is not None): print "-------------------------- LocalLength" - hypLength1 = smesh.CreateHypothesis("LocalLength", "libStdMeshersEngine.so") - hypLength1.SetLength(len) - print "Hypothesis type: ", hypLength1.GetName() - print "Hypothesis ID: ", hypLength1.GetId() + algoReg = mesh.Segment() + hypLength1 = algoReg.LocalLength(len) + print "Hypothesis type : ", hypLength1.GetName() + print "Hypothesis ID : ", hypLength1.GetId() print "Hypothesis Value: ", hypLength1.GetLength() if (nbseg is not None): print "-------------------------- NumberOfSegments" - hypNbSeg1 = smesh.CreateHypothesis("NumberOfSegments", "libStdMeshersEngine.so") - hypNbSeg1.SetNumberOfSegments(nbseg) - print "Hypothesis type: ", hypNbSeg1.GetName() - print "Hypothesis ID: ", hypNbSeg1.GetId() + algoReg = mesh.Segment() + hypNbSeg1 = algoReg.NumberOfSegments(nbseg) + print "Hypothesis type : ", hypNbSeg1.GetName() + print "Hypothesis ID : ", hypNbSeg1.GetId() print "Hypothesis Value: ", hypNbSeg1.GetNumberOfSegments() - if (area == "LengthFromEdges"): print "-------------------------- LengthFromEdges" - hypLengthFromEdges = smesh.CreateHypothesis("LengthFromEdges", "libStdMeshersEngine.so") - hypLengthFromEdges.SetMode(1) - print "Hypothesis type: ", hypLengthFromEdges.GetName() - print "Hypothesis ID: ", hypLengthFromEdges.GetId() + algoMef = mesh.Triangle() + hypLengthFromEdges = algoMef.LengthFromEdges(1) + print "Hypothesis type : ", hypLengthFromEdges.GetName() + print "Hypothesis ID : ", hypLengthFromEdges.GetId() print "LengthFromEdges Mode: ", hypLengthFromEdges.GetMode() else: print "-------------------------- MaxElementArea" - hypArea1 = smesh.CreateHypothesis("MaxElementArea", "libStdMeshersEngine.so") - hypArea1.SetMaxElementArea(area) - print "Hypothesis type: ", hypArea1.GetName() - print "Hypothesis ID: ", hypArea1.GetId() + algoMef = mesh.Triangle() + hypArea1 = algoMef.MaxElementArea(area) + print "Hypothesis type : ", hypArea1.GetName() + print "Hypothesis ID : ", hypArea1.GetId() print "Hypothesis Value: ", hypArea1.GetMaxElementArea() - print "-------------------------- Regular_1D" - algoReg = smesh.CreateHypothesis("Regular_1D", "libStdMeshersEngine.so") - listHyp = algoReg.GetCompatibleHypothesis() for hyp in listHyp: print hyp print "Algo name: ", algoReg.GetName() - print "Algo ID: ", algoReg.GetId() + print "Algo ID : ", algoReg.GetId() print "-------------------------- MEFISTO_2D" - algoMef = smesh.CreateHypothesis("MEFISTO_2D", "libStdMeshersEngine.so") - - listHyp=algoMef.GetCompatibleHypothesis() + listHyp = algoMef.GetCompatibleHypothesis() for hyp in listHyp: print hyp print "Algo name: ", algoMef.GetName() - print "Algo ID: ", algoMef.GetId() - + print "Algo ID : ", algoMef.GetId() # ---- add hypothesis to shape - print "-------------------------- add hypothesis to shape" - mesh = smesh.CreateMesh(shape_mesh) - - ret = mesh.AddHypothesis(shape_mesh, algoReg) - print "Add Regular_1D algo .... ", - print ret - - if (nbseg is not None): - ret=mesh.AddHypothesis(shape_mesh, hypNbSeg1) - print "Add Number Of Segements algo .... ", - print ret - - if (len is not None): - ret=mesh.AddHypothesis(shape_mesh,hypLength1) - print "Add Local Length algo .... ", - print ret - - ret=mesh.AddHypothesis(shape_mesh, algoMef) - print "Add MEFISTO_2D algo....", - print ret - - if (area == "LengthFromEdges"): - ret = mesh.AddHypothesis( shape_mesh, hypLengthFromEdges) # length from edge - print "Add Length From Edges algo .... ", - print ret - else: - ret=mesh.AddHypothesis(shape_mesh, hypArea1) - print "Add Max Triangle Area algo .... ", - print ret - - # ---- compute mesh - print "-------------------------- compute mesh" - ret=smesh.Compute(mesh,shape_mesh) + ret = mesh.Compute() print "Compute Mesh .... ", print ret - log=mesh.GetLog(0); # no erase trace + log = mesh.GetLog(0); # no erase trace #for linelog in log: # print linelog print "------------ INFORMATION ABOUT MESH ------------" - print "Number of nodes: ", mesh.NbNodes() - print "Number of edges: ", mesh.NbEdges() - print "Number of faces: ", mesh.NbFaces() + print "Number of nodes : ", mesh.NbNodes() + print "Number of edges : ", mesh.NbEdges() + print "Number of faces : ", mesh.NbFaces() print "Number of triangles: ", mesh.NbTriangles() return mesh