~vcs-imports/escript-finley/trunk

« back to all changes in this revision

Viewing changes to dudley/test/python/rectangle.py

  • Committer: jfenwick
  • Date: 2010-10-11 01:48:14 UTC
  • Revision ID: svn-v4:77569008-7704-0410-b7a0-a92fef0b09fd:trunk:3259
Merging dudley and scons updates from branches

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
 
 
2
########################################################
 
3
#
 
4
# Copyright (c) 2003-2010 by University of Queensland
 
5
# Earth Systems Science Computational Center (ESSCC)
 
6
# http://www.uq.edu.au/esscc
 
7
#
 
8
# Primary Business: Queensland, Australia
 
9
# Licensed under the Open Software License version 3.0
 
10
# http://www.opensource.org/licenses/osl-3.0.php
 
11
#
 
12
########################################################
 
13
 
 
14
__copyright__="""Copyright (c) 2003-2010 by University of Queensland
 
15
Earth Systems Science Computational Center (ESSCC)
 
16
http://www.uq.edu.au/esscc
 
17
Primary Business: Queensland, Australia"""
 
18
__license__="""Licensed under the Open Software License version 3.0
 
19
http://www.opensource.org/licenses/osl-3.0.php"""
 
20
__url__="https://launchpad.net/escript-finley"
 
21
 
 
22
from esys.escript import *
 
23
from esys.pycad import *
 
24
from esys.pycad.gmsh import Design
 
25
from esys.dudley import MakeDomain
 
26
 
 
27
 
 
28
p0=Point(0.,0.)
 
29
p1=Point(1.,0.)
 
30
p2=Point(1.,1.)
 
31
p3=Point(0.,1.)
 
32
 
 
33
l01=Line(p0,p1)
 
34
l12=Line(p1,p2)
 
35
l23=Line(p2,p3)
 
36
l30=Line(p3,p0)
 
37
 
 
38
s=PlaneSurface(CurveLoop(l01,l12,l23,l30))
 
39
des=Design(dim=2, order=1, element_size = 1, keep_files=True)
 
40
des.setMeshFileName("rec.geo")
 
41
des.addItems(s)
 
42
 
 
43
dom=MakeDomain(des)
 
44
dom.write("rec.fly")