Go to the source code of this file.
Initial value:= sigma1*( -w -vx +
uxx )
wttt = sigma1*( -wtt -vxtt + uxxtt )
wtttt = 0. ! **** sigma1*( -wttt -vxttt + uxxttt )
! wan(
i1,
i2,
i3,
m)=wam(i1,i2,i3,
m)+(2.*dt)*( wt + (dt**2/6.)*wttt )
wa
n(i1,i2,i3,
m)=wa(i1,i2,i3,
m)+(dt)*( wt + dt*( .5*wtt + dt*( (1./6.)*wttt + dt*( (1./24.)*wtttt ) ) ) )
! ********
This file generated from pmlUpdate.
h using pml.
maple *****
c **
NOTE: run pml.
maple to take this file and generate
the pml.
h file
c restart