/[escript]/trunk/ripley/test/python/run_linearPDEsOnRipley1.py
ViewVC logotype

Contents of /trunk/ripley/test/python/run_linearPDEsOnRipley1.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3792 - (show annotations)
Wed Feb 1 06:16:25 2012 UTC (7 years, 1 month ago) by caltinay
File MIME type: text/x-python
File size: 3189 byte(s)
Merged ripley rectangular domain into trunk.

1
2 ########################################################
3 #
4 # Copyright (c) 2003-2012 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-2012 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 """
23 Test suite for the linearPDE and pdetools on ripley
24
25 :remark:
26
27 :var __author__: name of author
28 :var __licence__: licence agreement
29 :var __url__: url entry point on documentation
30 :var __version__: version
31 :var __date__: date of the version
32 """
33
34 __author__="Lutz Gross, l.gross@uq.edu.au"
35
36 import os
37
38 import unittest
39 from test_linearPDEs import Test_Poisson,Test_LinearPDE, Test_LinearPDE_noLumping, Test_TransportPDE
40 from test_assemblage import Test_assemblage_2Do1, Test_assemblage_2Do2, Test_assemblage_3Do1, Test_assemblage_3Do2, \
41 Test_assemblage_2Do1_Contact,Test_assemblage_2Do2_Contact, Test_assemblage_3Do1_Contact, Test_assemblage_3Do2_Contact
42 from test_pdetools import Test_pdetools, Test_pdetools_noLumping
43 from esys.escript import *
44 from esys.ripley import Rectangle,Brick, ReadMesh
45 import sys
46
47
48 try:
49 RIPLEY_TEST_DATA=os.environ['RIPLEY_TEST_DATA']
50 except KeyError:
51 RIPLEY_TEST_DATA='.'
52
53 NE=10 # number of element in each spatial direction (must be even)
54 mpiSize=getMPISizeWorld()
55
56 class Test_LinearPDEOnRipleyRect(Test_LinearPDE,Test_pdetools,Test_assemblage_2Do1, Test_TransportPDE):
57 RES_TOL=1.e-7
58 ABS_TOL=1.e-8
59 def setUp(self):
60 for x in [int(sqrt(mpiSize)),2,3,5,7,1]:
61 NX=x
62 NY=mpiSize/x
63 if NX*NY == mpiSize:
64 break
65 self.domain=Rectangle(n0=NE*NX-1, n1=NE*NY-1, l0=1., l1=1., d0=NX, d1=NY)
66 self.order = 1
67 def tearDown(self):
68 del self.domain
69
70 class Test_LinearPDEOnRipleyBrick(Test_LinearPDE,Test_pdetools,Test_assemblage_3Do1, Test_TransportPDE):
71 RES_TOL=1.e-7
72 ABS_TOL=1.e-8
73 def setUp(self):
74 for x in [(int(mpiSize**(1/3.)),int(mpiSize**(1/3.))),(2,3),(2,2),(1,2),(1,1)]:
75 NX=x[0]
76 NY=x[1]
77 NZ=mpiSize/(x[0]*x[1])
78 if NX*NY*NZ == mpiSize:
79 break
80
81 self.domain = Brick(n0=NE*NX-1, n1=NE*NY-1, n2=NE*NZ-1, l0=1., l1=1., l2=1., d0=NX, d1=NY, d2=NZ)
82 self.order = 1
83 def tearDown(self):
84 del self.domain
85
86
87 if __name__ == '__main__':
88 suite = unittest.TestSuite()
89 suite.addTest(unittest.makeSuite(Test_LinearPDEOnRipleyRect))
90 suite.addTest(unittest.makeSuite(Test_LinearPDEOnRipleyBrick))
91 s=unittest.TextTestRunner(verbosity=2).run(suite)
92 if not s.wasSuccessful(): sys.exit(1)
93

  ViewVC Help
Powered by ViewVC 1.1.26