njoy / NJOY2016

Nuclear data processing with legacy NJOY
https://www.njoy21.io/NJOY2016
Other
96 stars 86 forks source link

Segmentation fault in 'error' module for Fe-54 ENDF file #152

Closed jlconlin closed 3 years ago

jlconlin commented 4 years ago

I’m trying to process multigroup covariance matrices of Fe-54. When I’m reading the ENDF file from ENDF8.0, NJOY21 crashes with a simple segmentation fault. I haven’t been able to reproduce the segmentation fault with ENDF7.1 or JEFF3.3 ENDF files with the same NJOY input.

I've included the NJOY input deck and the tape20

jlconlin commented 4 years ago

Unfortunately, this errors in different modules depending on whether you have compiled in debug mode or in release mode. You can see from captured output below that NJOY fails in ERRORR in release mode, but in GROUPR in debug mode.

release mode

<Lots of stuff here>
 ---message from lunion---xsec nonzero at threshold for mt=633
                          adjusted using jump in xsec

 broadr...                                                                1.5s

 groupr...                                                                2.4s

 moder...                                                                 3.4s

 errorr...                                                                7.5s
At line 56 of file /Users/jlconlin/NJOY21/Code/NJOY2016/src/endf.f90 (unit = 91, file = 'tape91')
Fortran runtime error: Bad value during floating point read

debug mode

 ---message from lunion---xsec nonzero at threshold for mt=633
                          adjusted using jump in xsec

 broadr...                                                               11.0s

 groupr...                                                               14.2s

Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
nathangibson14 commented 3 years ago

There are many issues at play in this one, some of which have been identified/solved, so I am going to close it.

Issues here: