22 |
from esys.pycad import * |
from esys.pycad import * |
23 |
from esys.pycad.gmsh import Design |
from esys.pycad.gmsh import Design |
24 |
from esys.finley import MakeDomain |
from esys.finley import MakeDomain |
25 |
p0=Point(0.,0.,0.5) |
p0=Point(0.,0.,0.) |
26 |
p1=Point(1.,0.,0.5) |
p1=Point(1.,0.,0.) |
27 |
p2=Point(1.,1.,0.) |
p2=Point(1.,1.,0.) |
28 |
p3=Point(0.,1.,0.) |
p3=Point(0.,1.,0.) |
29 |
l01=Line(p0,p1) |
l01=Line(p0,p1) |
35 |
ps=PropertySet("The whole domain",s) |
ps=PropertySet("The whole domain",s) |
36 |
pl1=PropertySet("sides",l01,l23) |
pl1=PropertySet("sides",l01,l23) |
37 |
pl2=PropertySet("top and bottom",l12,l30) |
pl2=PropertySet("top and bottom",l12,l30) |
38 |
d=Design(element_size=0.2) |
d=Design(dim=2,element_size=0.2) |
39 |
d.addItems(pl1,pl2) |
d.addItems(pl1,pl2) |
40 |
d.addItems(ps) |
d.addItems(ps) |
41 |
d.setScriptFileName("quad.geo") |
d.setScriptFileName("quad.geo") |