test/python/run_linearPDEsOnFinley3.py


########################################################
#
# Copyright (c) 2003-2018 by The University of Queensland
# Earth Systems Science Computational Center (ESSCC)
# 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
#
########################################################

from __future__ import print_function, division

__copyright__="""Copyright (c) 2003-2018 by The University of Queensland
Earth Systems Science Computational Center (ESSCC)
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"

"""
Test suite for the linearPDE  and pdetools test on finley

:remark:

:var __author__: name of author
:var __licence__: licence agreement
:var __url__: url entry point on documentation
:var __version__: version
:var __date__: date of the version
"""

__author__="Lutz Gross, l.gross@uq.edu.au"

import os

import esys.escriptcore.utestselect as unittest
from esys.escriptcore.testing import *
from test_linearPDEs import Test_Helmholtz, Test_LameEquation, Test_Poisson
from test_assemblage import Test_assemblage_2Do1_Contact, Test_assemblage_2Do2_Contact, Test_assemblage_3Do1_Contact, Test_assemblage_3Do2_Contact
from esys.finley import Rectangle, ReadMesh

try:
     FINLEY_TEST_DATA=os.environ['FINLEY_TEST_DATA']
except KeyError:
     FINLEY_TEST_DATA='.'

FINLEY_TEST_MESH_PATH=os.path.join(FINLEY_TEST_DATA,"data_meshes")

NE=8 # number of element in each spatial direction (must be even)

class Test_LameOnFinley(Test_LameEquation):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
        self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
   def tearDown(self):
        del self.domain

class Test_PoissonOnFinley(Test_Poisson):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
        self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
   def tearDown(self):
        del self.domain

class Test_HelmholtzOnFinley(Test_Helmholtz):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
        self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
   def tearDown(self):
        del self.domain

class Test_AssemblePDEwithFinley_2Do1_Contact(Test_assemblage_2Do1_Contact):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
       self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do1_Contact.fly"))
   def tearDown(self):
        del self.domain

class Test_AssemblePDEwithFinley_2Do2_Contact(Test_assemblage_2Do2_Contact):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
       self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do2_Contact.fly"))
   def tearDown(self):
        del self.domain

class Test_AssemblePDEwithFinley_3Do1_Contact(Test_assemblage_3Do1_Contact):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
       self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do1_Contact.fly"))
   def tearDown(self):
        del self.domain

class Test_AssemblePDEwithFinley_3Do2_Contact(Test_assemblage_3Do2_Contact):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
       self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do2_Contact.fly"))
   def tearDown(self):
        del self.domain


class Test_AssemblePDEwithFinley_2Do1_Contact_withElementsOnFace(Test_assemblage_2Do1_Contact):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
       self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do1_Contact_withElementsOnFace.fly"))
   def tearDown(self):
        del self.domain

class Test_AssemblePDEwithFinley_2Do2_Contact_withElementsOnFace(Test_assemblage_2Do2_Contact):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
       self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do2_Contact_withElementsOnFace.fly"))
   def tearDown(self):
        del self.domain

class Test_AssemblePDEwithFinley_3Do1_Contact_withElementsOnFace(Test_assemblage_3Do1_Contact):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
       self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do1_Contact_withElementsOnFace.fly"))
   def tearDown(self):
        del self.domain

class Test_AssemblePDEwithFinley_3Do2_Contact_withElementsOnFace(Test_assemblage_3Do2_Contact):
   RES_TOL=1.e-7
   ABS_TOL=1.e-8
   def setUp(self):
       self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do2_Contact_withElementsOnFace.fly"))
   def tearDown(self):
        del self.domain

if __name__ == '__main__':
    run_tests(__name__, exit_on_failure=True)

1
2	########################################################
3	#
4	# Copyright (c) 2003-2018 by The University of Queensland
5	# Earth Systems Science Computational Center (ESSCC)
6	# http://www.uq.edu.au
7	#
8	# Primary Business: Queensland, Australia
9	# Licensed under the Apache License, version 2.0
10	# http://www.apache.org/licenses/LICENSE-2.0
11	#
12	########################################################
13
14	from __future__ import print_function, division
15
16	__copyright__="""Copyright (c) 2003-2018 by The University of Queensland
17	Earth Systems Science Computational Center (ESSCC)
18	http://www.uq.edu.au
19	Primary Business: Queensland, Australia"""
20	__license__="""Licensed under the Apache License, version 2.0
21	http://www.apache.org/licenses/LICENSE-2.0"""
22	__url__="https://launchpad.net/escript-finley"
23
24	"""
25	Test suite for the linearPDE and pdetools test on finley
26
27	:remark:
28
29	:var __author__: name of author
30	:var __licence__: licence agreement
31	:var __url__: url entry point on documentation
32	:var __version__: version
33	:var __date__: date of the version
34	"""
35
36	__author__="Lutz Gross, l.gross@uq.edu.au"
37
38	import os
39
40	import esys.escriptcore.utestselect as unittest
41	from esys.escriptcore.testing import *
42	from test_linearPDEs import Test_Helmholtz, Test_LameEquation, Test_Poisson
43	from test_assemblage import Test_assemblage_2Do1_Contact, Test_assemblage_2Do2_Contact, Test_assemblage_3Do1_Contact, Test_assemblage_3Do2_Contact
44	from esys.finley import Rectangle, ReadMesh
45
46	try:
47	FINLEY_TEST_DATA=os.environ['FINLEY_TEST_DATA']
48	except KeyError:
49	FINLEY_TEST_DATA='.'
50
51	FINLEY_TEST_MESH_PATH=os.path.join(FINLEY_TEST_DATA,"data_meshes")
52
53	NE=8 # number of element in each spatial direction (must be even)
54
55	class Test_LameOnFinley(Test_LameEquation):
56	RES_TOL=1.e-7
57	ABS_TOL=1.e-8
58	def setUp(self):
59	self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
60	def tearDown(self):
61	del self.domain
62
63	class Test_PoissonOnFinley(Test_Poisson):
64	RES_TOL=1.e-7
65	ABS_TOL=1.e-8
66	def setUp(self):
67	self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
68	def tearDown(self):
69	del self.domain
70
71	class Test_HelmholtzOnFinley(Test_Helmholtz):
72	RES_TOL=1.e-7
73	ABS_TOL=1.e-8
74	def setUp(self):
75	self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
76	def tearDown(self):
77	del self.domain
78
79	class Test_AssemblePDEwithFinley_2Do1_Contact(Test_assemblage_2Do1_Contact):
80	RES_TOL=1.e-7
81	ABS_TOL=1.e-8
82	def setUp(self):
83	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do1_Contact.fly"))
84	def tearDown(self):
85	del self.domain
86
87	class Test_AssemblePDEwithFinley_2Do2_Contact(Test_assemblage_2Do2_Contact):
88	RES_TOL=1.e-7
89	ABS_TOL=1.e-8
90	def setUp(self):
91	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do2_Contact.fly"))
92	def tearDown(self):
93	del self.domain
94
95	class Test_AssemblePDEwithFinley_3Do1_Contact(Test_assemblage_3Do1_Contact):
96	RES_TOL=1.e-7
97	ABS_TOL=1.e-8
98	def setUp(self):
99	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do1_Contact.fly"))
100	def tearDown(self):
101	del self.domain
102
103	class Test_AssemblePDEwithFinley_3Do2_Contact(Test_assemblage_3Do2_Contact):
104	RES_TOL=1.e-7
105	ABS_TOL=1.e-8
106	def setUp(self):
107	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do2_Contact.fly"))
108	def tearDown(self):
109	del self.domain
110
111
112	class Test_AssemblePDEwithFinley_2Do1_Contact_withElementsOnFace(Test_assemblage_2Do1_Contact):
113	RES_TOL=1.e-7
114	ABS_TOL=1.e-8
115	def setUp(self):
116	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do1_Contact_withElementsOnFace.fly"))
117	def tearDown(self):
118	del self.domain
119
120	class Test_AssemblePDEwithFinley_2Do2_Contact_withElementsOnFace(Test_assemblage_2Do2_Contact):
121	RES_TOL=1.e-7
122	ABS_TOL=1.e-8
123	def setUp(self):
124	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do2_Contact_withElementsOnFace.fly"))
125	def tearDown(self):
126	del self.domain
127
128	class Test_AssemblePDEwithFinley_3Do1_Contact_withElementsOnFace(Test_assemblage_3Do1_Contact):
129	RES_TOL=1.e-7
130	ABS_TOL=1.e-8
131	def setUp(self):
132	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do1_Contact_withElementsOnFace.fly"))
133	def tearDown(self):
134	del self.domain
135
136	class Test_AssemblePDEwithFinley_3Do2_Contact_withElementsOnFace(Test_assemblage_3Do2_Contact):
137	RES_TOL=1.e-7
138	ABS_TOL=1.e-8
139	def setUp(self):
140	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do2_Contact_withElementsOnFace.fly"))
141	def tearDown(self):
142	del self.domain
143
144	if __name__ == '__main__':
145	run_tests(__name__, exit_on_failure=True)
146