/[escript]/trunk/pycad/examples/quad.py
ViewVC logotype

Annotation of /trunk/pycad/examples/quad.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1809 - (hide annotations)
Thu Sep 25 06:43:44 2008 UTC (12 years, 11 months ago) by ksteube
File MIME type: text/x-python
File size: 1713 byte(s)
Copyright updated in all python files

1 ksteube 1809
2     ########################################################
3     #
4     # Copyright (c) 2003-2008 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-2008 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__="http://www.uq.edu.au/esscc/escript-finley"
21    
22 gross 1010 """
23     a simple 1x1 quad`
24    
25     @var __author__: name of author
26     @var __licence__: licence agreement
27     @var __url__: url entry point on documentation
28     @var __version__: version
29     @var __date__: date of the version
30     """
31    
32     __author__="Lutz Gross, l.gross@uq.edu.au"
33    
34     from esys.pycad import *
35     from esys.pycad.gmsh import Design
36     from esys.finley import MakeDomain
37 gross 1044 from esys.escript import getTagNames
38 gross 1034 p0=Point(0.,0.,0.)
39     p1=Point(1.,0.,0.)
40 gross 1010 p2=Point(1.,1.,0.)
41     p3=Point(0.,1.,0.)
42     l01=Line(p0,p1)
43     l12=Line(p1,p2)
44     l23=Line(p2,p3)
45     l30=Line(p3,p0)
46     c=CurveLoop(l01,l12,l23,l30)
47     s=PlaneSurface(c)
48 gross 1044 ps=PropertySet("The_whole_domain",s)
49 gross 1010 pl1=PropertySet("sides",l01,l23)
50 gross 1044 pl2=PropertySet("top_and_bottom",l12,l30)
51 gross 1716
52 gross 1338 d=Design(dim=2,element_size=0.005)
53 gross 1010 d.addItems(pl1,pl2)
54     d.addItems(ps)
55     d.setScriptFileName("quad.geo")
56     d.setMeshFileName("quad.msh")
57     dom=MakeDomain(d,integrationOrder=-1, reducedIntegrationOrder=-1, optimizeLabeling=True)
58 gross 1044 print getTagNames(dom)
59 gross 1010 dom.write("quad.fly")

  ViewVC Help
Powered by ViewVC 1.1.26