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

revision 4820 by jfenwick, Thu Feb 6 06:12:20 2014 UTC revision 4821 by sshaw, Tue Apr 1 04:58:33 2014 UTC
# Line 42  import numpy as np Line 42  import numpy as np
42
43  import pylab as pl  import pylab as pl
44  import matplotlib.cm as cm  import matplotlib.cm as cm
45  from esys.escript.linearPDEs import LinearPDE  from esys.escript.linearPDEs import LinearPDE, SolverOptions
46
47  ########################################################MPI WORLD CHECK  ########################################################MPI WORLD CHECK
48  if getMPISizeWorld() > 1:  if getMPISizeWorld() > 1:
# Line 127  x=domain.getX() # get the locations of t Line 127  x=domain.getX() # get the locations of t
127  mypde=LinearPDE(domain) # create pde  mypde=LinearPDE(domain) # create pde
128  mypde.setSymmetryOn() # turn symmetry on  mypde.setSymmetryOn() # turn symmetry on
129  # turn lumping on for more efficient solving  # turn lumping on for more efficient solving
130  mypde.getSolverOptions().setSolverMethod(mypde.getSolverOptions().HRZ_LUMPING)  mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
131  kmat = kronecker(domain) # create the kronecker delta function of the domain  kmat = kronecker(domain) # create the kronecker delta function of the domain
132  mypde.setValue(D=kmat*rho) #set the general form value D  mypde.setValue(D=kmat*rho) #set the general form value D
133
# Line 186  src_dir=numpy.array([0.,1.]) # defines d Line 186  src_dir=numpy.array([0.,1.]) # defines d
186  y=y*src_dir  y=y*src_dir
187  mypde.setValue(y=y) #set the source as a function on the boundary  mypde.setValue(y=y) #set the source as a function on the boundary
188  # turn lumping on for more efficient solving  # turn lumping on for more efficient solving
189  mypde.getSolverOptions().setSolverMethod(mypde.getSolverOptions().HRZ_LUMPING)  mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
190  # for first two time steps  # for first two time steps
191  u=[0.0,0.0]*wherePositive(x)  u=[0.0,0.0]*wherePositive(x)
192  u_m1=u  u_m1=u

Legend:
 Removed from v.4820 changed lines Added in v.4821