Mesh on 5 Edges ?

eXav wrote:

ADAM Erwan wrote:

eXav wrote:
Sorry...
Figure:
Point1 (0,0);
Point2 (2,0);
Point3 (2,1);
Point4 (1,2);
Point5 (0,2);
Point5 (0,0);
Like a house

Something like that ?

http://erwan.adam.free.fr/salome/scripts/Mesh_on_5_edges_00.png

Yes... And how it made???...

Save the following :

# -----------

# --
# Geom part ...
# --

import geompy

p1 = geompy.MakeVertex(0.0, 0.0, 0.0)
p2 = geompy.MakeVertex(2.0, 0.0, 0.0)
p3 = geompy.MakeVertex(2.0, 1.0, 0.0)
p4 = geompy.MakeVertex(1.0, 2.0, 0.0)
p5 = geompy.MakeVertex(0.0, 2.0, 0.0)

points = [p1, p2, p3, p4, p5]
edges = [geompy.MakeEdge(points[i], points[(i+1)%len(points)]) for i in range(len(points))]
wire = geompy.MakeWire(edges)
face = geompy.MakeFace(wire, 1)

shape = face

p12 = geompy.MakeVertex(1.0, 0.0, 0.0)

tools = [
    geompy.MakeEdge(p12, p4),
    ]

shape = geompy.MakePartition([shape], tools)

tool_p1 = geompy.MakeVertex(0.0, 1.0, 0.0)
tool_p2 = geompy.MakeVertex(1.0, 1.0, 0.0)
tool_p3 = geompy.MakeVertex(2.0, 0.5, 0.0)

tools = [
    geompy.MakeEdge(tool_p1, tool_p2),
    geompy.MakeEdge(tool_p2, tool_p3),
    ]

shape = geompy.MakePartition([shape], tools)
geompy.addToStudy(shape, "shape"

salome.sg.updateObjBrowser(0)

# --
# Mesh part ...
# CAUTION : I use here the corba API of SMESH
# which is not the "official" API ... I've got my own reasons !
# --

import salome
import SMESH
import StdMeshers

smesh = salome.lcc.FindOrLoadComponent("FactoryServer", "SMESH"
smesh.SetCurrentStudy(salome.myStudy)
mesh = smesh.CreateMesh(shape)

mesh.AddHypothesis(shape, smesh.CreateHypothesis("Regular_1D", "StdMeshersEngine")

hypo = smesh.CreateHypothesis("NumberOfSegments", "StdMeshersEngine"
hypo.SetNumberOfSegments(10)
mesh.AddHypothesis(shape, hypo)

mesh.AddHypothesis(shape, smesh.CreateHypothesis("Quadrangle_2D", "StdMeshersEngine")

smesh.Compute(mesh, shape)

salome.sg.updateObjBrowser(0)

# -----------

in foo.py and launch :

runSalome --execute=foo.py

it should do the job

E.A.