/[escript]/trunk/finley/test/python/FCT_test0.py
ViewVC logotype

Diff of /trunk/finley/test/python/FCT_test0.py

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

revision 3824 by gross, Mon Feb 6 06:17:47 2012 UTC revision 3825 by gross, Tue Feb 14 06:11:58 2012 UTC
# Line 26  from esys.finley import Rectangle, Brick Line 26  from esys.finley import Rectangle, Brick
26  from esys.weipa import saveVTK  from esys.weipa import saveVTK
27  from math import pi, ceil  from math import pi, ceil
28    
29  NE=5  NE=50
30  dom=Rectangle(NE,1,l1=1./NE)  dom=Rectangle(NE,1,l1=1./NE)
31    dom=Rectangle(NE,NE)
32  fc=TransportPDE(dom,numEquations=1)  fc=TransportPDE(dom,numEquations=1)
33  fc.getSolverOptions().setVerbosityOn()  fc.getSolverOptions().setVerbosityOn()
34    fc.getSolverOptions().setODESolver(fc.getSolverOptions().LINEAR_CRANK_NICOLSON)
35  fc.getSolverOptions().setODESolver(fc.getSolverOptions().BACKWARD_EULER)  fc.getSolverOptions().setODESolver(fc.getSolverOptions().BACKWARD_EULER)
36  fc.getSolverOptions().setODESolver(fc.getSolverOptions().CRANK_NICOLSON)  fc.getSolverOptions().setODESolver(fc.getSolverOptions().CRANK_NICOLSON)
 fc.getSolverOptions().setODESolver(fc.getSolverOptions().LINEAR_CRANK_NICOLSON)  
37  fc.setValue(M=1,C=[-1,0])  fc.setValue(M=1,C=[-1,0])
38  x=dom.getX()  x=dom.getX()
39  u0=whereNegative(x[0]-1./NE)  u0=whereNegative(x[0]-1./NE)

Legend:
Removed from v.3824  
changed lines
  Added in v.3825

  ViewVC Help
Powered by ViewVC 1.1.26