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

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

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

revision 3075 by ahallam, Wed Jul 28 02:51:20 2010 UTC revision 3089 by ahallam, Mon Aug 9 07:20:58 2010 UTC
# Line 112  for it in range(0,ls): Line 112  for it in range(0,ls):
112      tt = t-t0      tt = t-t0
113      dum1 = np.exp(-a * tt * tt)      dum1 = np.exp(-a * tt * tt)
114      source[it] = -2. * a * tt * dum1      source[it] = -2. * a * tt * dum1
 #   source[it] = exp(-a * tt * tt)    !gaussian  
115      if (abs(source[it]) > ampmax):      if (abs(source[it]) > ampmax):
116          ampmax = abs(source[it])          ampmax = abs(source[it])
     #source[t]=np.exp(g*t)*U0*np.sin(2.*np.pi*t/(0.75*ls))*(np.exp(-.1*g*t)-1)  
     #decay1[t]=np.exp(g*t)  
     #decay2[t]=(np.exp(-.1*g*t)-1)  
117      time[t]=t*h      time[t]=t*h
118    
119  ####################################################DOMAIN CONSTRUCTION  ####################################################DOMAIN CONSTRUCTION
# Line 196  mypde.setValue(D=rho*kmat) #set the gene Line 192  mypde.setValue(D=rho*kmat) #set the gene
192    
193  ##########################################################ESTABLISH ABC  ##########################################################ESTABLISH ABC
194  # Define where the boundary decay will be applied.  # Define where the boundary decay will be applied.
195  bn=50.  bn=20.
196  bleft=xstep*bn; bright=width-(xstep*bn); bbot=depth-(ystep*bn)  bleft=xstep*bn; bright=width-(xstep*bn); bbot=depth-(ystep*bn)
197  # btop=ystep*bn # don't apply to force boundary!!!  # btop=ystep*bn # don't apply to force boundary!!!
198    

Legend:
Removed from v.3075  
changed lines
  Added in v.3089

  ViewVC Help
Powered by ViewVC 1.1.26