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

[Andras V. <and...@ui...>]

I could already detect the error at T=30 (after several continuations).
There is also a hint of some bug, cf.
http://sourceforge.net/tracker/?func=detail&aid=3512988&group_id=187775&atid=922653


On Thu, Mar 29, 2012 at 9:16 PM, Andras Vukics <and...@ui...>wrote:

> Dear Raimar,
> I've started to experiment, and try to reproduce your error.
> Just a quick question: did you try in debug mode?
> András
>
>
>
> On Wed, Mar 21, 2012 at 6:59 PM, Raimar Sandner <rai...@ui...
> > wrote:
>
>> 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
>>
>> ------------------------------------------------------------------------------
>> This SF email is sponsosred by:
>> Try Windows Azure free for 90 days Click Here
>> http://p.sf.net/sfu/sfd2d-msazure
>> _______________________________________________
>> Cppqed-support mailing list
>> Cpp...@li...
>> https://lists.sourceforge.net/lists/listinfo/cppqed-support
>>
>>
>