/[escript]/trunk/dudley/test/python/fixme_run_generators.py
ViewVC logotype

Diff of /trunk/dudley/test/python/fixme_run_generators.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1388 by trankine, Fri Jan 11 07:45:58 2008 UTC revision 2722 by gross, Fri Oct 16 06:45:01 2009 UTC
# Line 1  Line 1 
1    
2    ########################################################
3  #  #
4  # $Id$  # Copyright (c) 2003-2009 by University of Queensland
5  #  # Earth Systems Science Computational Center (ESSCC)
6  #######################################################  # http://www.uq.edu.au/esscc
 #  
 #           Copyright 2003-2007 by ACceSS MNRF  
 #       Copyright 2007 by University of Queensland  
 #  
 #                http://esscc.uq.edu.au  
 #        Primary Business: Queensland, Australia  
 #  Licensed under the Open Software License version 3.0  
 #     http://www.opensource.org/licenses/osl-3.0.php  
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-2009 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  checks the mesh generators against the reference meshes in test_meshes and test the finley integration schemes.  checks the mesh generators against the reference meshes in test_meshes and test the finley integration schemes.
24  """  """
25    
 __copyright__="""  Copyright (c) 2006 by ACcESS MNRF  
                     http://www.access.edu.au  
                 Primary Business: Queensland, Australia"""  
 __license__="""Licensed under the Open Software License version 3.0  
              http://www.opensource.org/licenses/osl-3.0.php"""  
   
26  import sys  import sys
27  import os  import os
28  import unittest  import unittest
# Line 40  try: Line 39  try:
39  except KeyError:  except KeyError:
40       FINLEY_WORKDIR='.'       FINLEY_WORKDIR='.'
41    
42  FINLEY_TEST_MESH_PATH=FINLEY_TEST_DATA+"/data_meshes/"  FINLEY_TEST_MESH_PATH=os.path.join(FINLEY_TEST_DATA,"data_meshes")
43  if os.name == "nt":  #if os.name == "nt":
44     FINLEY_TEST_MESH_PATH = FINLEY_TEST_MESH_PATH+"win32/"  #   FINLEY_TEST_MESH_PATH = os.path.join(FINLEY_TEST_MESH_PATH,"win32")
45  FINLEY_WORKDIR_PATH=FINLEY_WORKDIR+"/"  FINLEY_WORKDIR_PATH=FINLEY_WORKDIR
46    
47  TEST_FILE_EXT=".test"  TEST_FILE_PRE="test_"
48  class Test_Generators(unittest.TestCase):  class Test_Generators(unittest.TestCase):
49    
50     def checker(self,dom,reference):     def checker(self,dom,reference):
51        dom_file=FINLEY_WORKDIR_PATH+TEST_FILE_EXT        dom_file=os.path.join(FINLEY_WORKDIR_PATH,TEST_FILE_PRE+reference)
52        dom.write(dom_file)        dom.write(dom_file)
53  # Uncomment this section to dump the files for regression testing  # Uncomment this section to dump the files for regression testing
54  #      if True:  #      if True:
55  #         dom.write(FINLEY_TEST_MESH_PATH+reference)  #         dom.write(os.path.join(FINLEY_TEST_MESH_PATH,reference))
56        dom_string=open(dom_file).read().splitlines()        dom_string=open(dom_file).read().splitlines()
57        ref_string=open(FINLEY_TEST_MESH_PATH+reference).read().splitlines()        ref_string=open(os.path.join(FINLEY_TEST_MESH_PATH,reference)).read().splitlines()
58        self.failUnlessEqual(len(dom_string),len(ref_string),"number of lines in mesh files does not match reference")        self.failUnlessEqual(len(dom_string),len(ref_string),"number of lines in mesh files does not match reference")
59        for l in range(1,len(ref_string)):        for l in range(1,len(ref_string)):
60           self.failUnlessEqual(dom_string[l].strip(),ref_string[l].strip(),"line %d (%s) in mesh files does not match reference (%s)"%(l,ref_string[l].strip(),dom_string[l].strip()))       line=dom_string[l].strip()
61         if os.name == "nt":
62               line=line.replace("e+00","e+0").replace("e-00","e-0")
63             self.failUnlessEqual(line,ref_string[l].strip(),"line %d (%s) in mesh files does not match reference (%s)"%(l,ref_string[l].strip(),line))
64    
65     def test_hex_2D_order1(self):     def test_hex_2D_order1(self):
66        file="hex_2D_order1.msh"        file="hex_2D_order1.msh"
# Line 75  class Test_Generators(unittest.TestCase) Line 77  class Test_Generators(unittest.TestCase)
77        my_dom=Rectangle(1,1,2)        my_dom=Rectangle(1,1,2)
78        self.checker(my_dom,file)        self.checker(my_dom,file)
79    
80       def test_hex_2D_order1_macro(self):
81          file="hex_2D_order1_macro.msh"
82          my_dom=Rectangle(1,1,-1)
83          self.checker(my_dom,file)
84    
85     def test_hex_2D_order2_onFace(self):     def test_hex_2D_order2_onFace(self):
86        file="hex_2D_order2_onFace.msh"        file="hex_2D_order2_onFace.msh"
87        my_dom=Rectangle(1,1,2,useElementsOnFace=1)        my_dom=Rectangle(1,1,2,useElementsOnFace=1)
# Line 95  class Test_Generators(unittest.TestCase) Line 102  class Test_Generators(unittest.TestCase)
102        my_dom=Brick(1,1,1,2)        my_dom=Brick(1,1,1,2)
103        self.checker(my_dom,file)        self.checker(my_dom,file)
104    
105       def test_hex_3D_order2(self):
106          file="hex_3D_order1_macro.msh"
107          my_dom=Brick(1,1,1,-1)
108          self.checker(my_dom,file)
109    
110     def test_hex_3D_order2_onFace(self):     def test_hex_3D_order2_onFace(self):
111        file="hex_3D_order2_onFace.msh"        file="hex_3D_order2_onFace.msh"
112        my_dom=Brick(1,1,1,2,useElementsOnFace=1)        my_dom=Brick(1,1,1,2,useElementsOnFace=1)
# Line 170  class Test_GMSHReader(unittest.TestCase) Line 182  class Test_GMSHReader(unittest.TestCase)
182        ref_string=open(reference_file).read().splitlines()        ref_string=open(reference_file).read().splitlines()
183        self.failUnlessEqual(len(dom_string),len(ref_string),"number of lines in mesh files does not match reference")        self.failUnlessEqual(len(dom_string),len(ref_string),"number of lines in mesh files does not match reference")
184        for l in range(1,len(ref_string)):        for l in range(1,len(ref_string)):
185           self.failUnlessEqual(dom_string[l].strip(),ref_string[l].strip(),"line %d (%s) in mesh files does not match reference (%s)"%(l,ref_string[l].strip(),dom_string[l].strip()))       line=dom_string[l].strip()
186         if os.name == "nt":
187               line=line.replace("e+00","e+0").replace("e-00","e-0")
188             self.failUnlessEqual(line,ref_string[l].strip(),"line %d (%s) in mesh files does not match reference (%s)"%(l,ref_string[l].strip(),line))
189    
190     def test_Tri3(self):     def test_Tri3(self):
191         file="tri3_gmsh.msh"           file="tri3_gmsh.msh"
192         ref ="tri3.fly"           ref ="tri3.fly"
193         test = FINLEY_WORKDIR+os.sep+"tri3_test.fly"           test = os.path.join(FINLEY_WORKDIR,"tri3_test.fly")
194         dom = ReadGmsh(FINLEY_TEST_MESH_PATH+os.sep+file,2,optimize=False)           dom = ReadGmsh(os.path.join(FINLEY_TEST_MESH_PATH,file),2,optimize=False)
195         dom.write(test)           dom.write(test)
196         self.compare(test, FINLEY_TEST_MESH_PATH+os.sep+ref)           self.compare(test, os.path.join(FINLEY_TEST_MESH_PATH,ref))
197    
198     def test_Tri6(self):     def test_Tri6(self):
199         file="tri6_gmsh.msh"           file="tri6_gmsh.msh"
200         ref="tri6.fly"           ref="tri6.fly"
201         test = FINLEY_WORKDIR+os.sep+"tri8_test.fly"           test = os.path.join(FINLEY_WORKDIR,"tri6_test.fly")
202         dom = ReadGmsh(FINLEY_TEST_MESH_PATH+os.sep+file,2,optimize=False)           dom = ReadGmsh(os.path.join(FINLEY_TEST_MESH_PATH,file),2,optimize=False)
203         dom.write(test)           dom.write(test)
204         self.compare(test, FINLEY_TEST_MESH_PATH+os.sep+ref)           self.compare(test, os.path.join(FINLEY_TEST_MESH_PATH,ref))
205    
206       def test_Tri6_macro(self):
207             file="tri6_gmsh.msh"
208             ref="tri6_macro.fly"
209             test = os.path.join(FINLEY_WORKDIR,"tri6_macro_test.fly")
210             dom = ReadGmsh(os.path.join(FINLEY_TEST_MESH_PATH,file),2,useMacroElements=True,optimize=False)
211             dom.write(test)
212             self.compare(test, os.path.join(FINLEY_TEST_MESH_PATH,ref))
213    
214     def test_Tet4(self):     def test_Tet4(self):
215         file="tet4_gmsh.msh"           file="tet4_gmsh.msh"
216         ref="tet4.fly"           ref="tet4.fly"
217         test = FINLEY_WORKDIR+os.sep+"tet4_test.fly"           test = os.path.join(FINLEY_WORKDIR,"tet4_test.fly")
218         dom = ReadGmsh(FINLEY_TEST_MESH_PATH+os.sep+file,3,optimize=False)           dom = ReadGmsh(os.path.join(FINLEY_TEST_MESH_PATH,file),3,optimize=False)
219         dom.write(test)           dom.write(test)
220         self.compare(test, FINLEY_TEST_MESH_PATH+os.sep+ref)           self.compare(test, os.path.join(FINLEY_TEST_MESH_PATH,ref))
221    
222     def test_Tet10(self):     def test_Tet10(self):
223         file="tet10_gmsh.msh"           file="tet10_gmsh.msh"
224         ref="tet10.fly"           ref="tet10.fly"
225         test = FINLEY_WORKDIR+os.sep+"tet10_test.fly"           test = os.path.join(FINLEY_WORKDIR,"tet10_test.fly")
226         dom = ReadGmsh(FINLEY_TEST_MESH_PATH+os.sep+file,3,optimize=False)           dom = ReadGmsh(os.path.join(FINLEY_TEST_MESH_PATH,file),3,optimize=False)
227         dom.write(test)           dom.write(test)
228         self.compare(test, FINLEY_TEST_MESH_PATH+os.sep+ref)           self.compare(test, os.path.join(FINLEY_TEST_MESH_PATH,ref))
229    
230       def test_Tet10_macro(self):
231             file="tet10_gmsh.msh"
232             ref="tet10_macro.fly"
233             test = os.path.join(FINLEY_WORKDIR,"tet10_macro_test.fly")
234             dom = ReadGmsh(os.path.join(FINLEY_TEST_MESH_PATH,file),3,useMacroElements=True,optimize=False)
235             dom.write(test)
236             self.compare(test, os.path.join(FINLEY_TEST_MESH_PATH,ref))
237    
238  class Test_Reader(unittest.TestCase):  class Test_Reader(unittest.TestCase):
239     def test_ReadWriteTagNames(self):     def test_ReadWriteTagNames(self):
240         file="hex_2D_order2.msh"         file="hex_2D_order2.msh"
241         test = FINLEY_WORKDIR+os.sep+"test.fly"         test = os.path.join(FINLEY_WORKDIR,"test.fly")
242         dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+file,3,optimize=False)         dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,file),3,optimize=False)
243         insertTagNames(dom,A=1,B=2)         insertTagNames(dom,A=1,B=2)
244         dom.write(test)         dom.write(test)
245         dom2 = ReadMesh(test,3,optimize=False)         dom2 = ReadMesh(test,3,optimize=False)
# Line 222  class Test_Reader(unittest.TestCase): Line 253  class Test_Reader(unittest.TestCase):
253         self.failUnless(dom2.isValidTagName("B"))         self.failUnless(dom2.isValidTagName("B"))
254    
255  class Test_Integration(unittest.TestCase):  class Test_Integration(unittest.TestCase):
256     TOL=EPSILON*100.     TOL=EPSILON*500.
257     def __test_2DQ(self,dom,order):     def __test_2DQ(self,dom,order):
258         x=Function(dom).getX()         x=Function(dom).getX()
259         x_bound=FunctionOnBoundary(dom).getX()         x_bound=FunctionOnBoundary(dom).getX()
# Line 281  class Test_Integration(unittest.TestCase Line 312  class Test_Integration(unittest.TestCase
312               res=integrate(x[0]**i*x[1]**j*x[2]**k)               res=integrate(x[0]**i*x[1]**j*x[2]**k)
313               ref=1./((i+1)*(j+1)*(k+1))               ref=1./((i+1)*(j+1)*(k+1))
314               error=abs(res-ref)/abs(ref)               error=abs(res-ref)/abs(ref)
315               self.failUnless(error<=self.TOL,"integration for order (%s,%s,%s) failed. True value = %s, calculated = %s"%(i,j,k,ref,res))               self.failUnless(error<=self.TOL,"integration for order (%s,%s,%s) failed. True value = %s, calculated = %s (error=%e)"%(i,j,k,ref,res, error))
316    
317               res=integrate(x_bound[0]**i*x_bound[1]**j*x_bound[2]**k)               res=integrate(x_bound[0]**i*x_bound[1]**j*x_bound[2]**k)
318               ref=0               ref=0
# Line 298  class Test_Integration(unittest.TestCase Line 329  class Test_Integration(unittest.TestCase
329               else:               else:
330                  ref += 1./((i+1)*(j+1))                  ref += 1./((i+1)*(j+1))
331               error=abs(res-ref)/abs(ref)               error=abs(res-ref)/abs(ref)
332               self.failUnless(error<=self.TOL,"surface integration for order (%s,%s,%s) failed. True value = %s, calculated = %s"%(i,j,k,ref,res))               self.failUnless(error<=self.TOL,"surface integration for order (%s,%s,%s) failed. True value = %s, calculated = %s (error=%e)"%(i,j,k,ref,res, error))
333    
334     def __test_3DT(self,dom,order,raise_tol=1.):     def __test_3DT(self,dom,order,raise_tol=1.):
335         x=Function(dom).getX()         x=Function(dom).getX()
# Line 310  class Test_Integration(unittest.TestCase Line 341  class Test_Integration(unittest.TestCase
341                  res=integrate(x[0]**i*x[1]**j*x[2]**k)                  res=integrate(x[0]**i*x[1]**j*x[2]**k)
342                  ref=1./((i+1)*(j+1)*(k+1))                  ref=1./((i+1)*(j+1)*(k+1))
343                  error=abs(res-ref)/abs(ref)                  error=abs(res-ref)/abs(ref)
344                  self.failUnless(error<=self.TOL*raise_tol,"integration for order (%s,%s,%s) failed. True value = %s, calculated = %s"%(i,j,k,ref,res))                  self.failUnless(error<=self.TOL*raise_tol,"integration for order (%s,%s,%s) failed. True value = %s, calculated = %s (error=%e)"%(i,j,k,ref,res,error))
345    
346                  res=integrate(x_bound[0]**i*x_bound[1]**j*x_bound[2]**k)                  res=integrate(x_bound[0]**i*x_bound[1]**j*x_bound[2]**k)
347                  ref=0                  ref=0
# Line 327  class Test_Integration(unittest.TestCase Line 358  class Test_Integration(unittest.TestCase
358                  else:                  else:
359                     ref += 1./((i+1)*(j+1))                     ref += 1./((i+1)*(j+1))
360                  error=abs(res-ref)/abs(ref)                  error=abs(res-ref)/abs(ref)
361                  self.failUnless(error<=self.TOL*raise_tol,"surface integration for order (%s,%s,%s) failed. True value = %s, calculated = %s"%(i,j,k,ref,res))                  self.failUnless(error<=self.TOL*raise_tol,"surface integration for order (%s,%s,%s) failed. True value = %s, calculated = %s (error=%e)"%(i,j,k,ref,res,error))
362    
363     def test_hex2D_order1(self):     #===================================================================================================
364        my_dom=Rectangle(1,1,integrationOrder=1)     def test_hex2D_order1_integorder1(self):
365          NE=getMPIRankWorld()
366          my_dom=Rectangle(NE,NE,integrationOrder=1)
367        self.__test_2DQ(my_dom,1)        self.__test_2DQ(my_dom,1)
368     def test_hex2D_order2(self):     def test_hex2D_order1_integorder2(self):
369        my_dom=Rectangle(1,1,integrationOrder=2)        NE=getMPIRankWorld()
370          my_dom=Rectangle(NE,NE,integrationOrder=2)
371        self.__test_2DQ(my_dom,2)        self.__test_2DQ(my_dom,2)
372     def test_hex2D_order3(self):     def test_hex2D_order1_integorder3(self):
373        my_dom=Rectangle(1,1,integrationOrder=3)        NE=getMPIRankWorld()
374          my_dom=Rectangle(NE,NE,integrationOrder=3)
375        self.__test_2DQ(my_dom,3)        self.__test_2DQ(my_dom,3)
376     def test_hex2D_order4(self):     def test_hex2D_order1_integorder4(self):
377        my_dom=Rectangle(1,1,integrationOrder=4)        NE=getMPIRankWorld()
378          my_dom=Rectangle(NE,NE,integrationOrder=4)
379        self.__test_2DQ(my_dom,4)        self.__test_2DQ(my_dom,4)
380     def test_hex2D_order5(self):     def test_hex2D_order1_integorder5(self):
381        my_dom=Rectangle(1,1,integrationOrder=5)        NE=getMPIRankWorld()
382          my_dom=Rectangle(NE,NE,integrationOrder=5)
383        self.__test_2DQ(my_dom,5)        self.__test_2DQ(my_dom,5)
384     def test_hex2D_order6(self):     def test_hex2D_order1_integorder6(self):
385        my_dom=Rectangle(1,1,integrationOrder=6)        NE=getMPIRankWorld()
386          my_dom=Rectangle(NE,NE,integrationOrder=6)
387        self.__test_2DQ(my_dom,6)        self.__test_2DQ(my_dom,6)
388     def test_hex2D_order7(self):     def test_hex2D_order1_integorder7(self):
389        my_dom=Rectangle(1,1,integrationOrder=7)        NE=getMPIRankWorld()
390          my_dom=Rectangle(NE,NE,integrationOrder=7)
391        self.__test_2DQ(my_dom,7)        self.__test_2DQ(my_dom,7)
392     def test_hex2D_order8(self):     def test_hex2D_order1_integorder8(self):
393        my_dom=Rectangle(1,1,integrationOrder=8)        NE=getMPIRankWorld()
394          my_dom=Rectangle(NE,NE,integrationOrder=8)
395        self.__test_2DQ(my_dom,8)        self.__test_2DQ(my_dom,8)
396     def test_hex2D_order9(self):     def test_hex2D_order1_integorder9(self):
397        my_dom=Rectangle(1,1,integrationOrder=9)        NE=getMPIRankWorld()
398          my_dom=Rectangle(NE,NE,integrationOrder=9)
399        self.__test_2DQ(my_dom,9)        self.__test_2DQ(my_dom,9)
400     def test_hex2D_order10(self):     def test_hex2D_order1_integorder10(self):
401        my_dom=Rectangle(1,1,integrationOrder=10)        NE=getMPIRankWorld()
402          my_dom=Rectangle(NE,NE,integrationOrder=10)
403        self.__test_2DQ(my_dom,10)        self.__test_2DQ(my_dom,10)
404       #===================================================================================================
405     def test_Tet2D_order1(self):     def test_hex2D_order2_integorder1(self):
406        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=1)        NE=getMPIRankWorld()
407          my_dom=Rectangle(NE,NE,order=2,integrationOrder=1)
408          self.__test_2DQ(my_dom,1)
409       def test_hex2D_order2_integorder2(self):
410          NE=getMPIRankWorld()
411          my_dom=Rectangle(NE,NE,order=2,integrationOrder=2)
412          self.__test_2DQ(my_dom,2)
413       def test_hex2D_order2_integorder3(self):
414          NE=getMPIRankWorld()
415          my_dom=Rectangle(NE,NE,order=2,integrationOrder=3)
416          self.__test_2DQ(my_dom,3)
417       def test_hex2D_order2_integorder4(self):
418          NE=getMPIRankWorld()
419          my_dom=Rectangle(NE,NE,order=2,integrationOrder=4)
420          self.__test_2DQ(my_dom,4)
421       def test_hex2D_order2_integorder5(self):
422          NE=getMPIRankWorld()
423          my_dom=Rectangle(NE,NE,order=2,integrationOrder=5)
424          self.__test_2DQ(my_dom,5)
425       def test_hex2D_order2_integorder6(self):
426          NE=getMPIRankWorld()
427          my_dom=Rectangle(NE,NE,order=2,integrationOrder=6)
428          self.__test_2DQ(my_dom,6)
429       def test_hex2D_order2_integorder7(self):
430          NE=getMPIRankWorld()
431          my_dom=Rectangle(NE,NE,order=2,integrationOrder=7)
432          self.__test_2DQ(my_dom,7)
433       def test_hex2D_order2_integorder8(self):
434          NE=getMPIRankWorld()
435          my_dom=Rectangle(NE,NE,order=2,integrationOrder=8)
436          self.__test_2DQ(my_dom,8)
437       def test_hex2D_order2_integorder9(self):
438          NE=getMPIRankWorld()
439          my_dom=Rectangle(NE,NE,order=2,integrationOrder=9)
440          self.__test_2DQ(my_dom,9)
441       def test_hex2D_order2_integorder10(self):
442          NE=getMPIRankWorld()
443          my_dom=Rectangle(NE,NE,order=2,integrationOrder=10)
444          self.__test_2DQ(my_dom,10)
445       #===================================================================================================
446       def test_hex2D_macro_integorder1(self):
447          NE=getMPIRankWorld()
448          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=1)
449          self.__test_2DQ(my_dom,1)
450       def test_hex2D_macro_integmacro(self):
451          NE=getMPIRankWorld()
452          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=2)
453          self.__test_2DQ(my_dom,2)
454       def test_hex2D_macro_integorder3(self):
455          NE=getMPIRankWorld()
456          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=3)
457          self.__test_2DQ(my_dom,3)
458       def test_hex2D_macro_integorder4(self):
459          NE=getMPIRankWorld()
460          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=4)
461          self.__test_2DQ(my_dom,4)
462       def test_hex2D_macro_integorder5(self):
463          NE=getMPIRankWorld()
464          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=5)
465          self.__test_2DQ(my_dom,5)
466       def test_hex2D_macro_integorder6(self):
467          NE=getMPIRankWorld()
468          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=6)
469          self.__test_2DQ(my_dom,6)
470       def test_hex2D_macro_integorder7(self):
471          NE=getMPIRankWorld()
472          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=7)
473          self.__test_2DQ(my_dom,7)
474       def test_hex2D_macro_integorder8(self):
475          NE=getMPIRankWorld()
476          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=8)
477          self.__test_2DQ(my_dom,8)
478       def test_hex2D_macro_integorder9(self):
479          NE=getMPIRankWorld()
480          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=9)
481          self.__test_2DQ(my_dom,9)
482       def test_hex2D_macro_integorder10(self):
483          NE=getMPIRankWorld()
484          my_dom=Rectangle(NE,NE,order=-1,integrationOrder=10)
485          self.__test_2DQ(my_dom,10)
486       #===================================================================================================
487       def test_Tet2D_order1_integorder1(self):
488          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=1)
489        self.__test_2DT(my_dom,1)        self.__test_2DT(my_dom,1)
490     def test_Tet2D_order2(self):     def test_Tet2D_order1_integorder2(self):
491        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=2)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=2)
492        self.__test_2DT(my_dom,2)        self.__test_2DT(my_dom,2)
493     def test_Tet2D_order3(self):     def test_Tet2D_order1_integorder3(self):
494        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=3)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=3)
495        self.__test_2DT(my_dom,3)        self.__test_2DT(my_dom,3)
496     def test_Tet2D_order4(self):     def test_Tet2D_order1_integorder4(self):
497        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=4)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=4)
498        self.__test_2DT(my_dom,4)        self.__test_2DT(my_dom,4)
499     def test_Tet2D_order5(self):     def test_Tet2D_order1_integorder5(self):
500        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=5)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=5)
501        self.__test_2DT(my_dom,5)        self.__test_2DT(my_dom,5)
502     def test_Tet2D_order6(self):     def test_Tet2D_order1_integorder6(self):
503        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=6)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=6)
504        self.__test_2DT(my_dom,6)        self.__test_2DT(my_dom,6)
505     def test_Tet2D_order7(self):     def test_Tet2D_order1_integorder7(self):
506        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=7)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=7)
507        self.__test_2DT(my_dom,7)        self.__test_2DT(my_dom,7)
508     def test_Tet2D_order8(self):     def test_Tet2D_order1_integorder8(self):
509        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=8)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=8)
510        self.__test_2DT(my_dom,8,1./sqrt(EPSILON))        self.__test_2DT(my_dom,8,1./sqrt(EPSILON))
511     def test_Tet2D_order9(self):     def test_Tet2D_order1_integorder9(self):
512        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=9)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=9)
513        self.__test_2DT(my_dom,9,1./sqrt(EPSILON))        self.__test_2DT(my_dom,9,1./sqrt(EPSILON))
514     def test_Tet2D_order10(self):     def test_Tet2D_order1_integorder10(self):
515        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tri3.fly",optimize=False,integrationOrder=10)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri3.fly"),optimize=False,integrationOrder=10)
516        self.__test_2DT(my_dom,10)        self.__test_2DT(my_dom,10)
517       #===================================================================================================
518     def test_hex3D_order1(self):     def test_Tet2D_order2_integorder2(self):
519        my_dom=Brick(1,1,1,integrationOrder=1)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=1)
520          self.__test_2DT(my_dom,1)
521       def test_Tet2D_order2_integorder2(self):
522          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=2)
523          self.__test_2DT(my_dom,2)
524       def test_Tet2D_order2_integorder3(self):
525          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=3)
526          self.__test_2DT(my_dom,3)
527       def test_Tet2D_order2_integorder4(self):
528          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=4)
529          self.__test_2DT(my_dom,4)
530       def test_Tet2D_order2_integorder5(self):
531          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=5)
532          self.__test_2DT(my_dom,5)
533       def test_Tet2D_order2_integorder6(self):
534          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=6)
535          self.__test_2DT(my_dom,6)
536       def test_Tet2D_order2_integorder7(self):
537          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=7)
538          self.__test_2DT(my_dom,7)
539       def test_Tet2D_order2_integorder8(self):
540          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=8)
541          self.__test_2DT(my_dom,8,1./sqrt(EPSILON))
542       def test_Tet2D_order2_integorder9(self):
543          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=9)
544          self.__test_2DT(my_dom,9,1./sqrt(EPSILON))
545       def test_Tet2D_order2_integorder10(self):
546          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6.fly"),optimize=False,integrationOrder=10)
547          self.__test_2DT(my_dom,10)
548       #===================================================================================================
549       def test_Tet2D_macro_integmacro(self):
550          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=1)
551          self.__test_2DT(my_dom,1)
552       def test_Tet2D_macro_integmacro(self):
553          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=2)
554          self.__test_2DT(my_dom,2)
555       def test_Tet2D_macro_integorder3(self):
556          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=3)
557          self.__test_2DT(my_dom,3)
558       def test_Tet2D_macro_integorder4(self):
559          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=4)
560          self.__test_2DT(my_dom,4)
561       def test_Tet2D_macro_integorder5(self):
562          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=5)
563          self.__test_2DT(my_dom,5)
564       def test_Tet2D_macro_integorder6(self):
565          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=6)
566          self.__test_2DT(my_dom,6)
567       def test_Tet2D_macro_integorder7(self):
568          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=7)
569          self.__test_2DT(my_dom,7)
570       def test_Tet2D_macro_integorder8(self):
571          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=8)
572          self.__test_2DT(my_dom,8,1./sqrt(EPSILON))
573       def test_Tet2D_macro_integorder9(self):
574          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=9)
575          self.__test_2DT(my_dom,9,1./sqrt(EPSILON))
576       def test_Tet2D_macro_integorder10(self):
577          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tri6_macro.fly"),optimize=False,integrationOrder=10)
578          self.__test_2DT(my_dom,10)
579       #===================================================================================================
580       def test_hex3D_order1_integorder1(self):
581          NE=getMPIRankWorld()
582          my_dom=Brick(NE,NE,NE,integrationOrder=1)
583        self.__test_3DQ(my_dom,1)        self.__test_3DQ(my_dom,1)
584     def test_hex3D_order2(self):     def test_hex3D_order1_integorder2(self):
585        my_dom=Brick(1,1,1,integrationOrder=2)        NE=getMPIRankWorld()
586          my_dom=Brick(NE,NE,NE,integrationOrder=2)
587        self.__test_3DQ(my_dom,2)        self.__test_3DQ(my_dom,2)
588     def test_hex3D_order3(self):     def test_hex3D_order1_integorder3(self):
589        my_dom=Brick(1,1,1,integrationOrder=3)        NE=getMPIRankWorld()
590          my_dom=Brick(NE,NE,NE,integrationOrder=3)
591        self.__test_3DQ(my_dom,3)        self.__test_3DQ(my_dom,3)
592     def test_hex3D_order4(self):     def test_hex3D_order1_integorder4(self):
593        my_dom=Brick(1,1,1,integrationOrder=4)        NE=getMPIRankWorld()
594          my_dom=Brick(NE,NE,NE,integrationOrder=4)
595        self.__test_3DQ(my_dom,4)        self.__test_3DQ(my_dom,4)
596     def test_hex3D_order5(self):     def test_hex3D_order1_integorder5(self):
597        my_dom=Brick(1,1,1,integrationOrder=5)        NE=getMPIRankWorld()
598          my_dom=Brick(NE,NE,NE,integrationOrder=5)
599        self.__test_3DQ(my_dom,5)        self.__test_3DQ(my_dom,5)
600     def test_hex3D_order6(self):     def test_hex3D_order1_integorder6(self):
601        my_dom=Brick(1,1,1,integrationOrder=6)        NE=getMPIRankWorld()
602          my_dom=Brick(NE,NE,NE,integrationOrder=6)
603        self.__test_3DQ(my_dom,6)        self.__test_3DQ(my_dom,6)
604     def test_hex3D_order7(self):     def test_hex3D_order1_integorder7(self):
605        my_dom=Brick(1,1,1,integrationOrder=7)        NE=getMPIRankWorld()
606          my_dom=Brick(NE,NE,NE,integrationOrder=7)
607        self.__test_3DQ(my_dom,7)        self.__test_3DQ(my_dom,7)
608     def test_hex3D_order8(self):     def test_hex3D_order1_integorder8(self):
609        my_dom=Brick(1,1,1,integrationOrder=8)        NE=getMPIRankWorld()
610          my_dom=Brick(NE,NE,NE,integrationOrder=8)
611        self.__test_3DQ(my_dom,8)        self.__test_3DQ(my_dom,8)
612     def test_hex3D_order9(self):     def test_hex3D_order1_integorder9(self):
613        my_dom=Brick(1,1,1,integrationOrder=9)        NE=getMPIRankWorld()
614          my_dom=Brick(NE,NE,NE,integrationOrder=9)
615        self.__test_3DQ(my_dom,9)        self.__test_3DQ(my_dom,9)
616     def test_hex3D_order10(self):     def test_hex3D_order1_integorder10(self):
617        my_dom=Brick(1,1,1,integrationOrder=10)        NE=getMPIRankWorld()
618          my_dom=Brick(NE,NE,NE,integrationOrder=10)
619        self.__test_3DQ(my_dom,10)        self.__test_3DQ(my_dom,10)
620       #===================================================================================================
621     def test_Tet3D_order1(self):     def test_hex3D_order2_integorder2(self):
622        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=1)        NE=getMPIRankWorld()
623          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=1)
624          self.__test_3DQ(my_dom,1)
625       def test_hex3D_order2_integorder2(self):
626          NE=getMPIRankWorld()
627          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=2)
628          self.__test_3DQ(my_dom,2)
629       def test_hex3D_order2_integorder3(self):
630          NE=getMPIRankWorld()
631          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=3)
632          self.__test_3DQ(my_dom,3)
633       def test_hex3D_order2_integorder4(self):
634          NE=getMPIRankWorld()
635          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=4)
636          self.__test_3DQ(my_dom,4)
637       def test_hex3D_order2_integorder5(self):
638          NE=getMPIRankWorld()
639          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=5)
640          self.__test_3DQ(my_dom,5)
641       def test_hex3D_order2_integorder6(self):
642          NE=getMPIRankWorld()
643          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=6)
644          self.__test_3DQ(my_dom,6)
645       def test_hex3D_order2_integorder7(self):
646          NE=getMPIRankWorld()
647          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=7)
648          self.__test_3DQ(my_dom,7)
649       def test_hex3D_order2_integorder8(self):
650          NE=getMPIRankWorld()
651          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=8)
652          self.__test_3DQ(my_dom,8)
653       def test_hex3D_order2_integorder9(self):
654          NE=getMPIRankWorld()
655          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=9)
656          self.__test_3DQ(my_dom,9)
657       def test_hex3D_order2_integorder10(self):
658          NE=getMPIRankWorld()
659          my_dom=Brick(NE,NE,NE,order=2,integrationOrder=10)
660          self.__test_3DQ(my_dom,10)
661       #===================================================================================================
662       def test_hex3D_macro_integmacro(self):
663          NE=getMPIRankWorld()
664          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=1)
665          self.__test_3DQ(my_dom,1)
666       def test_hex3D_macro_integmacro(self):
667          NE=getMPIRankWorld()
668          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=2)
669          self.__test_3DQ(my_dom,2)
670       def test_hex3D_macro_integorder3(self):
671          NE=getMPIRankWorld()
672          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=3)
673          self.__test_3DQ(my_dom,3)
674       def test_hex3D_macro_integorder4(self):
675          NE=getMPIRankWorld()
676          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=4)
677          self.__test_3DQ(my_dom,4)
678       def test_hex3D_macro_integorder5(self):
679          NE=getMPIRankWorld()
680          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=5)
681          self.__test_3DQ(my_dom,5)
682       def test_hex3D_macro_integorder6(self):
683          NE=getMPIRankWorld()
684          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=6)
685          self.__test_3DQ(my_dom,6)
686       def test_hex3D_macro_integorder7(self):
687          NE=getMPIRankWorld()
688          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=7)
689          self.__test_3DQ(my_dom,7)
690       def test_hex3D_macro_integorder8(self):
691          NE=getMPIRankWorld()
692          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=8)
693          self.__test_3DQ(my_dom,8)
694       def test_hex3D_macro_integorder9(self):
695          NE=getMPIRankWorld()
696          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=9)
697          self.__test_3DQ(my_dom,9)
698       def test_hex3D_macro_integorder10(self):
699          NE=getMPIRankWorld()
700          my_dom=Brick(NE,NE,NE,order=-1,integrationOrder=10)
701          self.__test_3DQ(my_dom,10)
702       #==========================================================================================
703       def test_Tet3D_order1_integorder1(self):
704          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=1)
705          self.__test_3DT(my_dom,1)
706       def test_Tet3D_order1_integorder2(self):
707          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=2)
708          self.__test_3DT(my_dom,2)
709       def test_Tet3D_order1_integorder3(self):
710          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=3)
711          self.__test_3DT(my_dom,3)
712       def test_Tet3D_order1_integorder4(self):
713          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=4)
714          self.__test_3DT(my_dom,4)
715       def test_Tet3D_order1_integorder5(self):
716          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=5)
717          self.__test_3DT(my_dom,5)
718       def test_Tet3D_order1_integorder6(self):
719          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=6)
720          self.__test_3DT(my_dom,6,1./sqrt(EPSILON))
721       def test_Tet3D_order1_integorder7(self):
722          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=7)
723          self.__test_3DT(my_dom,7,1./sqrt(EPSILON))
724       def test_Tet3D_order1_integorder8(self):
725          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=8)
726          self.__test_3DT(my_dom,8,1./sqrt(EPSILON))
727       def test_Tet3D_order1_integorder9(self):
728          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=9)
729          self.__test_3DT(my_dom,9,1./sqrt(EPSILON))
730       def test_Tet3D_order1_integorder10(self):
731          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=10)
732          self.__test_3DT(my_dom,10,1./sqrt(EPSILON))
733       #==========================================================================================
734       def test_Tet3D_order2_integorder2(self):
735          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=1)
736        self.__test_3DT(my_dom,1)        self.__test_3DT(my_dom,1)
737     def test_Tet3D_order2(self):     def test_Tet3D_order2_integorder2(self):
738        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=2)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=2)
739        self.__test_3DT(my_dom,2)        self.__test_3DT(my_dom,2)
740     def test_Tet3D_order3(self):     def test_Tet3D_order2_integorder3(self):
741        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=3)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=3)
742        self.__test_3DT(my_dom,3)        self.__test_3DT(my_dom,3)
743     def test_Tet3D_order4(self):     def test_Tet3D_order2_integorder4(self):
744        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=4)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=4)
745        self.__test_3DT(my_dom,4)        self.__test_3DT(my_dom,4)
746     def test_Tet3D_order5(self):     def test_Tet3D_order2_integorder5(self):
747        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=5)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=5)
748        self.__test_3DT(my_dom,5)        self.__test_3DT(my_dom,5)
749     def test_Tet3D_order6(self):     def test_Tet3D_order2_integorder6(self):
750        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=6)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=6)
751        self.__test_3DT(my_dom,6,1./sqrt(EPSILON))        self.__test_3DT(my_dom,6,1./sqrt(EPSILON))
752     def test_Tet3D_order7(self):     def test_Tet3D_order2_integorder7(self):
753        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=7)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=7)
754        self.__test_3DT(my_dom,7,1./sqrt(EPSILON))        self.__test_3DT(my_dom,7,1./sqrt(EPSILON))
755     def test_Tet3D_order8(self):     def test_Tet3D_order2_integorder8(self):
756        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=8)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=8)
757        self.__test_3DT(my_dom,8,1./sqrt(EPSILON))        self.__test_3DT(my_dom,8,1./sqrt(EPSILON))
758     def test_Tet3D_order9(self):     def test_Tet3D_order2_integorder9(self):
759        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=9)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=9)
760        self.__test_3DT(my_dom,9,1./sqrt(EPSILON))        self.__test_3DT(my_dom,9,1./sqrt(EPSILON))
761     def test_Tet3D_order10(self):     def test_Tet3D_order2_integorder10(self):
762        my_dom = ReadMesh(FINLEY_TEST_MESH_PATH+os.sep+"tet4.fly",optimize=False,integrationOrder=10)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=10)
763        self.__test_3DT(my_dom,10)        self.__test_3DT(my_dom,10,1./sqrt(EPSILON))
764       #==========================================================================================
765       def test_Tet3D_macro_integmacro(self):
766          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=1)
767          self.__test_3DT(my_dom,1)
768       def test_Tet3D_macro_integmacro(self):
769          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=2)
770          self.__test_3DT(my_dom,2)
771       def test_Tet3D_macro_integorder3(self):
772          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=3)
773          self.__test_3DT(my_dom,3)
774       def test_Tet3D_macro_integorder4(self):
775          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=4)
776          self.__test_3DT(my_dom,4)
777       def test_Tet3D_macro_integorder5(self):
778          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=5)
779          self.__test_3DT(my_dom,5)
780       def test_Tet3D_macro_integorder6(self):
781          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=6)
782          self.__test_3DT(my_dom,6,1./sqrt(EPSILON))
783       def test_Tet3D_macro_integorder7(self):
784          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=7)
785          self.__test_3DT(my_dom,7,1./sqrt(EPSILON))
786       def test_Tet3D_macro_integorder8(self):
787          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=8)
788          self.__test_3DT(my_dom,8,1./sqrt(EPSILON))
789       def test_Tet3D_macro_integorder9(self):
790          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=9)
791          self.__test_3DT(my_dom,9,1./sqrt(EPSILON))
792       def test_Tet3D_macro_integorder10(self):
793          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10_macro.fly"),optimize=False,integrationOrder=10)
794          self.__test_3DT(my_dom,10,1./sqrt(EPSILON))
795    
796  if __name__ == '__main__':  if __name__ == '__main__':
797     suite = unittest.TestSuite()     suite = unittest.TestSuite()
798     suite.addTest(unittest.makeSuite(Test_Generators))  
799     suite.addTest(unittest.makeSuite(Test_GMSHReader))     if getMPISizeWorld() == 1:
800     suite.addTest(unittest.makeSuite(Test_Reader))          suite.addTest(unittest.makeSuite(Test_Generators))
801       else:
802            print "Test_Generators is dropped as number of processors >1"
803    
804       if getMPISizeWorld() == 1:
805            suite.addTest(unittest.makeSuite(Test_GMSHReader))
806       else:
807            print "Test_GMSHReader is dropped as number of processors >1"
808    
809       if getMPISizeWorld() == 1:
810            suite.addTest(unittest.makeSuite(Test_Reader))
811       else:
812            print "Test_Reader is dropped as number of processors >1"
813    
814     suite.addTest(unittest.makeSuite(Test_Integration))     suite.addTest(unittest.makeSuite(Test_Integration))
815     s=unittest.TextTestRunner(verbosity=2).run(suite)     s=unittest.TextTestRunner(verbosity=2).run(suite)
816     if not s.wasSuccessful(): sys.exit(1)     if not s.wasSuccessful(): sys.exit(1)
   

Legend:
Removed from v.1388  
changed lines
  Added in v.2722

  ViewVC Help
Powered by ViewVC 1.1.26