Floating point exception

[ehsan]

Hello Rich,
I get the following error with GYRE 3.0.

Code: Select all

gyre_nad [3.0]
==============

Compiler         : GCC version 4.9.0 20131223 (experimental)
Compiler options : -I /Applications/mesasdk/include -march=sandybridge -mmmx -mno-3dnow -msse -msse2 -msse3 -mssse3 -mno-sse4a -mcx16 -msahf -mno-movbe -maes -mpclmul -mpopcnt -mno-abm -mno-lwp -mno-fma -mno-fma4 -mno-xop -mno-bmi -mno-tbm -mno-avx2 -msse4.2 -msse4.1 -mno-lzcnt -mno-rtm -mno-hle -mno-rdrnd -mno-f16c -mno-fsgsbase -mno-rdseed -mno-prfchw -mno-adx -mfxsr -mxsave -mxsaveopt -mno-avx512f -mno-avx512er -mno-avx512cd -mno-avx512pf --param l1-cache-size=32 --param l1-cache-line-size=64 --param l2-cache-size=6144 -mtune=sandybridge -fPIC -mmacosx-version-min=10.7.4 -mno-avx -mno-bmi2 -O2 -std=f2008 -ffpe-summary=none -finit-real=snan -ffpe-trap=invalid,overflow,zero -fbacktrace -fopenmp
OpenMP Threads   : 3
MPI Processors   : 1

Initialization
==============

Reading from MESA file LOGS/M0.80-Z0.2000E-03-00990.gyre
  Detected new-variant file
  Adding central point
Building omega grid
  omega points : 500
  omega range  :   0.5000000000000000E-01 ->   0.1000000000000000E+02

Building x grid
  x points : 903 [after CREATE_CLONE op]
  x range  :  0.0000000000000000E+00 ->   0.9999999999997321E+00

Root bracketing
  Time elapsed :     1.183 s

Root Solving
       l    n_pg     n_p     n_g                 Re(omega)                 Im(omega)                       chi  n_iter        n
       0       1       1       0    0.3168910995191177E+01    0.0000000000000000E+00    0.5459608487790016E-09       6      904
       0       2       2       0    0.4441956022826214E+01    0.0000000000000000E+00    0.2486252112898347E-09       5      904
       0       3       3       0    0.5812441523730745E+01    0.0000000000000000E+00    0.9289675516264947E-10       5      904
       0       4       4       0    0.7103718131769730E+01    0.0000000000000000E+00    0.3439882665091477E-10       5      904
       0       5       5       0    0.8424097565454742E+01    0.0000000000000000E+00    0.3139878367092182E-10       7      904
       0       6       6       0    0.9780817115955132E+01    0.0000000000000000E+00    0.3102075118079560E-10       7      904
  Time elapsed :      0.127 s
Root Solving
       l    n_pg     n_p     n_g                 Re(omega)                 Im(omega)                       chi  n_iter        n

Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:
#0  0x109f570f2
#1  0x109f578ae
#2  0x7fff9355ccf9
#3  0x109a86183
#4  0x1099d3ffa
#5  0x109a73212
#6  0x109a7336d
#7  0x109a72370
#8  0x1099cfbf5
#9  0x1099c7587
#10  0x109b1a0b8
Floating point exception: 8

Can you please capture my mistake?
I have attached my inlist and the input model file

Cheers
Ehsan.

[rhtownsend]

Hi Ehsan --

Thanks for the report. I'm working on a fix; in the meantime, try changing variables_type to DZIEM.

cheers,

Rich