[Cppqed-support] structure::InfiniteDetectedException when continuing a MCWF trajectory

[Raimar S. <rai...@ui...>]

Dear András,

I tried to continue a MCWF simulation by increasing --T and leaving the output 
file and .sv file in place. The simulation aborted with an 
InfiniteDetectedException.

Indeed: I looked into the average function of the mode, the first displayed 
step is fine, but in the second display step this is the LazyDensityOperator 
which is passed to the average method (std::cout of the corresponding 
DensityOperatorLow):

(0,4) x (0,4)
[ (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan)
  (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan)
  (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan)
  (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan)
  (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan) (-nan,-nan) ]

It is the same for the other frees of the system. Unfortunately I couldn't find 
the reason.

This only happens for longer trajectories, i.e. everything is fine if I 
continue a short trajectory with T=15, but I get the exception with T=120.

For the case that you don't know what is going on right away, I have attached 
the output file and .sv file as a test case. You would need to checkout the 
raimar/private bazaar branch for the script and the interactions I use 
therein. The command line is 

2ParticlesRingCavityFull --etaSin 9 --minitSin 0.1 --modeCavSin Sin \
   --kappaCos 5 --pinit1 "(0.5 0 0.09 0)" --pinit2 "(-0.5 0 0.09 0)" \
   --UnotSin -5 --Dt 1 --cutoffCos 3 --deltaCSin -20 --deltaCCos -20 \
   --evol single --fin1 5 --fin2 5 --precision 6 --dc 0 --minitCos 0.1 \
   --cutoffSin 5 --modeCavCos Cos --UnotCos -5 --kappaSin 5 --seed 1001 \
   --o 2ParticlesRingCavityFull.out.1001 --T 121

Best regards,
Raimar