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

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

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

revision 104 by jgs, Fri Dec 17 07:43:12 2004 UTC revision 108 by jgs, Thu Jan 27 06:21:59 2005 UTC
# Line 49  def Solve1(mydomain,height): Line 49  def Solve1(mydomain,height):
49               f_p*=x[j]*(x[j]-l[j])               f_p*=x[j]*(x[j]-l[j])
50         f+=f_p         f+=f_p
51    
52      mypde=Poisson(f=f,q=msk)      mypde=Poisson(mydomain)
53        mypde.setValue(f=f,q=msk)
54      u=mypde.getSolution()      u=mypde.getSolution()
55      error=Lsup(u-u_ex)/Lsup(u_ex)      error=Lsup(u-u_ex)/Lsup(u_ex)
56      print "error = ",error      print "error = ",error
# Line 80  def Solve2(mydomain,height): Line 81  def Solve2(mydomain,height):
81            else:            else:
82               f_p*=x[j]*(2*l[j]-x[j])               f_p*=x[j]*(2*l[j]-x[j])
83         f+=f_p         f+=f_p
84      mypde=Poisson(f=f,q=msk)      mypde=Poisson(mydomain)
85        mypde.setValue(f=f,q=msk)
86      u=mypde.getSolution()      u=mypde.getSolution()
87      error=Lsup(u-u_ex)/Lsup(u_ex)      error=Lsup(u-u_ex)/Lsup(u_ex)
88      print "error = ",error      print "error = ",error

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

  ViewVC Help
Powered by ViewVC 1.1.26