test/python/rectangle.py

# $Id$

__copyright__="""  Copyright (c) 2006, 2007 by ACcESS MNRF
                    http://www.access.edu.au
                Primary Business: Queensland, Australia"""
__license__="""Licensed under the Open Software License version 3.0
             http://www.opensource.org/licenses/osl-3.0.php"""

from esys.escript import *
from esys.pycad import *
from esys.pycad.gmsh import Design
from esys.finley import MakeDomain


p0=Point(0.,0.)
p1=Point(1.,0.)
p2=Point(1.,1.)
p3=Point(0.,1.)

l01=Line(p0,p1)
l12=Line(p1,p2)
l23=Line(p2,p3)
l30=Line(p3,p0)

s=PlaneSurface(CurveLoop(l01,l12,l23,l30))
des=Design(dim=2, order=1, element_size = 1, keep_files=True)
des.addItems(s)

dom=MakeDomain(des)
dom.write("rec.fly")
1	# $Id$
2
3	__copyright__=""" Copyright (c) 2006, 2007 by ACcESS MNRF
4	http://www.access.edu.au
5	Primary Business: Queensland, Australia"""
6	__license__="""Licensed under the Open Software License version 3.0
7	http://www.opensource.org/licenses/osl-3.0.php"""
8
9	from esys.escript import *
10	from esys.pycad import *
11	from esys.pycad.gmsh import Design
12	from esys.finley import MakeDomain
13
14
15	p0=Point(0.,0.)
16	p1=Point(1.,0.)
17	p2=Point(1.,1.)
18	p3=Point(0.,1.)
19
20	l01=Line(p0,p1)
21	l12=Line(p1,p2)
22	l23=Line(p2,p3)
23	l30=Line(p3,p0)
24
25	s=PlaneSurface(CurveLoop(l01,l12,l23,l30))
26	des=Design(dim=2, order=1, element_size = 1, keep_files=True)
27	des.addItems(s)
28
29	dom=MakeDomain(des)
30	dom.write("rec.fly")