/[escript]/trunk/finley/test/python/seismic_wave.py
ViewVC logotype

Diff of /trunk/finley/test/python/seismic_wave.py

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

revision 1964 by ksteube, Thu Sep 25 06:43:44 2008 UTC revision 1965 by caltinay, Wed Nov 5 03:24:09 2008 UTC
# Line 36  from esys.escript.linearPDEs import Line Line 36  from esys.escript.linearPDEs import Line
36  from esys.finley import Brick  from esys.finley import Brick
37  import time  import time
38    
39  WORKDIR="/raid2/lutz/waves/"  WORKDIR="./waves/"
40  output=True  output=True
41  n_end=10000  n_end=10000
42    
# Line 310  def wavePropagation(dom,rho,mu,lmbd,eta) Line 310  def wavePropagation(dom,rho,mu,lmbd,eta)
310     a2=Vector(0.,Solution(dom))     a2=Vector(0.,Solution(dom))
311     v=Vector(0.,Solution(dom))     v=Vector(0.,Solution(dom))
312    
313       if not os.path.isdir(WORKDIR): os.mkdir(WORKDIR)
314    
315     starttime = time.clock()     starttime = time.clock()
316     while t<t_end and n<n_end:     while t<t_end and n<n_end:
317       if output: print n+1,"-th time step t ",t+dt," max u and F: ",Lsup(u),       if output: print n+1,"-th time step t ",t+dt," max u and F: ",Lsup(u),

Legend:
Removed from v.1964  
changed lines
  Added in v.1965

  ViewVC Help
Powered by ViewVC 1.1.26