/[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 3891 by caltinay, Fri Nov 12 01:19:02 2010 UTC revision 3892 by jfenwick, Tue Apr 10 08:57:23 2012 UTC
# Line 1  Line 1 
1    
2    from __future__ import print_function
3    
4  ########################################################  ########################################################
5  #  #
6  # Copyright (c) 2003-2010 by University of Queensland  # Copyright (c) 2003-2010 by University of Queensland
# Line 54  T=T_D Line 57  T=T_D
57  N=0; t=0  N=0; t=0
58  # stop when t_end is reached:  # stop when t_end is reached:
59  while t<t_end:  while t<t_end:
60      print N,"-th time step t=",t      print(N,"-th time step t=",t)
61      # update PDE coefficient:      # update PDE coefficient:
62      mypde.setValue(Y=dt*rhocp*T+dt*Q)      mypde.setValue(Y=dt*rhocp*T+dt*Q)
63      # new temperature:      # new temperature:

Legend:
Removed from v.3891  
changed lines
  Added in v.3892

  ViewVC Help
Powered by ViewVC 1.1.26