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

Diff of /trunk/doc/examples/cookbook/example09a.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 40  import numpy as np Line 40  import numpy as np
40    
41  import pylab as pl  import pylab as pl
42  import matplotlib.cm as cm  import matplotlib.cm as cm
43  from esys.escript.linearPDEs import LinearPDE  from esys.escript.linearPDEs import LinearPDE, SolverOptions
44  from esys.finley import ReadMesh  from esys.finley import ReadMesh
45  from esys.weipa import saveVTK  from esys.weipa import saveVTK
46    
# Line 137  rho.setTaggedValue("vintfb",rho2) Line 137  rho.setTaggedValue("vintfb",rho2)
137  mypde=LinearPDE(domain) # create pde  mypde=LinearPDE(domain) # create pde
138  mypde.setSymmetryOn() # turn symmetry on  mypde.setSymmetryOn() # turn symmetry on
139  # turn lumping on for more efficient solving  # turn lumping on for more efficient solving
140  #mypde.getSolverOptions().setSolverMethod(mypde.getSolverOptions().HRZ_LUMPING)  #mypde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
141  kmat = kronecker(domain) # create the kronecker delta function of the domain  kmat = kronecker(domain) # create the kronecker delta function of the domain
142  mypde.setValue(D=rho*kmat) #set the general form value D  mypde.setValue(D=rho*kmat) #set the general form value D
143    

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

  ViewVC Help
Powered by ViewVC 1.1.26