2nd resolution Wave_t_ex
jil at rpi.edu
Wed Jun 6 15:03:50 CEST 2001
I did not see the big warning about the second scheme in your last email.
Now, I read about it. But, I still don't quite understand how to assemble the
source. I do express TimeFct and dfdt in the same order as in your
TimeFct = yo;
// TimeFct = ( (X<.03) )? yo : 0 ;
dfdt = TimeFct ;
SaveFct = !($TimeStep % 2) ;
instead of 'TimeFct = dfdt;'. It seems I have to set 'TimeFct =
dfdt/dfdt[t==t0]' in order to make 'dfdt[t==t0]*TimeFct' the correct
source as you explained in the email. Does that mean dfdt[t==t0] can not be
zero? Please help.
Dept. of Math., RPI
More information about the getdp