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

Diff of /trunk/doc/examples/geotutorial/forward_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 60  while t<t_end: Line 60  while t<t_end:
60      mypde.setValue(Y=rhocp*T+dt*Q, X=-k*dt*grad(T))      mypde.setValue(Y=rhocp*T+dt*Q, X=-k*dt*grad(T))
61      # new temperature:      # new temperature:
62      T=mypde.getSolution()      T=mypde.getSolution()
63      # save to vtk for mayavi:      # save as VTK for visualisation:
64      # saveVTK("u.%s.xml"%N,T=T)      # saveVTK("u.%s.vtu"%N,T=T)
65      # increase counter and marker:      # increase counter and marker:
66      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