/[escript]/trunk/doc/examples/geotutorial/backward_euler.py
ViewVC logotype

Diff of /trunk/doc/examples/geotutorial/backward_euler.py

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

revision 2344 by jfenwick, Mon Mar 30 02:13:58 2009 UTC revision 2534 by caltinay, Thu Jul 16 06:49:19 2009 UTC
# Line 58  while t<t_end: Line 58  while t<t_end:
58      mypde.setValue(Y=dt*rhocp*T+dt*Q)      mypde.setValue(Y=dt*rhocp*T+dt*Q)
59      # new temperature:      # new temperature:
60      T=mypde.getSolution()      T=mypde.getSolution()
61      # save to vtk for mayavi:      # save as VTK for visualisation:
62      saveVTK("u.%s.xml"%N,T=T)      saveVTK("u.%s.vtu"%N,T=T)
63      # increase counter and marker:      # increase counter and marker:
64      N+=1; t+=dt      N+=1; t+=dt

Legend:
Removed from v.2344  
changed lines
  Added in v.2534

  ViewVC Help
Powered by ViewVC 1.1.26