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

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

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

revision 3891 by caltinay, Wed Feb 1 06:16:25 2012 UTC revision 3892 by jfenwick, Tue Apr 10 08:57:23 2012 UTC
# Line 38  mypde.setValue(A=kappa*kronecker(mydomai Line 38  mypde.setValue(A=kappa*kronecker(mydomai
38                 d=eta,y=kappa*n[0]+eta*x[0])                 d=eta,y=kappa*n[0]+eta*x[0])
39  #... calculate error of the PDE solution ...  #... calculate error of the PDE solution ...
40  u=mypde.getSolution()  u=mypde.getSolution()
41  print "error is ",Lsup(u-x[0])  print("error is ",Lsup(u-x[0]))
42  # output should be similar to "error is 1.e-7"  # output should be similar to "error is 1.e-7"
43  saveVTK("x0.vtu",sol=u)  saveVTK("x0.vtu",sol=u)
44    

Legend:
Removed from v.3891  
changed lines
  Added in v.3892

  ViewVC Help
Powered by ViewVC 1.1.26