Computation failed with SIGFPE

[joelongjiamian]

When computing frequencies for some red giant models, I receive a SIGFPE signal (Floating-point exception - erroneous arithmetic operation). Any suggestions on how I should go around diagnosing why this is happening? For reference, my input namelists are as follows:

Code: Select all

&model
    model_type = 'EVOL'
    file = '../profile50.data.GYRE'
    file_format = 'MESA'
    ! file_format = 'FGONG'
    ! data_format = '(1P5E16.9,x)'
/

&constants
    G_GRAVITY = 6.67408d-8
/


&mode
    l = 0          ! Harmonic degree
/


&osc
    outer_bound = 'JCD'
    variables_set = 'JCD'
    inertia_norm = 'BOTH'
    !reduce_order = .FALSE.
    rotation_method = 'TAR'
/

&num
    diff_scheme = 'COLLOC_GL4'
/

&scan
    grid_type = 'LINEAR'
    freq_min_units = 'UHZ'
    freq_max_units = 'UHZ'
    freq_min = 0.0071211757296827105
    freq_max = 1
    n_freq = 100
/

&grid
    alpha_osc = 10            ! At least alpha points per oscillatory wavelength
    alpha_exp = 4            ! At least alpha points per exponential 'wavelength'
    n_inner = 20                ! At least n points in the evanescent region
/

&ad_output
    summary_file = 'profile50-freqs-0.dat'
    summary_file_format = 'TXT'
    summary_item_list = 'l,n_pg,n_p,n_g,freq,E_norm'
    freq_units = 'UHZ'
/

The output from GYRE is as follows:

Code: Select all

(expres) [jjo25@mo profile50-freqs-0]$ $GYRE_DIR/bin/gyre l02.in
gyre [5.2]
----------

OpenMP Threads   : 24
Input filename   : l02.in

Model Init
----------

Reading from MESA file
   File name ../profile50.data.GYRE
   File version 1.01
   Read 1975 points
   No need to add central point

Mode Search
-----------

Mode parameters
   l : 0
   m : 0

Building frequency scan
   added scan interval :  0.1362E+00 ->  0.1912E+02 (100 points, LINEAR)

Building x grid
   Found inner turning points, x range 0.0000 -> 0.0000
   Adding 286884 inner point(s)
   Adding 465957 global point(s) in iteration 1
   Adding 172330 global point(s) in iteration 2
   Adding 0 global point(s) in iteration 3
   Final grid has 1 segment(s) and 927146 point(s):
      Segment 1 : x range 0.0000 -> 1.1122 (1 -> 927146)

Starting search (adiabatic)

Root bracketing
 INFO at line 86 <gyre_atmos:atmos_beta_r_>:
 WARNING: Discarding imaginary part of atmospheric radial wavenumber
  Time elapsed :    63.039 s

Root Solving
   l    m    n_pg    n_p    n_g       Re(omega)       Im(omega)        chi n_iter

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

Backtrace for this error:
#0  0x7f9a9b802ebf in ???
#1  0x518bfa in ???
#2  0x519c73 in ???
#3  0x7f9a9be2cdf5 in gomp_thread_start
        at /root/sdk2/build/gcc/libgomp/team.c:123
#4  0x7f9a9b9994bf in ???
#5  0x7f9a9b8c7162 in ???
#6  0xffffffffffffffff in ???
Floating point exception (core dumped)

[rhtownsend]

Can you also post your model?

cheers,

Rich

[joelongjiamian]

Certainly, please see attached.

[rhtownsend]

It looks like you've run into a bug using the traditional approximation of rotation (TAR) with radial modes. Thankfully, the fix is simple: comment out (or remove) the "rotation_method = 'TAR'" line in your input namelist.

I'll fix the bug in the next development release of GYRE.

cheers,

Rich