/[escript]/trunk/doc/examples/cookbook/twodheatdiff001.py
ViewVC logotype

Diff of /trunk/doc/examples/cookbook/twodheatdiff001.py

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

revision 2401 by ahallam, Wed Apr 29 04:23:07 2009 UTC revision 2534 by caltinay, Thu Jul 16 06:49:19 2009 UTC
# Line 90  mypde.setValue(A=A*kronecker(model),D=D, Line 90  mypde.setValue(A=A*kronecker(model),D=D,
90  # ... set initial temperature ....  # ... set initial temperature ....
91    
92  T= Ti*whereNegative(bound)+Tc*wherePositive(bound) #defining the initial temperatures.  T= Ti*whereNegative(bound)+Tc*wherePositive(bound) #defining the initial temperatures.
93  saveVTK(os.path.join(save_path,"dataedge.xml"), sol=bound)  saveVTK(os.path.join(save_path,"dataedge.vtu"), sol=bound)
94  saveVTK(os.path.join(save_path,"data%03d.xml") %i,sol=T)  saveVTK(os.path.join(save_path,"data%03d.vtu") %i,sol=T)
95    
96  #... start iteration:  #... start iteration:
97  while t<=tend:  while t<=tend:
# Line 100  while t<=tend: Line 100  while t<=tend:
100        Y = T*D        Y = T*D
101        mypde.setValue(Y=Y)        mypde.setValue(Y=Y)
102        T=mypde.getSolution()        T=mypde.getSolution()
103        saveVTK(os.path.join(save_path,"data%03d.xml") %i,sol=T)        saveVTK(os.path.join(save_path,"data%03d.vtu") %i,sol=T)
104    

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

  ViewVC Help
Powered by ViewVC 1.1.26