RE: RE: RE: Re: R: RE: Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Dear Francesca,
the bug is already spotted by Alfredo and corrected. However for the respin I am afraid you have to wait at least another week, when Alfredo will be back to CERN.
Best Regards
Vasilis
________________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of francesca.belloni_at_libero.it [francesca.belloni_at_libero.it]
Sent: 22 May 2013 11:02
To: francesca.belloni_at_libero.it; Vasilis Vlachoudis; fasso_at_slac.stanford.edu; fluka-discuss_at_fluka.org
Subject: R: RE: RE: Re: R: RE: Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Dear all,
the run I had launched after adding a value for WHAT(2) in RANDOMIZE,
run a little longer than the others, but in end crashed.
I tried some different values for the WHAT(2). For every value the simulation,
sooner or later, crashes.
What can I do to temporally fix the problem while waiting for the respin,
please?
Thank you,
Francesca
>----Messaggio originale----
>Da: francesca.belloni_at_libero.it
>Data: 16/05/2013 12.04
>A: <Vasilis.Vlachoudis_at_cern.ch>, "fasso_at_slac.stanford.edu"<fasso_at_SLAC.
Stanford.EDU>, "fluka-discuss_at_fluka.org"<fluka-discuss_at_fluka.org>
>Ogg: R: RE: RE: Re: R: RE: Program received signal SIGFPE: Floating-point
exception - erroneous arithmetic operation.
>
>I added a WHAT(2) in randomize and as for now it is working.
>Thank you a lot!!!
>
>Francesca
>
>
>
>
>>----Messaggio originale----
>>Da: Vasilis.Vlachoudis_at_cern.ch
>>Data: 16/05/2013 11.51
>>A: "francesca.belloni_at_libero.it"<francesca.belloni_at_libero.it>, "fasso_at_slac.
>stanford.edu"<fasso_at_slac.stanford.edu>, "fluka-discuss_at_fluka.org"<fluka-
>discuss_at_fluka.org>
>>Ogg: RE: RE: Re: R: RE: Program received signal SIGFPE: Floating-point
>exception - erroneous arithmetic operation.
>>
>>The error is rare so probably changing the random seed will help while
>waiting for the final fix.
>>
>>Best Regards
>>Vasilis
>>
>>________________________________________
>>From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on
>behalf of francesca.belloni_at_libero.it [francesca.belloni_at_libero.it]
>>Sent: 15 May 2013 21:55
>>To: Vasilis Vlachoudis; fasso_at_slac.stanford.edu; fluka-discuss_at_fluka.org
>>Subject: R: RE: Re: R: RE: Program received signal SIGFPE: Floating-point
>exception - erroneous arithmetic operation.
>>
>>Dear Vasilis,
>> thank you a lot.
>>After reading your answer I looked for the routine enionf.f, but it is not
>>available as source
>>(cd $FLUPRO
>> grep -r "enionf" .
>>gives
>>
>>./flutil/usbrea.map:
>>/home/irfulx157/mnt/fbelloni/workspace/FLUKA/FLUKA_2012/libflukahp.a(enionf.
>o)
>>./flutil/usbrea.map:
>>/home/irfulx157/mnt/fbelloni/workspace/FLUKA/FLUKA_2012/libflukahp.a(enion.
o)
>>(enionf_)
>>./flutil/usbrea.map:
>>/home/irfulx157/mnt/fbelloni/workspace/FLUKA/FLUKA_2012/libflukahp.a(enionf.
>o)
>>(cmpips_)
>>./flutil/usbrea.map:
>>/home/irfulx157/mnt/fbelloni/workspace/FLUKA/FLUKA_2012/libflukahp.a(enionf.
>o)
>>(cmscio_)
>>./flutil/usbrea.map:
>>/home/irfulx157/mnt/fbelloni/workspace/FLUKA/FLUKA_2012/libflukahp.a(enionf.
>o)
>>(dedxfl_)
>>.....)
>>
>> What should I change in order to make the simulation to run up to the
>end,
>>please?
>>Does it help if I change the initial seed (that I am not sure I am able to
>>do), please?
>>I ask because a similar simulation (only the geometry is different.It
>includes
>>only the gas and the 10B), runs up to the end with just a couple of
warnings.
>>
>>Thanks a lot!!!
>>Regards,
>>Francesca
>>
>>>----Messaggio originale----
>>>Da: Vasilis.Vlachoudis_at_cern.ch
>>>Data: 15/05/2013 18.43
>>>A: "francesca.belloni_at_libero.it"<francesca.belloni_at_libero.it>, "fasso_at_slac.
>>stanford.edu"<fasso_at_slac.stanford.edu>, "fluka-discuss_at_fluka.org"<fluka-
>>discuss_at_fluka.org>
>>>Ogg: RE: Re: R: RE: Program received signal SIGFPE: Floating-point
>>exception - erroneous arithmetic operation.
>>>
>>>Dear Francesca,
>>>
>>>I could reproduce the error, which appeared to be a floating point
precision
>>problem. Now is in the hands of level 3 support (Alfredo Ferrari) see below
>his
>>reply:
>>>
>>>... it is precision mismatch apparently peculiar of gfortran.
>>>It comes down from the enionf routine which sets the momentum (PLA)
>>>to 0 and the energy (EKIN) as initial kinetic - (total - mass) which
>>>apart rounding should be zero as well. The code has some tolerance on the
>>>rounding but apparently in some cases with gfortran the rounding is larger
>>>than expected, EKIN larger than 'epsilon" and PLA still zero...
>>>
>>>I changed the logic in enionf so that to get *exactly* zero for EKIN as
>>>well, I'll make a respin later
>>>
>>>Best Regards
>>>Vasilis
>>>
>>>________________________________________
>>>From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on
>>behalf of francesca.belloni_at_libero.it [francesca.belloni_at_libero.it]
>>>Sent: 15 May 2013 00:26
>>>To: fasso_at_slac.stanford.edu; fluka-discuss_at_fluka.org
>>>Subject: R: Re: R: RE: Program received signal SIGFPE: Floating-point
>>exception - erroneous arithmetic operation.
>>>
>>>Dear Alberto,
>>> thank you really a lot for your answer, because it corrects
>>mistakes
>>>that I usually do.
>>>I modified the file source.f as you told, but the same mistake still shows
>>up.
>>>In attachment is the new source.f and the *.err, *.log and *.out files.
>>>
>>>Thank you,
>>>best regards,
>>>Francesca
>>>
>>>
>>>
>>>
>>
>>
>>
>>
>
>
>
>
Received on Mon May 27 2013 - 10:56:55 CEST
This archive was generated by hypermail 2.3.0
: Mon May 27 2013 - 10:56:55 CEST