/[escript]/trunk/doc/examples/helmholtz.py
ViewVC logotype

Diff of /trunk/doc/examples/helmholtz.py

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

revision 108 by jgs, Thu Jan 27 06:21:59 2005 UTC revision 113 by jgs, Mon Feb 28 07:06:33 2005 UTC
# Line 7  kappa=1. Line 7  kappa=1.
7  omega=0.1  omega=0.1
8  eta=10.  eta=10.
9  #... generate domain ...  #... generate domain ...
10  mydomain = esys.finley.Rectangle(l0=5.,l1=1.,n0=50, n1=10)  mydomain = Rectangle(l0=5.,l1=1.,n0=50, n1=10)
11  #... open PDE and set coefficients ...  #... open PDE and set coefficients ...
12  mypde=Helmholtz(mydomain)  mypde=Helmholtz(mydomain)
13  n=mydomain.getNormal()  n=mydomain.getNormal()
# Line 16  mypde.setValue(kappa,omega,omega*x[0],et Line 16  mypde.setValue(kappa,omega,omega*x[0],et
16  #... calculate error of the PDE solution ...  #... calculate error of the PDE solution ...
17  u=mypde.getSolution()  u=mypde.getSolution()
18  print "error is ",Lsup(u-x[0])  print "error is ",Lsup(u-x[0])
 # output should be similar to "error is 1.e-7"  
19    # output should be similar to "error is 1.e-7"

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

  ViewVC Help
Powered by ViewVC 1.1.26