test/python/rectangle.py


##############################################################################
#
# Copyright (c) 2003-2018 by The University of Queensland
# http://www.uq.edu.au
#
# Primary Business: Queensland, Australia
# Licensed under the Apache License, version 2.0
# http://www.apache.org/licenses/LICENSE-2.0
#
# Development until 2012 by Earth Systems Science Computational Center (ESSCC)
# Development 2012-2013 by School of Earth Sciences
# Development from 2014 by Centre for Geoscience Computing (GeoComp)
#
##############################################################################

from __future__ import print_function, division

__copyright__="""Copyright (c) 2003-2018 by The University of Queensland
http://www.uq.edu.au
Primary Business: Queensland, Australia"""
__license__="""Licensed under the Apache License, version 2.0
http://www.apache.org/licenses/LICENSE-2.0"""
__url__="https://launchpad.net/escript-finley"

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.setMeshFileName("rec.geo")
des.addItems(s)

dom=MakeDomain(des)
dom.write("rec.fly")
1
2	##############################################################################
3	#
4	# Copyright (c) 2003-2018 by The University of Queensland
5	# http://www.uq.edu.au
6	#
7	# Primary Business: Queensland, Australia
8	# Licensed under the Apache License, version 2.0
9	# http://www.apache.org/licenses/LICENSE-2.0
10	#
11	# Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12	# Development 2012-2013 by School of Earth Sciences
13	# Development from 2014 by Centre for Geoscience Computing (GeoComp)
14	#
15	##############################################################################
16
17	from __future__ import print_function, division
18
19	__copyright__="""Copyright (c) 2003-2018 by The University of Queensland
20	http://www.uq.edu.au
21	Primary Business: Queensland, Australia"""
22	__license__="""Licensed under the Apache License, version 2.0
23	http://www.apache.org/licenses/LICENSE-2.0"""
24	__url__="https://launchpad.net/escript-finley"
25
26	from esys.escript import *
27	from esys.pycad import *
28	from esys.pycad.gmsh import Design
29	from esys.finley import MakeDomain
30
31
32	p0=Point(0.,0.)
33	p1=Point(1.,0.)
34	p2=Point(1.,1.)
35	p3=Point(0.,1.)
36
37	l01=Line(p0,p1)
38	l12=Line(p1,p2)
39	l23=Line(p2,p3)
40	l30=Line(p3,p0)
41
42	s=PlaneSurface(CurveLoop(l01,l12,l23,l30))
43	des=Design(dim=2, order=1, element_size = 1, keep_files=True)
44	des.setMeshFileName("rec.geo")
45	des.addItems(s)
46
47	dom=MakeDomain(des)
48	dom.write("rec.fly")