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

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

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

revision 103 by jgs, Wed Dec 15 07:08:39 2004 UTC revision 104 by jgs, Fri Dec 17 07:43:12 2004 UTC
# Line 2  Line 2 
2  import sys  import sys
3  import os  import os
4  import unittest  import unittest
                                                                                                                                                             
 esys_root=os.getenv('ESYS_ROOT')  
 sys.path.append(esys_root+'/finley/lib')  
 sys.path.append(esys_root+'/escript/lib')  
 sys.path.append(esys_root+'/escript/py_src')  
                                                                                                                                                             
 from escript import *  
 from util import *  
 from linearPDEs import *  
5    
6  import finley  from esys.escript import *
7    from esys.linearPDEs import *
8    from esys.finley import *
9    
10  import numarray  import numarray
11    
12  #  #
# Line 30  def TheTest(msh,constraints,reduce): Line 24  def TheTest(msh,constraints,reduce):
24      u_ex=Scalar(1,what=n)      u_ex=Scalar(1,what=n)
25      for i in range(msh.getDim()):      for i in range(msh.getDim()):
26        u_ex*=sin(2*Pi*x[i])        u_ex*=sin(2*Pi*x[i])
27      mypde=LinearPDE(A=numarray.identity(msh.getDim()),D=sml,Y=(sml+4*Pi**2*msh.getDim())*u_ex,q=constraints,r=u_ex)      mypde=LinearPDE(msh)
28        mypde.setValue(A=numarray.identity(msh.getDim()),D=sml,Y=(sml+4*Pi**2*msh.getDim())*u_ex,q=constraints,r=u_ex)
29      mypde.setSymmetryOn()      mypde.setSymmetryOn()
30      mypde.setDebugOn()      mypde.setDebugOn()
31      mypde.setReducedOrderTo(reduce)      mypde.setReducedOrderTo(reduce)
# Line 55  for onElements in [False,True]: Line 50  for onElements in [False,True]:
50            #  for i1 in [True,False]:            #  for i1 in [True,False]:
51            for i0 in [True,True]:            for i0 in [True,True]:
52              for i1 in [True,True]:              for i1 in [True,True]:
53                msh=finley.Rectangle(numElements,numElements,order,periodic0=i0,periodic1=i1,useElementsOnFace=onElements)                msh=Rectangle(numElements,numElements,order,periodic0=i0,periodic1=i1,useElementsOnFace=onElements)
54                n=ContinuousFunction(msh)                n=ContinuousFunction(msh)
55                x=n.getX()                x=n.getX()
56                c=Scalar(0,what=n)                c=Scalar(0,what=n)
# Line 73  for onElements in [False,True]: Line 68  for onElements in [False,True]:
68            for i0 in [True,False]:            for i0 in [True,False]:
69              for i1 in [True,False]:              for i1 in [True,False]:
70                for i2 in [True,False]:                for i2 in [True,False]:
71                  msh=finley.Brick(numElements,numElements,numElements,order,periodic0=i0,periodic1=i1,periodic2=i2,useElementsOnFace=onElements)                  msh=Brick(numElements,numElements,numElements,order,periodic0=i0,periodic1=i1,periodic2=i2,useElementsOnFace=onElements)
72                  n=ContinuousFunction(msh)                  n=ContinuousFunction(msh)
73                  x=n.getX()                  x=n.getX()
74                  c=Scalar(0,what=n)                  c=Scalar(0,what=n)

Legend:
Removed from v.103  
changed lines
  Added in v.104

  ViewVC Help
Powered by ViewVC 1.1.26