Re: Fluka crashes for higher primaries for heavy ion

From: Jean-Emmanuel Groetz <jegroetz_at_univ-fcomte.fr>
Date: Tue, 23 Oct 2012 14:20:32 +0200

Users and experts,

For the last post of Mohammad concerning the problem with heavy ion, I
tried to run his input file, and I confirm that Fluka stops after a
while:
- after ~18000 primaries for NEW-DEFA
- after ~80000 primaries for HADROTHE

It seems that the flukadpm3 module crashed. It was generated with the
ldpmqmd shell script from the latest version of Fluka, i.e. 2011.2.15,
with the option FLUFOR=gfortran (gfortran 4.6.3).

Here is the error message in the log file:

Program received signal 8 (SIGFPE): Floating-point exception.

Backtrace for this error:
   + /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fe2d3a594a0]
   + /lib/x86_64-linux-gnu/libm.so.6(log+0x24) [0x7fe2d401e544]
   + function sgdedx_ (0xC756D5)
     at line 50 of file sgdedx.f
   + function electr_ (0xC467CD)
     at line 6 of file electr.f
   + function kasemf_ (0xAC3A6F)
     at line 183 of file kasemf.f
   + function kaskad_ (0x91240A)
     at line 3717 of file kaskad.f
   + function feeder_ (0x907725)
     at line 257 of file feeder.f
   + function flukam_ (0x8388C5)
     at line 3095 of file flukam.f
   + function fluka (0x40B8CC)
     at line 289 of file fluka.f
   + /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed)
[0x7fe2d3a4476d]
Quit (core dumped)

-- 
_______________________________________________
Jean-Emmanuel Groetz
Laboratoire Chrono-Environnement UMR CNRS 6249
Chimie-Physique & Rayonnements
Université de Franche-Comté
16 route de Gray
25030 Besançon Cedex
France
mailto : jegroetz_at_univ-fcomte.fr
tel : +33 (0) 3 81 66 65 07
fax : +33 (0) 3 81 66 65 22
_______________________________________________
Le dimanche 21 octobre 2012   19:15 -0500, Mohammad Islam a écrit :
> Dear fluka users,
>
>
> First of all thank you very much Jean-Emmanuel for your suggestion. I
> still have this problem.
>
>
> While I installed fluka, I followed the instructions as mentioned in
> the fluka website. I used the term $FLUPRO/flutil/ldpmqmd
> which I believe generates the flukadpm3 executable. I thought this
> executable enables nucleus-nucleus interactions and that's why I used
> flukadpm3 in falir as executable in this problem. In this setup, I was
> able to run for 5000, 10,000 or 20,000 primaries but it was random.
> Sometimes it crushes even for 2,000 primaries. I absolutely have no
> idea what is happening. Attached is the error file for 50000
> primaries .Someone please help me as I am not much familiar with
> fortran and linux. BTW, I am using fedora 15.
>
>
> Thank you very much.
>
> Sincerely
> Mohammad
> Dept of Physics, Oklahoma State Unversity
> USA
>
>
> On Fri, Oct 19, 2012 at 9:18 AM, Jean-Emmanuel Groetz
> <jegroetz_at_univ-fcomte.fr> wrote:
>         Mohammad,
>
>         On my workstation, I ran your input with 10000 primaries: see
>         the
>         attached plot.
>         I linked the DPMJET/RQMD modules with flair (Link ->ldpmqmd
>         and Exe
>         ->flukadpm).
>         Don't forget to first compile the ldpmqmd file located in the
>         flutil
>         directory.
>
>         Hope this helps.
>
>         Jean-Emmanuel
>
>         --
>         _______________________________________________
>         Jean-Emmanuel Groetz
>         Laboratoire Chrono-Environnement UMR CNRS 6249
>         Chimie-Physique & Rayonnements
>         Université de Franche-Comté
>         16 route de Gray
>         25030 Besançon Cedex
>         France
>
>         mailto : jegroetz_at_univ-fcomte.fr
>         tel : +33 (0) 3 81 66 65 07
>         fax : +33 (0) 3 81 66 65 22
>         _______________________________________________
>
>
>         Le jeudi 18 octobre 2012   21:04 -0500, Mohammad Islam a
>         écrit :
>         >
>         > Dear Fluka experts,
>         >
>         > I was trying to get an LET spectrum from 956 MeV/n Fe beam
>         behind
>         > aluminium shield using an gas filled detector. For heavy
>         ion, I used
>         > flukadpm3 executables (I am not sure if I am right). Using
>         FLAIR, I
>         > can run this input for 1000 primaries but if I increase the
>         primaries
>         > (e.g. 10,000) it finishes with error. If I don't use any
>         executables
>         > this input runs fine even for higher primaries but that is
>         not what I
>         > want. I tried this in several computers but I did not
>         succeed. Could
>         > any one please suggest me how can I run this input for
>         higher number
>         > of primaries (10^6 at least). I  Vastly appreciate any kind
>         of help.
>         > Attached is the input.
>         >
>         > Sincerely
>         > Mohammad
>         > Physics dept, Oklahoma State University
>         >  USA
>         >
>         >
>
>
>
Received on Tue Oct 23 2012 - 17:27:31 CEST

This archive was generated by hypermail 2.2.0 : Tue Oct 23 2012 - 17:27:32 CEST