--- /dev/null
+<!doctype HTML public "-//W3C//DTD HTML 4.0 Frameset//EN">\r
+\r
+<html>\r
+\r
+<head>\r
+<title>Defining hypotheses TUI</title>\r
+<meta http-equiv="content-type" content="text/html; charset=windows-1252">\r
+<meta name="generator" content="RoboHelp by eHelp Corporation www.ehelp.com"><style>\r
+<!--\r
+P { margin-top:0pt; margin-bottom:0pt; }\r
+LI.kadov-P { }\r
+-->\r
+</style><style type="text/css">\r
+<!--\r
+p.whs1 { font-weight:bold; }\r
+p.whs2 { font-family:'Lucida Console' , monospace; }\r
+p.whs3 { font-family:'Times New Roman' , serif; }\r
+p.whs4 { margin-top:0pt; margin-bottom:0pt; font-family:'Lucida Console' , monospace; }\r
+p.whs5 { margin-top:0pt; margin-bottom:0pt; }\r
+p.whs6 { margin-top:0pt; margin-bottom:0pt; font-family:Tahoma, sans-serif; }\r
+p.whs7 { margin-top:0pt; margin-bottom:0pt; font-family:'Times New Roman' , serif; }\r
+-->\r
+</style><script type="text/javascript" language="JavaScript">\r
+<!--\r
+if ((navigator.appName == "Netscape") && (parseInt(navigator.appVersion) == 4))\r
+{\r
+ var strNSS = "<style type='text/css'>";\r
+ strNSS += "p.whs4 {margin-top:1pt;margin-bottom:1pt; }";\r
+ strNSS += "p.whs5 {margin-top:1pt;margin-bottom:1pt; }";\r
+ strNSS += "p.whs6 {margin-top:1pt;margin-bottom:1pt; }";\r
+ strNSS += "p.whs7 {margin-top:1pt;margin-bottom:1pt; }";\r
+ strNSS +="</style>";\r
+ document.write(strNSS);\r
+}\r
+//-->\r
+</script>\r
+<script type="text/javascript" language="JavaScript" title="WebHelpInlineScript">\r
+<!--\r
+function reDo() {\r
+ if (innerWidth != origWidth || innerHeight != origHeight)\r
+ location.reload();\r
+}\r
+if ((parseInt(navigator.appVersion) == 4) && (navigator.appName == "Netscape")) {\r
+ origWidth = innerWidth;\r
+ origHeight = innerHeight;\r
+ onresize = reDo;\r
+}\r
+onerror = null; \r
+//-->\r
+</script>\r
+<style type="text/css">\r
+<!--\r
+div.WebHelpPopupMenu { position:absolute; left:0px; top:0px; z-index:4; visibility:hidden; }\r
+-->\r
+</style><script type="text/javascript" language="javascript1.2" src="whmsg.js"></script>\r
+<script type="text/javascript" language="javascript" src="whver.js"></script>\r
+<script type="text/javascript" language="javascript1.2" src="whproxy.js"></script>\r
+<script type="text/javascript" language="javascript1.2" src="whutils.js"></script>\r
+<script type="text/javascript" language="javascript1.2" src="whtopic.js"></script>\r
+</head>\r
+<body><script type="text/javascript" language="javascript1.2">\r
+<!--\r
+if (window.gbWhTopic)\r
+{\r
+ if (window.addTocInfo)\r
+ {\r
+ addTocInfo("MESH module\nTUI Scripts\nDefining Hypotheses");\r
+addButton("show",BTN_TEXT,"Show","","","","",0,0,"","","");\r
+\r
+ }\r
+ if (window.writeBtnStyle)\r
+ writeBtnStyle();\r
+\r
+ if (window.writeIntopicBar)\r
+ writeIntopicBar(1);\r
+\r
+ if (window.setRelStartPage)\r
+ {\r
+ setRelStartPage("smesh.htm");\r
+\r
+ autoSync(1);\r
+ sendSyncInfo();\r
+ sendAveInfoOut();\r
+ }\r
+}\r
+else\r
+ if (window.gbIE4)\r
+ document.location.reload();\r
+//-->\r
+</script>\r
+<h1>Defining Hypotheses and Algorithms</h1>\r
+\r
+<h3>Defining 1D Hypotheses</h3>\r
+\r
+<p class="whs1"><a name=bookmark>1D Arithmetic</a></p>\r
+\r
+<p> </p>\r
+\r
+<p class="whs2">from geompy import \r
+ *</p>\r
+\r
+<p class="whs2">import smesh</p>\r
+\r
+<p> </p>\r
+\r
+<p># create vertices</p>\r
+\r
+<p class="whs2">Point111 = MakeVertex( \r
+ 0, 0, 0)</p>\r
+\r
+<p class="whs2">Point211 = MakeVertex(10, \r
+ 0, 0)</p>\r
+\r
+<p class="whs2">Point121 = MakeVertex( \r
+ 0, 10, 0)</p>\r
+\r
+<p class="whs2">Point221 = MakeVertex(10, \r
+ 10, 0)</p>\r
+\r
+<p class="whs2">Point112 = MakeVertex( \r
+ 0, 0, 10)</p>\r
+\r
+<p class="whs2">Point212 = MakeVertex(10, \r
+ 0, 10)</p>\r
+\r
+<p class="whs2">Point122 = MakeVertex( \r
+ 0, 10, 10)</p>\r
+\r
+<p class="whs2">Point222 = MakeVertex(10, \r
+ 10, 10)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create edges</p>\r
+\r
+<p class="whs2">EdgeX111 = MakeEdge(Point111, \r
+ Point211)</p>\r
+\r
+<p class="whs2">EdgeX121 = MakeEdge(Point121, \r
+ Point221)</p>\r
+\r
+<p class="whs2">EdgeX112 = MakeEdge(Point112, \r
+ Point212)</p>\r
+\r
+<p class="whs2">EdgeX122 = MakeEdge(Point122, \r
+ Point222)</p>\r
+\r
+<p class="whs2">EdgeY11 = MakeEdge(Point111, \r
+ Point121)</p>\r
+\r
+<p class="whs2">EdgeY21 = MakeEdge(Point211, \r
+ Point221)</p>\r
+\r
+<p class="whs2">EdgeY12 = MakeEdge(Point112, \r
+ Point122)</p>\r
+\r
+<p class="whs2">EdgeY22 = MakeEdge(Point212, \r
+ Point222)</p>\r
+\r
+<p class="whs2">EdgeZ111 = MakeEdge(Point111, \r
+ Point112)</p>\r
+\r
+<p class="whs2">EdgeZ211 = MakeEdge(Point211, \r
+ Point212)</p>\r
+\r
+<p class="whs2">EdgeZ121 = MakeEdge(Point121, \r
+ Point122)</p>\r
+\r
+<p class="whs2">EdgeZ221 = MakeEdge(Point221, \r
+ Point222)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create faces</p>\r
+\r
+<p class="whs2">FaceX11 = MakeQuad(EdgeY11, \r
+ EdgeZ111, EdgeY12, EdgeZ121)</p>\r
+\r
+<p class="whs2">FaceX21 = MakeQuad(EdgeY21, \r
+ EdgeZ211, EdgeY22, EdgeZ221)</p>\r
+\r
+<p class="whs2">FaceY111 = MakeQuad(EdgeX111, \r
+ EdgeZ111, EdgeX112, EdgeZ211)</p>\r
+\r
+<p class="whs2">FaceY121 = MakeQuad(EdgeX121, \r
+ EdgeZ121, EdgeX122, EdgeZ221)</p>\r
+\r
+<p class="whs2">FaceZ11 = MakeQuad(EdgeX111, \r
+ EdgeY11, EdgeX121, EdgeY21)</p>\r
+\r
+<p class="whs2">FaceZ12 = MakeQuad(EdgeX112, \r
+ EdgeY12, EdgeX122, EdgeY22)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a solid</p>\r
+\r
+<p class="whs2">Block = MakeHexa(FaceX11, \r
+ FaceX21, FaceY111, FaceY121, FaceZ11, FaceZ12)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a compound</p>\r
+\r
+<p class="whs2">box = MakeCompound([Block])</p>\r
+\r
+<p> </p>\r
+\r
+<p># add in the study</p>\r
+\r
+<p class="whs2">box_id = addToStudy(box, \r
+ "Box compound")</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a hexahedral mesh on the box</p>\r
+\r
+<p class="whs2">hexa = smesh.Mesh(box, \r
+ "Box compound : hexahedrical mesh")</p>\r
+\r
+<p class="whs2">algo = hexa.Segment()</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "NumberOfSegments" hypothesis to cut an edge in a \r
+ fixed number of segments</p>\r
+\r
+<p class="whs2">algo.NumberOfSegments(4)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a quadrangle 2D algorithm for faces</p>\r
+\r
+<p class="whs2">hexa.Quadrangle()</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a hexahedron 3D algorithm for solids</p>\r
+\r
+<p class="whs2">hexa.Hexahedron()</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a local hypothesis</p>\r
+\r
+<p class="whs2">algo = hexa.Segment(EdgeX111)</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "Arithmetic1D" hypothesis to cut an edge in several \r
+ segments with arithmetic length increasing</p>\r
+\r
+<p class="whs2">algo.Arithmetic1D(1, \r
+ 4)</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "Propagation" hypothesis that propagates all other \r
+ hypothesis on all edges on the opposite side in case of quadrangular faces</p>\r
+\r
+<p class="whs2">algo.Propagation()</p>\r
+\r
+<p> </p>\r
+\r
+<p># compute the mesh</p>\r
+\r
+<p class="whs2">hexa.Compute() </p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<h4><a name=bookmark9>Deflection 1D and Number of Segments</a></h4>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs2">import smesh</p>\r
+\r
+<p class="whs2">import geompy</p>\r
+\r
+<p class="whs2">import salome</p>\r
+\r
+<p class="whs2">gg = salome.ImportComponentGUI("GEOM")</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create vertices</p>\r
+\r
+<p class="whs2">px = \r
+ geompy.MakeVertex(100., 0. , \r
+ 0. )</p>\r
+\r
+<p class="whs2">py = \r
+ geompy.MakeVertex(0. , \r
+ 100., 0. )</p>\r
+\r
+<p class="whs2">pz = \r
+ geompy.MakeVertex(0. , \r
+ 0. , 100.)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a vector from \r
+ two points</p>\r
+\r
+<p class="whs2">vxy = geompy.MakeVector(px, \r
+ py)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create an arc from \r
+ three points</p>\r
+\r
+<p class="whs2">arc = geompy.MakeArc(py, \r
+ pz, px)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a wire</p>\r
+\r
+<p class="whs2">wire = geompy.MakeWire([vxy, \r
+ arc])</p>\r
+\r
+<p class="whs2">isPlanarFace = 1</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a face from \r
+ the wire</p>\r
+\r
+<p class="whs2">face1 = geompy.MakeFace(wire, \r
+ isPlanarFace)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># add objects in the \r
+ study</p>\r
+\r
+<p class="whs2">id_face1 = geompy.addToStudy(face1,"Face1")</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># display faces</p>\r
+\r
+<p class="whs2">gg.createAndDisplayGO(id_face1)</p>\r
+\r
+<p class="whs2">gg.setDisplayMode(id_face1,1)</p>\r
+\r
+<p class="whs2">gg.setTransparency(id_face1,0.2)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create hexahedral mesh</p>\r
+\r
+<p class="whs2">hexa = smesh.Mesh(face1, \r
+ "Face compound : hexahedrical mesh")</p>\r
+\r
+<p class="whs2">algo = hexa.Triangle()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># define "MaxElementArea" \r
+ hypothesis to be applied to \r
+ each triangle</p>\r
+\r
+<p class="whs2">algo.MaxElementArea(30)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a quadrangle \r
+ 2D algorithm for faces</p>\r
+\r
+<p class="whs2">hexa.Quadrangle()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a local hypothesis</p>\r
+\r
+<p class="whs2">algo = hexa.Segment(wire)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># define "NumberOfSegments" \r
+ hypothesis to cut an edge in a fixed number of segments</p>\r
+\r
+<p class="whs2">algo.NumberOfSegments(6)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># define "Deflection1D" \r
+ hypothesis</p>\r
+\r
+<p class="whs2">algo.Deflection1D(1)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># compute the mesh</p>\r
+\r
+<p class="whs2">hexa.Compute() </p>\r
+\r
+<p> </p>\r
+\r
+<h4><a name=bookmark2>Start and End Length</a></h4>\r
+\r
+<p class="whs2">from geompy import \r
+ *</p>\r
+\r
+<p class="whs2">import smesh</p>\r
+\r
+<p> </p>\r
+\r
+<p># create vertices</p>\r
+\r
+<p class="whs2">Point111 = MakeVertex( \r
+ 0, 0, 0)</p>\r
+\r
+<p class="whs2">Point211 = MakeVertex(10, \r
+ 0, 0)</p>\r
+\r
+<p class="whs2">Point121 = MakeVertex( \r
+ 0, 10, 0)</p>\r
+\r
+<p class="whs2">Point221 = MakeVertex(10, \r
+ 10, 0)</p>\r
+\r
+<p class="whs2">Point112 = MakeVertex( \r
+ 0, 0, 10)</p>\r
+\r
+<p class="whs2">Point212 = MakeVertex(10, \r
+ 0, 10)</p>\r
+\r
+<p class="whs2">Point122 = MakeVertex( \r
+ 0, 10, 10)</p>\r
+\r
+<p class="whs2">Point222 = MakeVertex(10, \r
+ 10, 10)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create edges</p>\r
+\r
+<p class="whs2">EdgeX111 = MakeEdge(Point111, \r
+ Point211)</p>\r
+\r
+<p class="whs2">EdgeX121 = MakeEdge(Point121, \r
+ Point221)</p>\r
+\r
+<p class="whs2">EdgeX112 = MakeEdge(Point112, \r
+ Point212)</p>\r
+\r
+<p class="whs2">EdgeX122 = MakeEdge(Point122, \r
+ Point222)</p>\r
+\r
+<p class="whs2">EdgeY11 = MakeEdge(Point111, \r
+ Point121)</p>\r
+\r
+<p class="whs2">EdgeY21 = MakeEdge(Point211, \r
+ Point221)</p>\r
+\r
+<p class="whs2">EdgeY12 = MakeEdge(Point112, \r
+ Point122)</p>\r
+\r
+<p class="whs2">EdgeY22 = MakeEdge(Point212, \r
+ Point222)</p>\r
+\r
+<p class="whs2">EdgeZ111 = MakeEdge(Point111, \r
+ Point112)</p>\r
+\r
+<p class="whs2">EdgeZ211 = MakeEdge(Point211, \r
+ Point212)</p>\r
+\r
+<p class="whs2">EdgeZ121 = MakeEdge(Point121, \r
+ Point122)</p>\r
+\r
+<p class="whs2">EdgeZ221 = MakeEdge(Point221, \r
+ Point222)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create faces</p>\r
+\r
+<p class="whs2">FaceX11 = MakeQuad(EdgeY11, \r
+ EdgeZ111, EdgeY12, EdgeZ121)</p>\r
+\r
+<p class="whs2">FaceX21 = MakeQuad(EdgeY21, \r
+ EdgeZ211, EdgeY22, EdgeZ221)</p>\r
+\r
+<p class="whs2">FaceY111 = MakeQuad(EdgeX111, \r
+ EdgeZ111, EdgeX112, EdgeZ211)</p>\r
+\r
+<p class="whs2">FaceY121 = MakeQuad(EdgeX121, \r
+ EdgeZ121, EdgeX122, EdgeZ221)</p>\r
+\r
+<p class="whs2">FaceZ11 = MakeQuad(EdgeX111, \r
+ EdgeY11, EdgeX121, EdgeY21)</p>\r
+\r
+<p class="whs2">FaceZ12 = MakeQuad(EdgeX112, \r
+ EdgeY12, EdgeX122, EdgeY22)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a solid</p>\r
+\r
+<p class="whs2">Block = MakeHexa(FaceX11, \r
+ FaceX21, FaceY111, FaceY121, FaceZ11, FaceZ12)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a compound</p>\r
+\r
+<p class="whs2">box = MakeCompound([Block])</p>\r
+\r
+<p> </p>\r
+\r
+<p># add in the study</p>\r
+\r
+<p class="whs2">box_id = addToStudy(box, \r
+ "Box compound")</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a hexahedral mesh on the box</p>\r
+\r
+<p class="whs2">hexa = smesh.Mesh(box, \r
+ "Box compound : hexahedrical mesh")</p>\r
+\r
+<p class="whs2">algo = hexa.Segment()</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "NumberOfSegments" hypothesis to cut an edge in a \r
+ fixed number of segments</p>\r
+\r
+<p class="whs2">algo.NumberOfSegments(4)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a quadrangle 2D algorithm for faces</p>\r
+\r
+<p class="whs2">hexa.Quadrangle()</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a hexahedron 3D algorithm for solids</p>\r
+\r
+<p class="whs2">hexa.Hexahedron()</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a local hypothesis</p>\r
+\r
+<p class="whs2">algo = hexa.Segment(EdgeX111)</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "StartEndLength" hypothesis to cut an edge in several \r
+ segments with increasing geometric length </p>\r
+\r
+<p class="whs2">algo.StartEndLength(1, \r
+ 6)</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "Propagation" hypothesis that propagates all other \r
+ hypothesis on all edges on the opposite side in case of quadrangular faces</p>\r
+\r
+<p class="whs2">algo.Propagation()</p>\r
+\r
+<p> </p>\r
+\r
+<p># compute the mesh</p>\r
+\r
+<p class="whs2">hexa.Compute() </p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<h4><a name=bookmark3>Average Length</a></h4>\r
+\r
+<p class="whs2">from geompy import \r
+ *</p>\r
+\r
+<p class="whs2">import smesh</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create vertices</p>\r
+\r
+<p class="whs2">Point111 = MakeVertex( \r
+ 0, 0, 0)</p>\r
+\r
+<p class="whs2">Point211 = MakeVertex(10, \r
+ 0, 0)</p>\r
+\r
+<p class="whs2">Point121 = MakeVertex( \r
+ 0, 10, 0)</p>\r
+\r
+<p class="whs2">Point221 = MakeVertex(10, \r
+ 10, 0)</p>\r
+\r
+<p class="whs2">Point112 = MakeVertex( \r
+ 0, 0, 10)</p>\r
+\r
+<p class="whs2">Point212 = MakeVertex(10, \r
+ 0, 10)</p>\r
+\r
+<p class="whs2">Point122 = MakeVertex( \r
+ 0, 10, 10)</p>\r
+\r
+<p class="whs2">Point222 = MakeVertex(10, \r
+ 10, 10)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create edges</p>\r
+\r
+<p class="whs2">EdgeX111 = MakeEdge(Point111, \r
+ Point211)</p>\r
+\r
+<p class="whs2">EdgeX121 = MakeEdge(Point121, \r
+ Point221)</p>\r
+\r
+<p class="whs2">EdgeX112 = MakeEdge(Point112, \r
+ Point212)</p>\r
+\r
+<p class="whs2">EdgeX122 = MakeEdge(Point122, \r
+ Point222)</p>\r
+\r
+<p class="whs2">EdgeY11 = MakeEdge(Point111, \r
+ Point121)</p>\r
+\r
+<p class="whs2">EdgeY21 = MakeEdge(Point211, \r
+ Point221)</p>\r
+\r
+<p class="whs2">EdgeY12 = MakeEdge(Point112, \r
+ Point122)</p>\r
+\r
+<p class="whs2">EdgeY22 = MakeEdge(Point212, \r
+ Point222)</p>\r
+\r
+<p class="whs2">EdgeZ111 = MakeEdge(Point111, \r
+ Point112)</p>\r
+\r
+<p class="whs2">EdgeZ211 = MakeEdge(Point211, \r
+ Point212)</p>\r
+\r
+<p class="whs2">EdgeZ121 = MakeEdge(Point121, \r
+ Point122)</p>\r
+\r
+<p class="whs2">EdgeZ221 = MakeEdge(Point221, \r
+ Point222)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create faces</p>\r
+\r
+<p class="whs2">FaceX11 = MakeQuad(EdgeY11, \r
+ EdgeZ111, EdgeY12, EdgeZ121)</p>\r
+\r
+<p class="whs2">FaceX21 = MakeQuad(EdgeY21, \r
+ EdgeZ211, EdgeY22, EdgeZ221)</p>\r
+\r
+<p class="whs2">FaceY111 = MakeQuad(EdgeX111, \r
+ EdgeZ111, EdgeX112, EdgeZ211)</p>\r
+\r
+<p class="whs2">FaceY121 = MakeQuad(EdgeX121, \r
+ EdgeZ121, EdgeX122, EdgeZ221)</p>\r
+\r
+<p class="whs2">FaceZ11 = MakeQuad(EdgeX111, \r
+ EdgeY11, EdgeX121, EdgeY21)</p>\r
+\r
+<p class="whs2">FaceZ12 = MakeQuad(EdgeX112, \r
+ EdgeY12, EdgeX122, EdgeY22)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a solid</p>\r
+\r
+<p class="whs2">Block = MakeHexa(FaceX11, \r
+ FaceX21, FaceY111, FaceY121, FaceZ11, FaceZ12)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a compound</p>\r
+\r
+<p class="whs2">box = MakeCompound([Block])</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># add in the study</p>\r
+\r
+<p class="whs2">box_id = addToStudy(box, \r
+ "Box compound")</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a hexahedral \r
+ mesh on the box</p>\r
+\r
+<p class="whs2">hexa = smesh.Mesh(box, \r
+ "Box compound : hexahedrical mesh")</p>\r
+\r
+<p class="whs2">algo = hexa.Segment()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># define "NumberOfSegments" \r
+ hypothesis to cut an edge in a fixed number of segments</p>\r
+\r
+<p class="whs2">algo.NumberOfSegments(4)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a quadrangle \r
+ 2D algorithm for faces</p>\r
+\r
+<p class="whs2">hexa.Quadrangle()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a hexahedron \r
+ 3D algorithm for solids</p>\r
+\r
+<p class="whs2">hexa.Hexahedron()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a local hypothesis</p>\r
+\r
+<p class="whs2">algo = hexa.Segment(EdgeX111)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># define "LocalLength" \r
+ hypothesis to cut an edge in several segments with the same length</p>\r
+\r
+<p class="whs2">algo.LocalLength(2)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># define "Propagation" \r
+ hypothesis that propagates all other hypothesis on all edges on<span style="font-family: 'Times New Roman', serif;"> \r
+ the opposite side in case of quadrangular faces</span></p>\r
+\r
+<p class="whs2">algo.Propagation()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># compute the mesh</p>\r
+\r
+<p class="whs2">hexa.Compute() </p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<h3>Defining 2D and 3D hypotheses</h3>\r
+\r
+<h4><a name=bookmark4>Maximum Element Area</a></h4>\r
+\r
+<p class="whs4">import smesh</p>\r
+\r
+<p class="whs4">import geompy</p>\r
+\r
+<p class="whs4">import salome</p>\r
+\r
+<p class="whs4">gg = salome.ImportComponentGUI("GEOM")</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create vertices</p>\r
+\r
+<p class="whs4">px = \r
+ geompy.MakeVertex(100., 0. , \r
+ 0. )</p>\r
+\r
+<p class="whs4">py = \r
+ geompy.MakeVertex(0. , \r
+ 100., 0. )</p>\r
+\r
+<p class="whs4">pz = \r
+ geompy.MakeVertex(0. , \r
+ 0. , 100.)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create a vector from \r
+ two points</p>\r
+\r
+<p class="whs4">vxy = geompy.MakeVector(px, \r
+ py)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create an arc from three \r
+ points</p>\r
+\r
+<p class="whs4">arc = geompy.MakeArc(py, \r
+ pz, px)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create a wire</p>\r
+\r
+<p class="whs4">wire = geompy.MakeWire([vxy, \r
+ arc])</p>\r
+\r
+<p class="whs4">isPlanarFace = 1</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create a face from the \r
+ wire</p>\r
+\r
+<p class="whs4">face1 = geompy.MakeFace(wire, \r
+ isPlanarFace)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># add objects in the study</p>\r
+\r
+<p class="whs4">id_face1 = geompy.addToStudy(face1,"Face1")</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># display faces</p>\r
+\r
+<p class="whs4">gg.createAndDisplayGO(id_face1)</p>\r
+\r
+<p class="whs4">gg.setDisplayMode(id_face1,1)</p>\r
+\r
+<p class="whs4">gg.setTransparency(id_face1,0.2)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create a hexahedral mesh</p>\r
+\r
+<p class="whs4">hexa = smesh.Mesh(face1, \r
+ "Face compound : hexahedrical mesh")</p>\r
+\r
+<p class="whs4">algo = hexa.Triangle()</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># define "MaxElementArea" \r
+ hypothesis to be applied to each triangle</p>\r
+\r
+<p class="whs4">algo.MaxElementArea(7)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create a quadrangle 2D \r
+ algorithm for faces</p>\r
+\r
+<p class="whs4">hexa.Quadrangle()</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create a local hypothesis</p>\r
+\r
+<p class="whs6">algo = hexa.Segment(wire)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># define "NumberOfSegments" \r
+ hypothesis to cut an edge in a fixed number of segments</p>\r
+\r
+<p class="whs4">algo.NumberOfSegments(10)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># compute the mesh</p>\r
+\r
+<p class="whs4">hexa.Compute() </p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<h4><a name=bookmark5>Maximum Element Volume</a></h4>\r
+\r
+<p class="whs4"><span style="font-family: 'Lucida Console', monospace;">import \r
+ salome</span></p>\r
+\r
+<p class="whs4">import geompy</p>\r
+\r
+<p class="whs4">import StdMeshers</p>\r
+\r
+<p class="whs4">smesh = salome.lcc.FindOrLoadComponent("FactoryServer", \r
+ "SMESH")</p>\r
+\r
+<p class="whs4">smeshgui = salome.ImportComponentGUI("SMESH")</p>\r
+\r
+<p class="whs4">smeshgui.Init(salome.myStudyId);</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># create a box</p>\r
+\r
+<p class="whs4">box = geompy.MakeCylinderRH(30, \r
+ 50) #MakeBox(0., 0., 0., 100., 200., 300.)</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># add the box to the \r
+ study</p>\r
+\r
+<p class="whs4">idbox = geompy.addToStudy(box, \r
+ "box")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># create vertices</p>\r
+\r
+<p class="whs4">px = \r
+ geompy.MakeVertex(100., 0. , \r
+ 0. )</p>\r
+\r
+<p class="whs4">py = \r
+ geompy.MakeVertex(0. , \r
+ 100., 0. )</p>\r
+\r
+<p class="whs4">pz = \r
+ geompy.MakeVertex(0. , \r
+ 0. , 100.)</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># create a vector from \r
+ two points</p>\r
+\r
+<p class="whs4">vxy = geompy.MakeVector(px, \r
+ py)</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># create an arc from \r
+ three points</p>\r
+\r
+<p class="whs4">arc = geompy.MakeArc(py, \r
+ pz, px)</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># create a wire</p>\r
+\r
+<p class="whs4">wire = geompy.MakeWire([vxy, \r
+ arc])</p>\r
+\r
+<p class="whs4">isPlanarFace = 1</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># create a face from \r
+ the wire</p>\r
+\r
+<p class="whs4">#face1</p>\r
+\r
+<p class="whs4">box = geompy.MakeFace(wire, \r
+ isPlanarFace)</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># add objects in the \r
+ study</p>\r
+\r
+<p class="whs4">id_face1 = geompy.addToStudy(box,"Face1")</p>\r
+\r
+<p class="whs4">#geompy.addToStudy(face1,"Face1")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># display faces</p>\r
+\r
+<p class="whs4">gg.createAndDisplayGO(id_face1)</p>\r
+\r
+<p class="whs4">gg.setDisplayMode(id_face1,1)</p>\r
+\r
+<p class="whs4">gg.setTransparency(id_face1,0.2)</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># create a hypothesis</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ create Hypothesis"</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ NumberOfSegments"</p>\r
+\r
+<p class="whs4">numberOfSegments \r
+ = 7</p>\r
+\r
+<p class="whs4">hypNbSeg = smesh.CreateHypothesis("NumberOfSegments", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">hypNbSeg.SetNumberOfSegments(numberOfSegments)</p>\r
+\r
+<p class="whs4">print hypNbSeg.GetName()</p>\r
+\r
+<p class="whs4">print hypNbSeg.GetId()</p>\r
+\r
+<p class="whs4">print hypNbSeg.GetNumberOfSegments()</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(hypNbSeg), \r
+ "NumberOfSegments_10")</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ MaxElementArea"</p>\r
+\r
+<p class="whs4">maxElementArea = \r
+ 800</p>\r
+\r
+<p class="whs4">hypArea = smesh.CreateHypothesis("MaxElementArea", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">hypArea.SetMaxElementArea(maxElementArea)</p>\r
+\r
+<p class="whs4">print hypArea.GetName()</p>\r
+\r
+<p class="whs4">print hypArea.GetId()</p>\r
+\r
+<p class="whs4">print hypArea.GetMaxElementArea()</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(hypArea), \r
+ "MaxElementArea_800")</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ MaxElementVolume"</p>\r
+\r
+<p class="whs4">maxElementVolume \r
+ = 900</p>\r
+\r
+<p class="whs4">hypVolume = smesh.CreateHypothesis("MaxElementVolume", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">hypVolume.SetMaxElementVolume(maxElementVolume)</p>\r
+\r
+<p class="whs4">print hypVolume.GetName()</p>\r
+\r
+<p class="whs4">print hypVolume.GetId()</p>\r
+\r
+<p class="whs4">print hypVolume.GetMaxElementVolume()</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(hypVolume), \r
+ "MaxElementVolume_900")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># create algorithms</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ create Algorithms"</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ Regular_1D"</p>\r
+\r
+<p class="whs4">regular1D = smesh.CreateHypothesis("Regular_1D", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(regular1D), \r
+ "Wire Discretisation")</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ MEFISTO_2D"</p>\r
+\r
+<p class="whs4">mefisto2D = smesh.CreateHypothesis("MEFISTO_2D", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(mefisto2D), \r
+ "MEFISTO_2D")</p>\r
+\r
+<p class="whs4">#print "-------------------------- \r
+ Hexa_3D (Hexahedron meshing algorithm)"</p>\r
+\r
+<p class="whs4">hexa3D = smesh.CreateHypothesis("Hexa_3D", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(hexa3D), \r
+ "HEXA_3D")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># initialize a mesh \r
+ with the box</p>\r
+\r
+<p class="whs4">mesh = smesh.CreateMesh(box)</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(mesh), \r
+ "MeshBox")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># add a hypothesis \r
+ to the box</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ add hypothesis to the box"</p>\r
+\r
+<p class="whs4">mesh.AddHypothesis(box,hypNbSeg)</p>\r
+\r
+<p class="whs4">mesh.AddHypothesis(box,hypArea)</p>\r
+\r
+<p class="whs4">mesh.AddHypothesis(box,hypVolume)</p>\r
+\r
+<p class="whs4">mesh.AddHypothesis(box,regular1D)</p>\r
+\r
+<p class="whs4">mesh.AddHypothesis(box,mefisto2D)</p>\r
+\r
+<p class="whs4">mesh.AddHypothesis(box,hexa3D)</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs7"># compute the mesh</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ compute the mesh of the box"</p>\r
+\r
+<p class="whs4">ret = smesh.Compute(mesh,box)</p>\r
+\r
+<p class="whs4">print ret</p>\r
+\r
+<p class="whs4">if ret == 0:</p>\r
+\r
+<p class="whs4"> print \r
+ "probleme when computing the mesh"</p>\r
+\r
+<p class="whs4">salome.sg.updateObjBrowser(1)</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<h3>Defining Additional Hypotheses</h3>\r
+\r
+<h4><a name=bookmark6>Length from Edges</a></h4>\r
+\r
+<p class="whs2">import smesh</p>\r
+\r
+<p class="whs2">import geompy</p>\r
+\r
+<p class="whs2">import salome</p>\r
+\r
+<p class="whs2">gg = salome.ImportComponentGUI("GEOM")</p>\r
+\r
+<p> </p>\r
+\r
+<p># create sketchers</p>\r
+\r
+<p class="whs2">sketcher1 = geompy.MakeSketcher("Sketcher:F \r
+ 0 0:TT 70 0:TT 70 70:TT 0 70:WW")</p>\r
+\r
+<p class="whs2">sketcher2 = geompy.MakeSketcher("Sketcher:F \r
+ 20 20:TT 50 20:TT 50 50:TT 20 50:WW")</p>\r
+\r
+<p class="whs2">isPlanarFace = 1</p>\r
+\r
+<p> </p>\r
+\r
+<p># create faces from two wires</p>\r
+\r
+<p class="whs2">face1 = geompy.MakeFaces([sketcher1, \r
+ sketcher2],isPlanarFace)</p>\r
+\r
+<p> </p>\r
+\r
+<p># add objects in the study</p>\r
+\r
+<p class="whs2">id_face1 = geompy.addToStudy(face1,"Face1")</p>\r
+\r
+<p> </p>\r
+\r
+<p># display faces</p>\r
+\r
+<p class="whs2">gg.createAndDisplayGO(id_face1)</p>\r
+\r
+<p class="whs2">gg.setDisplayMode(id_face1,1)</p>\r
+\r
+<p class="whs2">gg.setTransparency(id_face1,0.2)</p>\r
+\r
+<p> </p>\r
+\r
+<p># create a mesh</p>\r
+\r
+<p class="whs2">hexa = smesh.Mesh(face1, \r
+ "Face : triangle 2D mesh")</p>\r
+\r
+<p class="whs2">algo = hexa.Triangle()</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "MaxElementArea" hypothesis to be applied to each \r
+ triangle</p>\r
+\r
+<p class="whs2">algo.MaxElementArea(30)</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "LengthFromEdges" hypothesis to build triangles based \r
+ on the length of the edges taken from the wire</p>\r
+\r
+<p class="whs2">algo.LengthFromEdges()</p>\r
+\r
+<p> </p>\r
+\r
+<p># create local hypothesis</p>\r
+\r
+<p class="whs2">algo = hexa.Segment()</p>\r
+\r
+<p> </p>\r
+\r
+<p># define "NumberOfSegments" hypothesis to cut an edge in a \r
+ fixed number of segments</p>\r
+\r
+<p class="whs2">algo.NumberOfSegments(2)</p>\r
+\r
+<p> </p>\r
+\r
+<p># compute the mesh</p>\r
+\r
+<p class="whs2">hexa.Compute() </p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<h4><a name=bookmark7>Propagation</a></h4>\r
+\r
+<p class="whs2">from geompy import \r
+ *</p>\r
+\r
+<p class="whs2">import smesh</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create vertices</p>\r
+\r
+<p class="whs2">Point111 = MakeVertex( \r
+ 0, 0, 0)</p>\r
+\r
+<p class="whs2">Point211 = MakeVertex(10, \r
+ 0, 0)</p>\r
+\r
+<p class="whs2">Point121 = MakeVertex( \r
+ 0, 10, 0)</p>\r
+\r
+<p class="whs2">Point221 = MakeVertex(10, \r
+ 10, 0)</p>\r
+\r
+<p class="whs2">Point112 = MakeVertex( \r
+ 0, 0, 10)</p>\r
+\r
+<p class="whs2">Point212 = MakeVertex(10, \r
+ 0, 10)</p>\r
+\r
+<p class="whs2">Point122 = MakeVertex( \r
+ 0, 10, 10)</p>\r
+\r
+<p class="whs2">Point222 = MakeVertex(10, \r
+ 10, 10)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create edges</p>\r
+\r
+<p class="whs2">EdgeX111 = MakeEdge(Point111, \r
+ Point211)</p>\r
+\r
+<p class="whs2">EdgeX121 = MakeEdge(Point121, \r
+ Point221)</p>\r
+\r
+<p class="whs2">EdgeX112 = MakeEdge(Point112, \r
+ Point212)</p>\r
+\r
+<p class="whs2">EdgeX122 = MakeEdge(Point122, \r
+ Point222)</p>\r
+\r
+<p class="whs2">EdgeY11 = MakeEdge(Point111, \r
+ Point121)</p>\r
+\r
+<p class="whs2">EdgeY21 = MakeEdge(Point211, \r
+ Point221)</p>\r
+\r
+<p class="whs2">EdgeY12 = MakeEdge(Point112, \r
+ Point122)</p>\r
+\r
+<p class="whs2">EdgeY22 = MakeEdge(Point212, \r
+ Point222)</p>\r
+\r
+<p class="whs2">EdgeZ111 = MakeEdge(Point111, \r
+ Point112)</p>\r
+\r
+<p class="whs2">EdgeZ211 = MakeEdge(Point211, \r
+ Point212)</p>\r
+\r
+<p class="whs2">EdgeZ121 = MakeEdge(Point121, \r
+ Point122)</p>\r
+\r
+<p class="whs2">EdgeZ221 = MakeEdge(Point221, \r
+ Point222)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create faces</p>\r
+\r
+<p class="whs2">FaceX11 = MakeQuad(EdgeY11, \r
+ EdgeZ111, EdgeY12, EdgeZ121)</p>\r
+\r
+<p class="whs2">FaceX21 = MakeQuad(EdgeY21, \r
+ EdgeZ211, EdgeY22, EdgeZ221)</p>\r
+\r
+<p class="whs2">FaceY111 = MakeQuad(EdgeX111, \r
+ EdgeZ111, EdgeX112, EdgeZ211)</p>\r
+\r
+<p class="whs2">FaceY121 = MakeQuad(EdgeX121, \r
+ EdgeZ121, EdgeX122, EdgeZ221)</p>\r
+\r
+<p class="whs2">FaceZ11 = MakeQuad(EdgeX111, \r
+ EdgeY11, EdgeX121, EdgeY21)</p>\r
+\r
+<p class="whs2">FaceZ12 = MakeQuad(EdgeX112, \r
+ EdgeY12, EdgeX122, EdgeY22)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a solid</p>\r
+\r
+<p class="whs2">Block = MakeHexa(FaceX11, \r
+ FaceX21, FaceY111, FaceY121, FaceZ11, FaceZ12)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a compound</p>\r
+\r
+<p class="whs2">box = MakeCompound([Block])</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># add in the study</p>\r
+\r
+<p class="whs2">box_id = addToStudy(box, \r
+ "Box compound")</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a hexahedral \r
+ mesh on the box</p>\r
+\r
+<p class="whs2">hexa = smesh.Mesh(box, \r
+ "Box compound : hexahedrical mesh")</p>\r
+\r
+<p class="whs2">algo = hexa.Segment()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># define "NumberOfSegments" \r
+ hypothesis to cut an edge in a fixed number of segments</p>\r
+\r
+<p class="whs2">algo.NumberOfSegments(4)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a quadrangle \r
+ 2D algorithm for faces</p>\r
+\r
+<p class="whs2">hexa.Quadrangle()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a hexahedron \r
+ 3D algorithm for solids</p>\r
+\r
+<p class="whs2">hexa.Hexahedron()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># create a local hypothesis</p>\r
+\r
+<p class="whs2">algo = hexa.Segment(EdgeX111)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs2"><span style="font-family: 'Times New Roman', serif;"># \r
+ define "Arithmetic1D" hypothesis to cut an edge in several segments \r
+ with arithmetic length</span> </p>\r
+\r
+<p class="whs2">increasing</p>\r
+\r
+<p class="whs2">algo.Arithmetic1D(1, \r
+ 4)</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># define "Propagation" \r
+ hypothesis that propagatea all other hypothesis on all edges on the opposite \r
+ side in case of quadrangular faces</p>\r
+\r
+<p class="whs2">algo.Propagation()</p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<p class="whs3"># compute the mesh</p>\r
+\r
+<p class="whs2">hexa.Compute() </p>\r
+\r
+<p class="whs2"> </p>\r
+\r
+<h3><a name=bookmark8>Defining Meshing Algorithms</a></h3>\r
+\r
+<p class="whs4"><span style="font-family: 'Lucida Console', monospace;">import \r
+ salome</span></p>\r
+\r
+<p class="whs4">import StdMeshers</p>\r
+\r
+<p class="whs4">import NETGENPlugin</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs4">smesh = salome.lcc.FindOrLoadComponent("FactoryServer", \r
+ "SMESH")</p>\r
+\r
+<p class="whs4">smeshgui = salome.ImportComponentGUI("SMESH")</p>\r
+\r
+<p class="whs4">smeshgui.Init(salome.myStudyId);</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<p class="whs5"># create algorithms</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ create Algorithms"</p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ Regular_1D (Wire discretisation)"</p>\r
+\r
+<p class="whs4">regular1D = smesh.CreateHypothesis("Regular_1D", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(regular1D), \r
+ "Wire Discretisation")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ MEFISTO_2D (Triangle meshing algorithm)"</p>\r
+\r
+<p class="whs4">mefisto2D = smesh.CreateHypothesis("MEFISTO_2D", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(mefisto2D), \r
+ "MEFISTO_2D")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ Quadrangle_2D (Quadrangle meshing algorithm)"</p>\r
+\r
+<p class="whs4">quad2D = smesh.CreateHypothesis( \r
+ "Quadrangle_2D", "libStdMeshersEngine.so" )</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(quad2D), \r
+ "Quadrangle_2D")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ Hexa_3D (Hexahedron meshing algorithm)"</p>\r
+\r
+<p class="whs4">hexa3D = smesh.CreateHypothesis("Hexa_3D", \r
+ "libStdMeshersEngine.so")</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(hexa3D), \r
+ "HEXA_3D")</p>\r
+\r
+<p class="whs4"> </p>\r
+\r
+<p class="whs4">print "-------------------------- \r
+ NETGEN_3D (Tetrahedron meshing algorithm)"</p>\r
+\r
+<p class="whs4">netgen3D = smesh.CreateHypothesis("NETGEN_3D", \r
+ "libNETGENEngine.so")</p>\r
+\r
+<p class="whs4">smeshgui.SetName(salome.ObjectToID(netgen3D), \r
+ "NETGEN_3D")</p>\r
+\r
+<p class="whs4">salome.sg.updateObjBrowser(1)</p>\r
+\r
+<p class="whs5"> </p>\r
+\r
+<script type="text/javascript" language="javascript1.2">\r
+<!--\r
+if (window.writeIntopicBar)\r
+ writeIntopicBar(0);\r
+//-->\r
+</script>\r
+</body>\r
+</html>\r