From: Andras V. <and...@ui...> - 2012-04-21 08:29:25
|
Dear Raimar, Great! Thanks a lot for tracking this down. But this means that the first timestep after a continuation was always erroneous so far (that is, if the system time development had *Exact* part)? And this eventually produced an overflow when the trajectory was continued after long enough a time? A little correction is yet necessary, however: if you look at * MCWF_Trajectory<RANK>::coherentTimeDevelopment*, you can see that *tIntPic0_ * is expected to equal the beginning of the present timestep only if *ex_*is true. Otherwise, it is expected to remain zero. (You can find a tentative explanation for this here<http://cppqed.sourceforge.net/structure.html#hamiltonian> .) This condition has to be included here as well. I have merged, corrected and pushed revision #212. Could you please have a look at * MCWF_Trajectory.tcc:128–130* once more and see whether you agree? (You can also remove the raimar/continue_fix branch). Thank again, and best regards, András Dr. Andras Vukics Institute for Theoretical Physics University of Innsbruck |