/[escript]/trunk/finley/test/python/run_generators.py
ViewVC logotype

Diff of /trunk/finley/test/python/run_generators.py

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

revision 2549 by jfenwick, Mon Jul 20 06:43:47 2009 UTC revision 2722 by gross, Fri Oct 16 06:45:01 2009 UTC
# Line 44  FINLEY_TEST_MESH_PATH=os.path.join(FINLE Line 44  FINLEY_TEST_MESH_PATH=os.path.join(FINLE
44  #   FINLEY_TEST_MESH_PATH = os.path.join(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=os.path.join(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:
# Line 77  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 97  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 188  class Test_GMSHReader(unittest.TestCase) Line 198  class Test_GMSHReader(unittest.TestCase)
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 = os.path.join(FINLEY_WORKDIR,"tri8_test.fly")           test = os.path.join(FINLEY_WORKDIR,"tri6_test.fly")
202           dom = ReadGmsh(os.path.join(FINLEY_TEST_MESH_PATH,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, os.path.join(FINLEY_TEST_MESH_PATH,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"
# Line 209  class Test_GMSHReader(unittest.TestCase) Line 227  class Test_GMSHReader(unittest.TestCase)
227           dom.write(test)           dom.write(test)
228           self.compare(test, os.path.join(FINLEY_TEST_MESH_PATH,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"
# Line 227  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 286  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 303  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 315  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 332  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       def test_hex2D_order1_integorder1(self):
365        NE=getMPIRankWorld()        NE=getMPIRankWorld()
366        my_dom=Rectangle(NE,NE,integrationOrder=1)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
370        my_dom=Rectangle(NE,NE,integrationOrder=2)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
374        my_dom=Rectangle(NE,NE,integrationOrder=3)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
378        my_dom=Rectangle(NE,NE,integrationOrder=4)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
382        my_dom=Rectangle(NE,NE,integrationOrder=5)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
386        my_dom=Rectangle(NE,NE,integrationOrder=6)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
390        my_dom=Rectangle(NE,NE,integrationOrder=7)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
394        my_dom=Rectangle(NE,NE,integrationOrder=8)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
398        my_dom=Rectangle(NE,NE,integrationOrder=9)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
402        my_dom=Rectangle(NE,NE,integrationOrder=10)        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          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)        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(os.path.join(FINLEY_TEST_MESH_PATH,"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(os.path.join(FINLEY_TEST_MESH_PATH,"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(os.path.join(FINLEY_TEST_MESH_PATH,"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(os.path.join(FINLEY_TEST_MESH_PATH,"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(os.path.join(FINLEY_TEST_MESH_PATH,"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(os.path.join(FINLEY_TEST_MESH_PATH,"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(os.path.join(FINLEY_TEST_MESH_PATH,"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(os.path.join(FINLEY_TEST_MESH_PATH,"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(os.path.join(FINLEY_TEST_MESH_PATH,"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 = 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()        NE=getMPIRankWorld()
582        my_dom=Brick(NE,NE,NE,integrationOrder=1)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
586        my_dom=Brick(NE,NE,NE,integrationOrder=2)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
590        my_dom=Brick(NE,NE,NE,integrationOrder=3)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
594        my_dom=Brick(NE,NE,NE,integrationOrder=4)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
598        my_dom=Brick(NE,NE,NE,integrationOrder=5)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
602        my_dom=Brick(NE,NE,NE,integrationOrder=6)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
606        my_dom=Brick(NE,NE,NE,integrationOrder=7)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
610        my_dom=Brick(NE,NE,NE,integrationOrder=8)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
614        my_dom=Brick(NE,NE,NE,integrationOrder=9)        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        NE=getMPIRankWorld()        NE=getMPIRankWorld()
618        my_dom=Brick(NE,NE,NE,integrationOrder=10)        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          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)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=1)
705        self.__test_3DT(my_dom,1)        self.__test_3DT(my_dom,1)
706     def test_Tet3D_order2(self):     def test_Tet3D_order1_integorder2(self):
707        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=2)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=2)
708        self.__test_3DT(my_dom,2)        self.__test_3DT(my_dom,2)
709     def test_Tet3D_order3(self):     def test_Tet3D_order1_integorder3(self):
710        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=3)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=3)
711        self.__test_3DT(my_dom,3)        self.__test_3DT(my_dom,3)
712     def test_Tet3D_order4(self):     def test_Tet3D_order1_integorder4(self):
713        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=4)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=4)
714        self.__test_3DT(my_dom,4)        self.__test_3DT(my_dom,4)
715     def test_Tet3D_order5(self):     def test_Tet3D_order1_integorder5(self):
716        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=5)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=5)
717        self.__test_3DT(my_dom,5)        self.__test_3DT(my_dom,5)
718     def test_Tet3D_order6(self):     def test_Tet3D_order1_integorder6(self):
719        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=6)        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))        self.__test_3DT(my_dom,6,1./sqrt(EPSILON))
721     def test_Tet3D_order7(self):     def test_Tet3D_order1_integorder7(self):
722        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=7)        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))        self.__test_3DT(my_dom,7,1./sqrt(EPSILON))
724     def test_Tet3D_order8(self):     def test_Tet3D_order1_integorder8(self):
725        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=8)        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))        self.__test_3DT(my_dom,8,1./sqrt(EPSILON))
727     def test_Tet3D_order9(self):     def test_Tet3D_order1_integorder9(self):
728        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=9)        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))        self.__test_3DT(my_dom,9,1./sqrt(EPSILON))
730     def test_Tet3D_order10(self):     def test_Tet3D_order1_integorder10(self):
731        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=10)        my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet4.fly"),optimize=False,integrationOrder=10)
732        self.__test_3DT(my_dom,10)        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)
737       def test_Tet3D_order2_integorder2(self):
738          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=2)
739          self.__test_3DT(my_dom,2)
740       def test_Tet3D_order2_integorder3(self):
741          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=3)
742          self.__test_3DT(my_dom,3)
743       def test_Tet3D_order2_integorder4(self):
744          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=4)
745          self.__test_3DT(my_dom,4)
746       def test_Tet3D_order2_integorder5(self):
747          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=5)
748          self.__test_3DT(my_dom,5)
749       def test_Tet3D_order2_integorder6(self):
750          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))
752       def test_Tet3D_order2_integorder7(self):
753          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))
755       def test_Tet3D_order2_integorder8(self):
756          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))
758       def test_Tet3D_order2_integorder9(self):
759          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))
761       def test_Tet3D_order2_integorder10(self):
762          my_dom = ReadMesh(os.path.join(FINLEY_TEST_MESH_PATH,"tet10.fly"),optimize=False,integrationOrder=10)
763          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()

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

  ViewVC Help
Powered by ViewVC 1.1.26