[fluka-discuss]: Fwd: FLUKA hangup

From: Zig Hampel-Arias <zig.hampel-arias_at_icecube.wisc.edu>
Date: Thu, 6 Mar 2014 23:52:43 -0600

Dear all,

        Does anyone have any advice or suggestions on how to understand this infinite loop? Thanks,

Zigfried Hampel-Arias
Physics PhD Candidate
UW-Madison

Inicio del mensaje reenviado:

> De: Zig Hampel-Arias <zhampel_at_icecube.wisc.edu>
> Asunto: FLUKA hangup
> Fecha: March 3, 2014 at 4:52:02 PM CST
> Para: fluka-discuss_at_fluka.org
>
> Hello all,
>
> I seem to have an issue with the latest fluka release (2011.2b.5) and it’s use in CORSIKA 7.400.
>
> When running scripts simulating 1000 proton showers on an E^-2 spectrum from 50 GeV to 200 TeV, more than 70% of my jobs run into an infinite loop. I have left these jobs unattended and they can run for days without advancing to the next simulated event.
> According to the GDB, the loop is in the umfadj.f file, around lines 94-102. However, seeing as how umfadj.f is not an included source file, I’m not quite sure how to proceed. I’ve included the printout from using the GDB having run several steps through this loop. I’ve also included the printout from the CORSIKA output, to show that this simulated shower is not particularly high in energy. Indeed, given the same input card, with the same RN seeds, corsika will hang on the same shower (#245). Of course, changing the random seeds does change where the hangup occurs, yet ~70% of my jobs will get stuck in this loop.
> Please provide some guidance how to fix this issue.
>
> Thanks for your help,
>
> Zigfried Hampel-Arias
> Physics PhD Candidate
> UW-Madison
>
> GDB Output:
>
> GNU gdb (GDB) Red Hat Enterprise Linux (7.2-48.el6)
> Copyright (C) 2010 Free Software Foundation, Inc.
> License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>
> This is free software: you are free to change and redistribute it.
> There is NO WARRANTY, to the extent permitted by law. Type "show copying"
> and "show warranty" for details.
> This GDB was configured as "x86_64-redhat-linux-gnu".
> For bug reporting instructions, please see:
> <http://www.gnu.org/software/gdb/bugs/>.
> Attaching to process 28232
> Reading symbols from /data/user/zhampel/hawc/software/corsika-v740/v740f/run/corsika74000Linux_QGSJET_fluka...done.
> Reading symbols from /usr/lib/libnuma.so.1...(no debugging symbols found)...done.
> Loaded symbols for /usr/lib/libnuma.so.1
> Reading symbols from /lib/libpthread.so.0...(no debugging symbols found)...done.
> [Thread debugging using libthread_db enabled]
> Loaded symbols for /lib/libpthread.so.0
> Reading symbols from /opt/pgi/linux86/13.3/lib/libpgc.so...(no debugging symbols found)...done.
> Loaded symbols for /opt/pgi/linux86/13.3/lib/libpgc.so
> Reading symbols from /lib/librt.so.1...(no debugging symbols found)...done.
> Loaded symbols for /lib/librt.so.1
> Reading symbols from /usr/lib/libg2c.so.0...(no debugging symbols found)...done.
> Loaded symbols for /usr/lib/libg2c.so.0
> Reading symbols from /lib/libm.so.6...(no debugging symbols found)...done.
> Loaded symbols for /lib/libm.so.6
> Reading symbols from /lib/libc.so.6...(no debugging symbols found)...done.
> Loaded symbols for /lib/libc.so.6
> Reading symbols from /lib/ld-linux.so.2...(no debugging symbols found)...done.
> Loaded symbols for /lib/ld-linux.so.2
> 0x0818dbe4 in umfadj_ (etcms=0xc52b8e0, amcms=0xc615ee0, pxcms=0xc566260, pycms=0xc5a0be0, pzcms=0xc5db560, nttprt=0xc6508e0, uma=0xc6508a8, umo=0xc650918) at umfadj.f:99
> 99 umfadj.f: No such file or directory.
> in umfadj.f
> Missing separate debuginfos, use: debuginfo-install compat-libf2c-34-3.4.6-19.el6.i686 glibc-2.12-1.132.el6.i686 numactl-2.0.3-9.el6.i686
> (gdb) where
> #0 0x0818dbe4 in umfadj_ (etcms=0xc52b8e0, amcms=0xc615ee0, pxcms=0xc566260, pycms=0xc5a0be0, pzcms=0xc5db560, nttprt=0xc6508e0, uma=0xc6508a8, umo=0xc650918) at umfadj.f:99
> #1 0x080ffcc0 in fixkin_ () at fixkin.f:158
> #2 0x080f8cb3 in eventv_ (ijij=0xc1051ac, pooo=0xc100668, ekee=0xc100660, txx=0xc100658, tyy=0xc100650, tzz=0xc100648, we=0xc100640, imat=0xc1051a8) at eventv.f:2713
> #3 0x080b668e in flulnk () at ./corsika.F:48015
> #4 0x080857c9 in nucint () at ./corsika.F:22519
> #5 0x0804e1b7 in aamain () at ./corsika.F:3243
> #6 0x0804ad4d in main ()
> #7 0x004fad26 in __libc_start_main () from /lib/libc.so.6
> #8 0x0804ac51 in _start ()
> (gdb) stop
> (gdb) setp
> Undefined command: "setp". Try "help".
> (gdb) step
> 102 in umfadj.f
> (gdb) step
> 94 in umfadj.f
> (gdb) step
> 97 in umfadj.f
> (gdb) step
> 99 in umfadj.f
> (gdb) step
> 97 in umfadj.f
> (gdb) step
> 99 in umfadj.f
> (gdb) step
> 98 in umfadj.f
> (gdb) step
> 97 in umfadj.f
> (gdb) step
> 98 in umfadj.f
> (gdb) step
> 99 in umfadj.f
> (gdb) step
> 102 in umfadj.f
> (gdb) step
> 94 in umfadj.f
> (gdb) step
> 97 in umfadj.f
>
>
>
>
>
>
>
> CORSIKA Output:
>
>
> ========== SHOWER NO 245 ===============================================
>
> PRESENT TIME : 03.03.2014 22:35:10 UTC
> AND RANDOM NUMBER GENERATOR AT BEGIN OF EVENT : 245
> SEQUENCE = 1 SEED = 2185549 CALLS = 18146828 BILLIONS = 0
> SEQUENCE = 2 SEED = 2250073 CALLS = 197774857 BILLIONS = 0
> SEQUENCE = 3 SEED = 2332032 CALLS = 0 BILLIONS = 0
> PRIMARY ENERGY = 186.1482057405328 GEV
> PRIMARY ANGLES ARE: THETA = 0.6592 RAD, PHI = 1.0683 RAD
>
> TRACKING STARTS AT FIRST INTERACTION
> PRIMARY PARAMETERS AT FIRST INTERACTION POINT
> 1.400E+01 1.984E+02 7.905E-01 2.950E-01 5.368E-01 1.526E+06 0.000E+00 0.000E+00 0.000E+00
>
> TYPE GAMMA COSTHETA PHIX PHIY HEIGHT TIME X-CM Y-CM GEN/CHI ALEVEL E ON STACK
>
>
>
>
>
>
Received on Fri Mar 07 2014 - 07:55:15 CET

This archive was generated by hypermail 2.3.0 : Fri Mar 07 2014 - 07:55:16 CET