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

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

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

revision 149 by jgs, Thu Sep 1 03:31:39 2005 UTC revision 153 by jgs, Tue Oct 25 01:51:20 2005 UTC
# Line 1  Line 1 
1  # $Id$  # $Id$
2    
3  import sys  import sys
4  import os  import os
5  import unittest  import unittest
# Line 8  from esys.escript.linearPDEs import * Line 9  from esys.escript.linearPDEs import *
9  from esys import finley  from esys import finley
10    
11  import numarray  import numarray
12    
13  Pi=numarray.pi  Pi=numarray.pi
14    numElements=15
15    sml=0.1
16    
17  #  #
18  # this test the assemblage of problems with periodic boundary conditions:  # this test the assemblage of problems with periodic boundary conditions:
19  #  #
20  numElements=15  
 sml=0.1  
21  #  #
22  #    test solution is u_ex=sin(2*Pi*n*x0)*...*sin(2*Pi*n*x_dim)  #    test solution is u_ex=sin(2*Pi*n*x0)*...*sin(2*Pi*n*x_dim)
23  #  #
# Line 47  for onElements in [False,True]: Line 50  for onElements in [False,True]:
50         redtext=""         redtext=""
51      for dim in [2,3]:      for dim in [2,3]:
52          if dim==2:          if dim==2:
           # for i0 in [True,False]:  
           #  for i1 in [True,False]:  
53            for i0 in [True,True]:            for i0 in [True,True]:
54              for i1 in [True,True]:              for i1 in [True,True]:
55                msh=Rectangle(numElements,numElements,order,periodic0=i0,periodic1=i1,useElementsOnFace=onElements)                msh=finley.Rectangle(numElements,numElements,order,periodic0=i0,periodic1=i1,useElementsOnFace=onElements)
56                n=ContinuousFunction(msh)                n=ContinuousFunction(msh)
57                x=n.getX()                x=n.getX()
58                c=Scalar(0,what=n)                c=Scalar(0,what=n)
# Line 69  for onElements in [False,True]: Line 70  for onElements in [False,True]:
70            for i0 in [True,False]:            for i0 in [True,False]:
71              for i1 in [True,False]:              for i1 in [True,False]:
72                for i2 in [True,False]:                for i2 in [True,False]:
73                  msh=Brick(numElements,numElements,numElements,order,periodic0=i0,periodic1=i1,periodic2=i2,useElementsOnFace=onElements)                  msh=finley.Brick(numElements,numElements,numElements,order,periodic0=i0,periodic1=i1,periodic2=i2,useElementsOnFace=onElements)
74                  n=ContinuousFunction(msh)                  n=ContinuousFunction(msh)
75                  x=n.getX()                  x=n.getX()
76                  c=Scalar(0,what=n)                  c=Scalar(0,what=n)

Legend:
Removed from v.149  
changed lines
  Added in v.153

  ViewVC Help
Powered by ViewVC 1.1.26