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

Diff of /trunk/doc/examples/cookbook/example09b.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 41  import numpy as np Line 41  import numpy as np
41    
42  import pylab as pl  import pylab as pl
43  import matplotlib.cm as cm  import matplotlib.cm as cm
44  from esys.escript.linearPDEs import LinearPDE  from esys.escript.linearPDEs import LinearPDE, SolverOptions
45  from esys.finley import ReadMesh  from esys.finley import ReadMesh
46    
47  ########################################################MPI WORLD CHECK  ########################################################MPI WORLD CHECK
# Line 138  for i in range(0,nlayers): Line 138  for i in range(0,nlayers):
138  mypde=LinearPDE(domain) # create pde  mypde=LinearPDE(domain) # create pde
139  mypde.setSymmetryOn() # turn symmetry on  mypde.setSymmetryOn() # turn symmetry on
140  # turn lumping on for more efficient solving  # turn lumping on for more efficient solving
141  #mypde.getSolverOptions().setSolverMethod(mypde.getSolverOptions().HRZ_LUMPING)  #mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
142  kmat = kronecker(domain) # create the kronecker delta function of the domain  kmat = kronecker(domain) # create the kronecker delta function of the domain
143  mypde.setValue(D=rho*kmat) #set the general form value D  mypde.setValue(D=rho*kmat) #set the general form value D
144    

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

  ViewVC Help
Powered by ViewVC 1.1.26