/[escript]/trunk/doc/examples/mount.py
ViewVC logotype

Diff of /trunk/doc/examples/mount.py

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

revision 2305 by artak, Tue Jan 13 03:52:13 2009 UTC revision 2306 by artak, Mon Mar 16 00:23:57 2009 UTC
# Line 6  from math import pi, ceil Line 6  from math import pi, ceil
6  NE=16  NE=16
7  DIM=3  DIM=3
8  H=1.  H=1.
9  L=1*H  L=2*H
10  TOL=1.e-4  TOL=1.e-4
11  OMEGA=10  OMEGA=10
12  EPS=0.01  EPS=0.01
# Line 44  dt=0. Line 44  dt=0.
44  while t<T_END:  while t<T_END:
45      print "STEP ", t      print "STEP ", t
46      u=v*cos(OMEGA*t)      u=v*cos(OMEGA*t)
47      u,Z=mts.update(u=u,H_t=H_t)      u,Z=mts.update(u=u,H_t=H_t,verbose=False)
48            
49      saveVTK("state.%d.vtu"%n,sol=Z)      saveVTK("state.%d.vtu"%n,sol=Z)
50      print "Integral(Z)=",integrate(Z),Lsup(u[DIM-1])      print "Integral(Z)=",integrate(Z),Lsup(u[DIM-1])
Legend:
Removed from v.2305  
changed lines
  Added in v.2306
  ViewVC Help
Powered by ViewVC 1.1.26