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

Diff of /trunk/doc/examples/cookbook/example08c.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 os #file path tool Line 40  import os #file path tool
40  from math import * # math package  from math import * # math package
41  from esys.escript import *  from esys.escript import *
42  from esys.escript.unitsSI import *  from esys.escript.unitsSI import *
43  from esys.escript.linearPDEs import LinearPDE  from esys.escript.linearPDEs import LinearPDE, SolverOptions
44  from esys.escript.pdetools import Projector  from esys.escript.pdetools import Projector
45  from esys.weipa import saveVTK  from esys.weipa import saveVTK
46  from cblib import toRegGrid, subsample  from cblib import toRegGrid, subsample
# Line 199  rho.setTaggedValue("bottom",rho2) Line 199  rho.setTaggedValue("bottom",rho2)
199  mypde=LinearPDE(domain) # create pde  mypde=LinearPDE(domain) # create pde
200  mypde.setSymmetryOn() # turn symmetry on  mypde.setSymmetryOn() # turn symmetry on
201  # turn lumping on for more efficient solving  # turn lumping on for more efficient solving
202  #mypde.getSolverOptions().setSolverMethod(mypde.getSolverOptions().LUMPING)  #mypde.getSolverOptions().setSolverMethod(SolverOptions.LUMPING)
203  kmat = kronecker(domain) # create the kronecker delta function of the domain  kmat = kronecker(domain) # create the kronecker delta function of the domain
204  mypde.setValue(D=rho*kmat) #set the general form value D  mypde.setValue(D=rho*kmat) #set the general form value D
205    

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

  ViewVC Help
Powered by ViewVC 1.1.26