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

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