# Diff of /branches/arrayview_from_1695_trunk/doc/user/diffusion.tex

revision 1780 by jfenwick, Mon Aug 11 03:33:40 2008 UTC revision 1781 by jfenwick, Thu Sep 11 05:03:14 2008 UTC
# Line 267  u=mypde.getSolution() Line 267  u=mypde.getSolution()
267  print "error is ",Lsup(u-x[0])  print "error is ",Lsup(u-x[0])
268  saveVTK("x0.xml",sol=u)  saveVTK("x0.xml",sol=u)
269  \end{python}  \end{python}
270  To visualise the solution x0.~xml' just use the command  To visualize the solution x0.~xml' just use the command
271  \begin{python}  \begin{python}
272  mayavi -d u.xml -m SurfaceMap &  mayavi -d u.xml -m SurfaceMap &
273  \end{python}  \end{python}
274  and it is easy to see that the solution $T=x\hackscore{0}$ is calculated.  and it is easy to see that the solution $T=x\hackscore{0}$ is calculated.
275
276  The script is similar to the script \file{poisson.py} dicussed in \Chap{FirstSteps}.  The script is similar to the script \file{poisson.py} discussed in \Chap{FirstSteps}.
277  \code{mydomain.getNormal()} returns the outer normal field on the surface of the domain. The function \function{Lsup}  \code{mydomain.getNormal()} returns the outer normal field on the surface of the domain. The function \function{Lsup}
278  imported by the \code{from esys.escript import *} statement and returns the maximum absolute value of its argument.  imported by the \code{from esys.escript import *} statement and returns the maximum absolute value of its argument.
279  The error shown by the print statement should be in the order of $10^{-7}$. As piecewise bi-linear interpolation is  The error shown by the print statement should be in the order of $10^{-7}$. As piecewise bi-linear interpolation is

Legend:
 Removed from v.1780 changed lines Added in v.1781