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

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

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

revision 1465 by gross, Tue Apr 1 05:20:57 2008 UTC revision 1466 by artak, Wed Apr 2 03:33:43 2008 UTC
# Line 80  while t<T_END: Line 80  while t<T_END:
80      viscosity=exp(A*(1./(1+T.interpolate(Function(dom)))-0.5))      viscosity=exp(A*(1./(1+T.interpolate(Function(dom)))-0.5))
81      print "viscosity range :", inf(viscosity), sup(viscosity)      print "viscosity range :", inf(viscosity), sup(viscosity)
82      sp.initialize(f=T*(RA*unitVector(DIM-1,DIM)),eta=viscosity,fixed_u_mask=fixed_v_mask)      sp.initialize(f=T*(RA*unitVector(DIM-1,DIM)),eta=viscosity,fixed_u_mask=fixed_v_mask)
83      v,p=sp.solve(v,p,show_details=VERBOSE, verbose=True,max_iter=500)      v,p=sp.solve(v,p,show_details=VERBOSE, verbose=True,max_iter=500,solver='GMRES')
84    
85      for d in range(DIM):      for d in range(DIM):
86           print "range %d-velocity"%d,inf(v[d]),sup(v[d])           print "range %d-velocity"%d,inf(v[d]),sup(v[d])

Legend:
Removed from v.1465  
changed lines
  Added in v.1466

  ViewVC Help
Powered by ViewVC 1.1.26