/[escript]/trunk/doc/user/examples/wave.py
ViewVC logotype

Diff of /trunk/doc/user/examples/wave.py

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

revision 110 by jgs, Mon Feb 14 04:14:42 2005 UTC revision 113 by jgs, Mon Feb 28 07:06:33 2005 UTC
# Line 48  def wavePropagation(domain,h,tend,lam,mu Line 48  def wavePropagation(domain,h,tend,lam,mu
48       print "a=",inf(a),sup(a)       print "a=",inf(a),sup(a)
49       print "u=",inf(u),sup(u)       print "u=",inf(u),sup(u)
50       # ... save current acceleration in units of gravity       # ... save current acceleration in units of gravity
51       if n%10==0 : (length(a)/9.81).saveDX("/tmp/res/u.%i.dx"%(n/10)       if n%10==0: (length(a)/9.81).saveDX("u.%i.dx"%(n/10))
52    
53    print int(width/depth)
54  mydomain=Brick(ne,int(width/depth)*ne,int(width/depth)*ne,l0=depth,l1=width,l2=width)  mydomain=Brick(ne,int(width/depth)*ne,int(width/depth)*ne,l0=depth,l1=width,l2=width)
55  wavePropagation(mydomain,h,tend,lam,mu,rho,s_tt)  wavePropagation(mydomain,h,tend,lam,mu,rho,s_tt)
56    
Legend:
Removed from v.110  
changed lines
  Added in v.113
  ViewVC Help
Powered by ViewVC 1.1.26