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

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

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

revision 3792 by jfenwick, Wed Jan 18 04:19:33 2012 UTC revision 3793 by gross, Wed Feb 1 07:39:43 2012 UTC
# Line 167  def XXX(dim,tend,dt, s, h,b,c,d,c_dir="x Line 167  def XXX(dim,tend,dt, s, h,b,c,d,c_dir="x
167          print("initial location = ",x0)          print("initial location = ",x0)
168      print("XX", interpolate(uRef(dom,0.,E,s,v,x0), FunctionOnBoundary(dom)))      print("XX", interpolate(uRef(dom,0.,E,s,v,x0), FunctionOnBoundary(dom)))
169            
170      fc_BE=TransportPDE(dom,numEquations=1,useBackwardEuler=True)      fc_BE=TransportPDE(dom,numEquations=1)
171      fc_BE.setValue(M=a, A=-b*kronecker(dom), B=-v_d*a, C=-v_c*a)      fc_BE.setValue(M=a, A=-b*kronecker(dom), B=-v_d*a, C=-v_c*a)
172      fc_BE.getSolverOptions().setVerbosity(VERBOSITY)      fc_BE.getSolverOptions().setVerbosity(VERBOSITY)
173      fc_BE.getSolverOptions().setTolerance(TOL)      fc_BE.getSolverOptions().setTolerance(TOL)
# Line 176  def XXX(dim,tend,dt, s, h,b,c,d,c_dir="x Line 176  def XXX(dim,tend,dt, s, h,b,c,d,c_dir="x
176      fc_BE.getSolverOptions().setNumSweeps(5)        fc_BE.getSolverOptions().setNumSweeps(5)  
177      if VERBOSITY: print("Backward Euler Transport created")      if VERBOSITY: print("Backward Euler Transport created")
178    
179      fc_CN=TransportPDE(dom,numEquations=1,useBackwardEuler=False)      fc_CN=TransportPDE(dom,numEquations=1)
180      fc_CN.setValue(M=a, A=-b*kronecker(dom), B=-v_d*a, C=-v_c*a)      fc_CN.setValue(M=a, A=-b*kronecker(dom), B=-v_d*a, C=-v_c*a)
181      fc_CN.getSolverOptions().setVerbosity(VERBOSITY)      fc_CN.getSolverOptions().setVerbosity(VERBOSITY)
182      fc_CN.getSolverOptions().setTolerance(TOL)      fc_CN.getSolverOptions().setTolerance(TOL)

Legend:
Removed from v.3792  
changed lines
  Added in v.3793

  ViewVC Help
Powered by ViewVC 1.1.26