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	phornby	1894
2			########################################################
3			#
4	jfenwick	6651	# Copyright (c) 2003-2018 by The University of Queensland
5	phornby	1894	# Earth Systems Science Computational Center (ESSCC)
6	jfenwick	3981	# http://www.uq.edu.au
7	phornby	1894	#
8			# Primary Business: Queensland, Australia
9	jfenwick	6112	# Licensed under the Apache License, version 2.0
10			# http://www.apache.org/licenses/LICENSE-2.0
11	phornby	1894	#
12			########################################################
13
14	sshaw	5706	from __future__ import print_function, division
15
16	jfenwick	6651	__copyright__="""Copyright (c) 2003-2018 by The University of Queensland
17	phornby	1894	Earth Systems Science Computational Center (ESSCC)
18	jfenwick	3981	http://www.uq.edu.au
19	phornby	1894	Primary Business: Queensland, Australia"""
20	jfenwick	6112	__license__="""Licensed under the Apache License, version 2.0
21			http://www.apache.org/licenses/LICENSE-2.0"""
22	jfenwick	2344	__url__="https://launchpad.net/escript-finley"
23	phornby	1894
24			"""
25			Test suite for the linearPDE and pdetools test on finley
26
27	jfenwick	2625	:remark:
28	phornby	1894
29	jfenwick	2625	: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	phornby	1894	"""
35
36			__author__="Lutz Gross, l.gross@uq.edu.au"
37
38			import os
39
40	jfenwick	4938	import esys.escriptcore.utestselect as unittest
41	sshaw	4984	from esys.escriptcore.testing import *
42	caltinay	6282	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	phornby	1894
46			try:
47			FINLEY_TEST_DATA=os.environ['FINLEY_TEST_DATA']
48			except KeyError:
49			FINLEY_TEST_DATA='.'
50
51	gross	2315	FINLEY_TEST_MESH_PATH=os.path.join(FINLEY_TEST_DATA,"data_meshes")
52	phornby	1894
53	gross	2315	NE=8 # number of element in each spatial direction (must be even)
54	phornby	1894
55	gross	2325	class Test_LameOnFinley(Test_LameEquation):
56			RES_TOL=1.e-7
57			ABS_TOL=1.e-8
58			def setUp(self):
59	caltinay	6282	self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
60	gross	2325	def tearDown(self):
61			del self.domain
62	caltinay	6282
63	phornby	1894	class Test_PoissonOnFinley(Test_Poisson):
64			RES_TOL=1.e-7
65			ABS_TOL=1.e-8
66			def setUp(self):
67	caltinay	6282	self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
68	phornby	1894	def tearDown(self):
69			del self.domain
70
71	gross	2323	class Test_HelmholtzOnFinley(Test_Helmholtz):
72			RES_TOL=1.e-7
73			ABS_TOL=1.e-8
74			def setUp(self):
75	caltinay	6282	self.domain = Rectangle(NE, NE, 2, useElementsOnFace=0, useFullElementOrder=True)
76	gross	2323	def tearDown(self):
77			del self.domain
78	phornby	1894
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	caltinay	6282	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do1_Contact.fly"))
84	phornby	1894	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	caltinay	6282	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do2_Contact.fly"))
92	phornby	1894	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	caltinay	6282	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do1_Contact.fly"))
100	phornby	1894	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	caltinay	6282	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do2_Contact.fly"))
108	phornby	1894	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	caltinay	6282	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do1_Contact_withElementsOnFace.fly"))
117	phornby	1894	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	caltinay	6282	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_2Do2_Contact_withElementsOnFace.fly"))
125	phornby	1894	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	caltinay	6282	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do1_Contact_withElementsOnFace.fly"))
133	phornby	1894	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	caltinay	6282	self.domain = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"mesh_3Do2_Contact_withElementsOnFace.fly"))
141	phornby	1894	def tearDown(self):
142			del self.domain
143
144			if __name__ == '__main__':
145	sshaw	4984	run_tests(__name__, exit_on_failure=True)
146	phornby	1894