/[escript]/trunk/doc/examples/cookbook/example08b.py
ViewVC logotype

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

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

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

  ViewVC Help
Powered by ViewVC 1.1.26