# Diff of /trunk/doc/examples/cookbook/example11b.py

revision 3981 by jfenwick, Fri Sep 21 02:47:54 2012 UTC revision 4087 by caltinay, Thu Nov 22 22:28:01 2012 UTC
# Line 25  Author: Antony Hallam antony.hallam@uqco Line 25  Author: Antony Hallam antony.hallam@uqco
25  """  """
26
28  # example10a.py  # example11b.py
29  # Model of gravitational Potential.  # Model of gravitational Potential.
30
31  #######################################################EXTERNAL MODULES  #######################################################EXTERNAL MODULES
32  # To solve the problem it is necessary to import the modules we require.  # To solve the problem it is necessary to import the modules we require.
33  from esys.escript import * # This imports everything from the escript library  import matplotlib
34  from esys.escript.unitsSI import *  matplotlib.use('agg') #It's just here for automated testing
from esys.escript.linearPDEs import LinearPDE # This defines LinearPDE as LinearPDE
from esys.weipa import saveVTK
from esys.finley import Rectangle # This imports the rectangle domain function from finley
35  import os, sys #This package is necessary to handle saving our data.  import os, sys #This package is necessary to handle saving our data.
36  from math import pi, sqrt, sin, cos  from math import pi, sqrt, sin, cos
37    import numpy as np
38    import pylab as pl #Plotting package
39
40    from esys.escript import * # This imports everything from the escript library
41    from esys.escript.linearPDEs import LinearPDE # This defines LinearPDE as LinearPDE
42  from esys.escript.pdetools import Projector  from esys.escript.pdetools import Projector
43    from esys.escript.unitsSI import *
44  import matplotlib  from esys.finley import Rectangle # This imports the rectangle domain function from finley
45  matplotlib.use('agg') #It's just here for automated testing  from esys.weipa import saveVTK
46
47  from cblib import toRegGrid  from cblib import toRegGrid
48

import pylab as pl #Plotting package
import numpy as np

49  ########################################################MPI WORLD CHECK  ########################################################MPI WORLD CHECK
50  if getMPISizeWorld() > 1:  if getMPISizeWorld() > 1:
51      import sys      import sys

Legend:
 Removed from v.3981 changed lines Added in v.4087