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

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