Error with gyre_nad and use_banded = .true.

Bug/problem reports for any of the GYRE executables (gyre_ad, gyre_nad, etc)
Post Reply
ehsan
Posts: 88
Joined: Sun Jun 16, 2013 11:31 am

Error with gyre_nad and use_banded = .true.

Post by ehsan » Tue Sep 03, 2013 1:15 am

Hello Rich,

I previously posted on this 2.5 Msun model. The input model, and the inlist are both attached.
I use gyre_nad with use_banded = .true., and get the following error:

Code: Select all

Scan Mode Search
================

       l    n_pg     n_p     n_g                 Re(omega)                 Im(omega)                       |D|  n_iter
       1     -96       0      96    0.1010334107095588E+00    0.0000000000000000E+00    0.2868200198504395E-12       5
       1     -95       0      95    0.1020024396757933E+00    0.0000000000000000E+00    0.4072925360993182E-12       5
       1     -94       0      94    0.1031808385382246E+00    0.0000000000000000E+00    0.6343933419524848E-12       5
       1     -93       0      93    0.1042019366414670E+00    0.0000000000000000E+00    0.1089227027367001E-12       5
       1     -92       0      92    0.1053993144404898E+00    0.0000000000000000E+00    0.4717216629870301E-12       5
       1     -91       0      91    0.1065174532507716E+00    0.0000000000000000E+00    0.1369557793327378E-11       5
       1     -90       0      90    0.1076884671190581E+00    0.0000000000000000E+00    0.2591394182583956E-12       5
       1     -89       0      89    0.1089037826070608E+00    0.0000000000000000E+00    0.3154708198645855E-12       5
       1     -88       0      88    0.1100476952878871E+00    0.0000000000000000E+00    0.1929570048936593E-12       5
       1     -87       0      87    0.1113220675462222E+00    0.0000000000000000E+00    0.9350208851949989E-12       5
       1     -86       0      86    0.1125046781018162E+00    0.0000000000000000E+00    0.6649230932944280E-12       5
       1     -85       0      85    0.1138980533991844E+00    0.0000000000000000E+00    0.3335861047834109E-12       5
       1     -84       0      84    0.1151761543001349E+00    0.0000000000000000E+00    0.2488762590345333E-12       5
       1     -83       0      83    0.1166384881063459E+00    0.0000000000000000E+00    0.1123509150193240E-11       5
       1     -82       0      82    0.1180156396017830E+00    0.0000000000000000E+00    0.4556006344677055E-12       5
       1     -81       0      81    0.1194956031575932E+00    0.0000000000000000E+00    0.3048415991940241E-13       5
       1     -80       0      80    0.1209805980435179E+00    0.0000000000000000E+00    0.3932233045350215E-12       5
       1     -79       0      79    0.1224764696601710E+00    0.0000000000000000E+00    0.1365194951170305E-12       5
       1     -78       0      78    0.1241209095134377E+00    0.0000000000000000E+00    0.5405485680222143E-13       5
       1     -77       0      77    0.1256217952175543E+00    0.0000000000000000E+00    0.3460959643354208E-12       5
       1     -76       0      76    0.1273538976195898E+00    0.0000000000000000E+00    0.9027012565592134E-12       5
       1     -75       0      75    0.1289041696434456E+00    0.0000000000000000E+00    0.1966500990541240E-12       5
       1     -74       0      74    0.1307015320669421E+00    0.0000000000000000E+00    0.2005917774355025E-12       5
       1     -73       0      73    0.1323821188759239E+00    0.0000000000000000E+00    0.2314816228884164E-12       5
       1     -72       0      72    0.1342573281999979E+00    0.0000000000000000E+00    0.2875263536262997E-12       5
       1     -71       0      71    0.1361454201980188E+00    0.0000000000000000E+00    0.4139485374055941E-13       5
       1     -70       0      70    0.1379928029758126E+00    0.0000000000000000E+00    0.2514975416197152E-12       5
       1     -69       0      69    0.1400926220543450E+00    0.0000000000000000E+00    0.3841298255035267E-16       5
       1     -68       0      68    0.1418983174948014E+00    0.0000000000000000E+00    0.1807423584970049E-13       5
       1     -67       0      67    0.1441826673534052E+00    0.0000000000000000E+00    0.5359599143865390E-13       5
       1     -66       0      66    0.1460190890787532E+00    0.0000000000000000E+00    0.5862547600898894E-13       6
       1     -65       0      65    0.1484466764598170E+00    0.0000000000000000E+00    0.1799318318797885E-13       5
       1     -64       0      64    0.1504260407913539E+00    0.0000000000000000E+00    0.5074583587300739E-13       5
       1     -63       0      63    0.1529844246990274E+00    0.0000000000000000E+00    0.2045258966533510E-12       5
       1     -62       0      62    0.1552145099721096E+00    0.0000000000000000E+00    0.3470499111766122E-12       5
       1     -61       0      61    0.1579603475865887E+00    0.0000000000000000E+00    0.2643479581310385E-13       5
       1     -60       0      60    0.1606030465513400E+00    0.0000000000000000E+00    0.3129306737053830E-13       5
       1     -59       0      59    0.1632342465684345E+00    0.0000000000000000E+00    0.3523055203979516E-12       5
       1     -58       0      58    0.1662693707677630E+00    0.0000000000000000E+00    0.9534631562836093E-12       5
       1     -57       0      57    0.1687221723105393E+00    0.0000000000000000E+00    0.1452647737116862E-12       5
       1     -56       0      56    0.1720592566438506E+00    0.0000000000000000E+00    0.1003043222706699E-11       5
       1     -55       0      55    0.1747183980177924E+00    0.0000000000000000E+00    0.2672793794797909E-12       5
       1     -54       0      54    0.1783671718852502E+00    0.0000000000000000E+00    0.3099348518764686E-12       5
       1     -53       0      53    0.1814403314915445E+00    0.0000000000000000E+00    0.1286674442855791E-12       5
       1     -52       0      52    0.1850973647157437E+00    0.0000000000000000E+00    0.4995704158901601E-12       5
       1     -51       0      51    0.1884200798340558E+00    0.0000000000000000E+00    0.2099664617069536E-12       5
       1     -50       0      50    0.1920602592755871E+00    0.0000000000000000E+00    0.7449526643347541E-13       5
       1     -49       0      49    0.1956046731417846E+00    0.0000000000000000E+00    0.3275490384147385E-12       5
       1     -48       0      48    0.1994632717051332E+00    0.0000000000000000E+00    0.1450195368547549E-12       5
       1     -47       0      47    0.2035610265107031E+00    0.0000000000000000E+00    0.1617788695635283E-13       5
       1     -46       0      46    0.2075362048732808E+00    0.0000000000000000E+00    0.9103957203956975E-14       5
       1     -45       0      45    0.2122421883663053E+00    0.0000000000000000E+00    0.5947789960697779E-13       5
       1     -44       0      44    0.2166091503944479E+00    0.0000000000000000E+00    0.1717152039879002E-12       5
       1     -43       0      43    0.2219573123460899E+00    0.0000000000000000E+00    0.2535126132140070E-12       5
       1     -42       0      42    0.2272019879931686E+00    0.0000000000000000E+00    0.9011709603020188E-13       5
       1     -41       0      41    0.2328398991108290E+00    0.0000000000000000E+00    0.7008200654639780E-13       5
       1     -40       0      40    0.2384152588738759E+00    0.0000000000000000E+00    0.2472127036650190E-12       5
       1     -39       0      39    0.2447939129531485E+00    0.0000000000000000E+00    0.2526871431912499E-13       5
       1     -38       0      38    0.2509183729291993E+00    0.0000000000000000E+00    0.2258214690422311E-13       5
       1     -37       0      37    0.2576737182409021E+00    0.0000000000000000E+00    0.1006207464892843E-12       5
       1     -36       0      36    0.2648979060312922E+00    0.0000000000000000E+00    0.4990079674446882E-13       5
       1     -35       0      35    0.2723357848791827E+00    0.0000000000000000E+00    0.1869774535653349E-12       5
       1     -34       0      34    0.2801522077955372E+00    0.0000000000000000E+00    0.8309411217145275E-13       5
       1     -33       0      33    0.2877540401047703E+00    0.0000000000000000E+00    0.8949618294264683E-13       5
       1     -32       0      32    0.2964721079380839E+00    0.0000000000000000E+00    0.2818683290124336E-13       5
       1     -31       0      31    0.3049946249040792E+00    0.0000000000000000E+00    0.8020764710449647E-13       5
       1     -30       0      30    0.3156406235324624E+00    0.0000000000000000E+00    0.1133852633680451E-12       5
       1     -29       0      29    0.3258921904054345E+00    0.0000000000000000E+00    0.3491832320865445E-12       5
       1     -28       0      28    0.3369538123164588E+00    0.0000000000000000E+00    0.2170593941345804E-12       5
       1     -27       0      27    0.3503445573742661E+00    0.0000000000000000E+00    0.2209534764782067E-13       6
       1     -26       0      26    0.3603555511798301E+00    0.0000000000000000E+00    0.8810069169346502E-13       6
       1     -25       0      25    0.3778502024402491E+00    0.0000000000000000E+00    0.5674634987792747E-13       5
       1     -24       0      24    0.3910996161649664E+00    0.0000000000000000E+00    0.2566852498640671E-12       5
       1     -23       0      23    0.4075390410935487E+00    0.0000000000000000E+00    0.1614520687390480E-12       5
       1     -22       0      22    0.4240919798021389E+00    0.0000000000000000E+00    0.1878605910602095E-13       5
       1     -21       0      21    0.4424557385913993E+00    0.0000000000000000E+00    0.4013776820726995E-13       5
       1     -20       0      20    0.4629310415299980E+00    0.0000000000000000E+00    0.4260860369734359E-13       5
       1     -19       0      19    0.4844929603034605E+00    0.0000000000000000E+00    0.7864725466794669E-13       6
       1     -18       0      18    0.5155741826733703E+00    0.0000000000000000E+00    0.1158202831013596E-12       5
       1     -17       0      17    0.5369718431328410E+00    0.0000000000000000E+00    0.3989860200486664E-13       6
       1     -16       0      16    0.5818347419705715E+00    0.0000000000000000E+00    0.7360505825524523E-13       6
       1     -15       0      15    0.6093234552340577E+00    0.0000000000000000E+00    0.1461649881373517E-12       6
       1     -14       0      14    0.6584004493843595E+00    0.0000000000000000E+00    0.5243629948270637E-13       5
       1     -13       0      13    0.7156859816444593E+00    0.0000000000000000E+00    0.8442692891314031E-13       5
       1     -12       0      12    0.7573155565086770E+00    0.0000000000000000E+00    0.3492681100834697E-13       6
       1     -11       0      11    0.8533620629884339E+00    0.0000000000000000E+00    0.1137279318715299E-12       6
       1     -10       0      10    0.9166058666189153E+00    0.0000000000000000E+00    0.2240136858395817E-13       5
       1      -9       0       9    0.1018826371588847E+01    0.0000000000000000E+00    0.1841986903920394E-13       6
       1      -8       0       8    0.1099616864884529E+01    0.0000000000000000E+00    0.1375429382321489E-13       6
       1      -7       0       7    0.1286316920764203E+01    0.0000000000000000E+00    0.4269791749379402E-13       6
       1      -6       0       6    0.1410462340047458E+01    0.0000000000000000E+00    0.6839903165304637E-13       7
       1      -5       0       5    0.1727785042039837E+01    0.0000000000000000E+00    0.5908796548106761E-14       6
       1      -4       0       4    0.2019586602182124E+01    0.0000000000000000E+00    0.1675225559192458E-13       6
       1      -3       0       3    0.2569988394587036E+01    0.0000000000000000E+00    0.4983138202902331E-13       5
       1      -2       0       2    0.3593337974145171E+01    0.0000000000000000E+00    0.1609728046319349E-12       8
       1      -1       0       1    0.3712427183384085E+01    0.0000000000000000E+00    0.4847513563274167E-13       6
       1       1       1       1    0.4740254221251930E+01    0.0000000000000000E+00    0.8139195684436946E-13       5
       1       2       2       1    0.5855951590329085E+01    0.0000000000000000E+00    0.1416329431539033E-13       5
       1       3       3       1    0.6667769226648659E+01    0.0000000000000000E+00    0.3581617709203667E-12       8
       1       4       4       1    0.7271644928292038E+01    0.0000000000000000E+00    0.1401054011053063E-12       5
       1       5       5       1    0.8377471557568768E+01    0.0000000000000000E+00    0.1140640864120857E-12       5
       1       6       5       0    0.9522030696892838E+01    0.0000000000000000E+00    0.2148936914938237E-13       5

Completed mode search; time elapsed:      6.278 s

Proximity Mode Search
=====================

       l    n_pg     n_p     n_g                Re(omega)                Im(omega)                      |D|  n_it
 INFO at line 679 <gyre_sysmtx:null_vector_slu_r>:
 Smallest element not in outer block

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

Backtrace for this error:
#0  0x107fb1092
#1  0x107fb185e
#2  0x7fff95a0fcf9
#3  0x7fff8c681da5
#4  0x107b5f4d7
#5  0x107ab9c28
#6  0x107aae072
#7  0x107bd69c8
Floating point exception: 8
I also noticed that in the p-mode regime, the net order, n_pg is no longer n_p - n_g, but rather has an excess of one order.
Can you please comment on this?

Best,
Ehsan.
Attachments
namelist.txt
Namelist
(1.51 KiB) Downloaded 188 times
model.gyre
Input model
(370.29 KiB) Downloaded 258 times

User avatar
rhtownsend
Site Admin
Posts: 397
Joined: Sun Mar 31, 2013 4:22 pm

Re: Error with gyre_nad and use_banded = .true.

Post by rhtownsend » Tue Sep 03, 2013 7:59 am

ehsan wrote:Hello Rich,

I previously posted on this 2.5 Msun model. The input model, and the inlist are both attached.
I use gyre_nad with use_banded = .true., and get the following error:

Code: Select all

Scan Mode Search
================

       l    n_pg     n_p     n_g                 Re(omega)                 Im(omega)                       |D|  n_iter
       1     -96       0      96    0.1010334107095588E+00    0.0000000000000000E+00    0.2868200198504395E-12       5
       1     -95       0      95    0.1020024396757933E+00    0.0000000000000000E+00    0.4072925360993182E-12       5
       1     -94       0      94    0.1031808385382246E+00    0.0000000000000000E+00    0.6343933419524848E-12       5
       1     -93       0      93    0.1042019366414670E+00    0.0000000000000000E+00    0.1089227027367001E-12       5
       1     -92       0      92    0.1053993144404898E+00    0.0000000000000000E+00    0.4717216629870301E-12       5
       1     -91       0      91    0.1065174532507716E+00    0.0000000000000000E+00    0.1369557793327378E-11       5
       1     -90       0      90    0.1076884671190581E+00    0.0000000000000000E+00    0.2591394182583956E-12       5
       1     -89       0      89    0.1089037826070608E+00    0.0000000000000000E+00    0.3154708198645855E-12       5
       1     -88       0      88    0.1100476952878871E+00    0.0000000000000000E+00    0.1929570048936593E-12       5
       1     -87       0      87    0.1113220675462222E+00    0.0000000000000000E+00    0.9350208851949989E-12       5
       1     -86       0      86    0.1125046781018162E+00    0.0000000000000000E+00    0.6649230932944280E-12       5
       1     -85       0      85    0.1138980533991844E+00    0.0000000000000000E+00    0.3335861047834109E-12       5
       1     -84       0      84    0.1151761543001349E+00    0.0000000000000000E+00    0.2488762590345333E-12       5
       1     -83       0      83    0.1166384881063459E+00    0.0000000000000000E+00    0.1123509150193240E-11       5
       1     -82       0      82    0.1180156396017830E+00    0.0000000000000000E+00    0.4556006344677055E-12       5
       1     -81       0      81    0.1194956031575932E+00    0.0000000000000000E+00    0.3048415991940241E-13       5
       1     -80       0      80    0.1209805980435179E+00    0.0000000000000000E+00    0.3932233045350215E-12       5
       1     -79       0      79    0.1224764696601710E+00    0.0000000000000000E+00    0.1365194951170305E-12       5
       1     -78       0      78    0.1241209095134377E+00    0.0000000000000000E+00    0.5405485680222143E-13       5
       1     -77       0      77    0.1256217952175543E+00    0.0000000000000000E+00    0.3460959643354208E-12       5
       1     -76       0      76    0.1273538976195898E+00    0.0000000000000000E+00    0.9027012565592134E-12       5
       1     -75       0      75    0.1289041696434456E+00    0.0000000000000000E+00    0.1966500990541240E-12       5
       1     -74       0      74    0.1307015320669421E+00    0.0000000000000000E+00    0.2005917774355025E-12       5
       1     -73       0      73    0.1323821188759239E+00    0.0000000000000000E+00    0.2314816228884164E-12       5
       1     -72       0      72    0.1342573281999979E+00    0.0000000000000000E+00    0.2875263536262997E-12       5
       1     -71       0      71    0.1361454201980188E+00    0.0000000000000000E+00    0.4139485374055941E-13       5
       1     -70       0      70    0.1379928029758126E+00    0.0000000000000000E+00    0.2514975416197152E-12       5
       1     -69       0      69    0.1400926220543450E+00    0.0000000000000000E+00    0.3841298255035267E-16       5
       1     -68       0      68    0.1418983174948014E+00    0.0000000000000000E+00    0.1807423584970049E-13       5
       1     -67       0      67    0.1441826673534052E+00    0.0000000000000000E+00    0.5359599143865390E-13       5
       1     -66       0      66    0.1460190890787532E+00    0.0000000000000000E+00    0.5862547600898894E-13       6
       1     -65       0      65    0.1484466764598170E+00    0.0000000000000000E+00    0.1799318318797885E-13       5
       1     -64       0      64    0.1504260407913539E+00    0.0000000000000000E+00    0.5074583587300739E-13       5
       1     -63       0      63    0.1529844246990274E+00    0.0000000000000000E+00    0.2045258966533510E-12       5
       1     -62       0      62    0.1552145099721096E+00    0.0000000000000000E+00    0.3470499111766122E-12       5
       1     -61       0      61    0.1579603475865887E+00    0.0000000000000000E+00    0.2643479581310385E-13       5
       1     -60       0      60    0.1606030465513400E+00    0.0000000000000000E+00    0.3129306737053830E-13       5
       1     -59       0      59    0.1632342465684345E+00    0.0000000000000000E+00    0.3523055203979516E-12       5
       1     -58       0      58    0.1662693707677630E+00    0.0000000000000000E+00    0.9534631562836093E-12       5
       1     -57       0      57    0.1687221723105393E+00    0.0000000000000000E+00    0.1452647737116862E-12       5
       1     -56       0      56    0.1720592566438506E+00    0.0000000000000000E+00    0.1003043222706699E-11       5
       1     -55       0      55    0.1747183980177924E+00    0.0000000000000000E+00    0.2672793794797909E-12       5
       1     -54       0      54    0.1783671718852502E+00    0.0000000000000000E+00    0.3099348518764686E-12       5
       1     -53       0      53    0.1814403314915445E+00    0.0000000000000000E+00    0.1286674442855791E-12       5
       1     -52       0      52    0.1850973647157437E+00    0.0000000000000000E+00    0.4995704158901601E-12       5
       1     -51       0      51    0.1884200798340558E+00    0.0000000000000000E+00    0.2099664617069536E-12       5
       1     -50       0      50    0.1920602592755871E+00    0.0000000000000000E+00    0.7449526643347541E-13       5
       1     -49       0      49    0.1956046731417846E+00    0.0000000000000000E+00    0.3275490384147385E-12       5
       1     -48       0      48    0.1994632717051332E+00    0.0000000000000000E+00    0.1450195368547549E-12       5
       1     -47       0      47    0.2035610265107031E+00    0.0000000000000000E+00    0.1617788695635283E-13       5
       1     -46       0      46    0.2075362048732808E+00    0.0000000000000000E+00    0.9103957203956975E-14       5
       1     -45       0      45    0.2122421883663053E+00    0.0000000000000000E+00    0.5947789960697779E-13       5
       1     -44       0      44    0.2166091503944479E+00    0.0000000000000000E+00    0.1717152039879002E-12       5
       1     -43       0      43    0.2219573123460899E+00    0.0000000000000000E+00    0.2535126132140070E-12       5
       1     -42       0      42    0.2272019879931686E+00    0.0000000000000000E+00    0.9011709603020188E-13       5
       1     -41       0      41    0.2328398991108290E+00    0.0000000000000000E+00    0.7008200654639780E-13       5
       1     -40       0      40    0.2384152588738759E+00    0.0000000000000000E+00    0.2472127036650190E-12       5
       1     -39       0      39    0.2447939129531485E+00    0.0000000000000000E+00    0.2526871431912499E-13       5
       1     -38       0      38    0.2509183729291993E+00    0.0000000000000000E+00    0.2258214690422311E-13       5
       1     -37       0      37    0.2576737182409021E+00    0.0000000000000000E+00    0.1006207464892843E-12       5
       1     -36       0      36    0.2648979060312922E+00    0.0000000000000000E+00    0.4990079674446882E-13       5
       1     -35       0      35    0.2723357848791827E+00    0.0000000000000000E+00    0.1869774535653349E-12       5
       1     -34       0      34    0.2801522077955372E+00    0.0000000000000000E+00    0.8309411217145275E-13       5
       1     -33       0      33    0.2877540401047703E+00    0.0000000000000000E+00    0.8949618294264683E-13       5
       1     -32       0      32    0.2964721079380839E+00    0.0000000000000000E+00    0.2818683290124336E-13       5
       1     -31       0      31    0.3049946249040792E+00    0.0000000000000000E+00    0.8020764710449647E-13       5
       1     -30       0      30    0.3156406235324624E+00    0.0000000000000000E+00    0.1133852633680451E-12       5
       1     -29       0      29    0.3258921904054345E+00    0.0000000000000000E+00    0.3491832320865445E-12       5
       1     -28       0      28    0.3369538123164588E+00    0.0000000000000000E+00    0.2170593941345804E-12       5
       1     -27       0      27    0.3503445573742661E+00    0.0000000000000000E+00    0.2209534764782067E-13       6
       1     -26       0      26    0.3603555511798301E+00    0.0000000000000000E+00    0.8810069169346502E-13       6
       1     -25       0      25    0.3778502024402491E+00    0.0000000000000000E+00    0.5674634987792747E-13       5
       1     -24       0      24    0.3910996161649664E+00    0.0000000000000000E+00    0.2566852498640671E-12       5
       1     -23       0      23    0.4075390410935487E+00    0.0000000000000000E+00    0.1614520687390480E-12       5
       1     -22       0      22    0.4240919798021389E+00    0.0000000000000000E+00    0.1878605910602095E-13       5
       1     -21       0      21    0.4424557385913993E+00    0.0000000000000000E+00    0.4013776820726995E-13       5
       1     -20       0      20    0.4629310415299980E+00    0.0000000000000000E+00    0.4260860369734359E-13       5
       1     -19       0      19    0.4844929603034605E+00    0.0000000000000000E+00    0.7864725466794669E-13       6
       1     -18       0      18    0.5155741826733703E+00    0.0000000000000000E+00    0.1158202831013596E-12       5
       1     -17       0      17    0.5369718431328410E+00    0.0000000000000000E+00    0.3989860200486664E-13       6
       1     -16       0      16    0.5818347419705715E+00    0.0000000000000000E+00    0.7360505825524523E-13       6
       1     -15       0      15    0.6093234552340577E+00    0.0000000000000000E+00    0.1461649881373517E-12       6
       1     -14       0      14    0.6584004493843595E+00    0.0000000000000000E+00    0.5243629948270637E-13       5
       1     -13       0      13    0.7156859816444593E+00    0.0000000000000000E+00    0.8442692891314031E-13       5
       1     -12       0      12    0.7573155565086770E+00    0.0000000000000000E+00    0.3492681100834697E-13       6
       1     -11       0      11    0.8533620629884339E+00    0.0000000000000000E+00    0.1137279318715299E-12       6
       1     -10       0      10    0.9166058666189153E+00    0.0000000000000000E+00    0.2240136858395817E-13       5
       1      -9       0       9    0.1018826371588847E+01    0.0000000000000000E+00    0.1841986903920394E-13       6
       1      -8       0       8    0.1099616864884529E+01    0.0000000000000000E+00    0.1375429382321489E-13       6
       1      -7       0       7    0.1286316920764203E+01    0.0000000000000000E+00    0.4269791749379402E-13       6
       1      -6       0       6    0.1410462340047458E+01    0.0000000000000000E+00    0.6839903165304637E-13       7
       1      -5       0       5    0.1727785042039837E+01    0.0000000000000000E+00    0.5908796548106761E-14       6
       1      -4       0       4    0.2019586602182124E+01    0.0000000000000000E+00    0.1675225559192458E-13       6
       1      -3       0       3    0.2569988394587036E+01    0.0000000000000000E+00    0.4983138202902331E-13       5
       1      -2       0       2    0.3593337974145171E+01    0.0000000000000000E+00    0.1609728046319349E-12       8
       1      -1       0       1    0.3712427183384085E+01    0.0000000000000000E+00    0.4847513563274167E-13       6
       1       1       1       1    0.4740254221251930E+01    0.0000000000000000E+00    0.8139195684436946E-13       5
       1       2       2       1    0.5855951590329085E+01    0.0000000000000000E+00    0.1416329431539033E-13       5
       1       3       3       1    0.6667769226648659E+01    0.0000000000000000E+00    0.3581617709203667E-12       8
       1       4       4       1    0.7271644928292038E+01    0.0000000000000000E+00    0.1401054011053063E-12       5
       1       5       5       1    0.8377471557568768E+01    0.0000000000000000E+00    0.1140640864120857E-12       5
       1       6       5       0    0.9522030696892838E+01    0.0000000000000000E+00    0.2148936914938237E-13       5

Completed mode search; time elapsed:      6.278 s

Proximity Mode Search
=====================

       l    n_pg     n_p     n_g                Re(omega)                Im(omega)                      |D|  n_it
 INFO at line 679 <gyre_sysmtx:null_vector_slu_r>:
 Smallest element not in outer block

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

Backtrace for this error:
#0  0x107fb1092
#1  0x107fb185e
#2  0x7fff95a0fcf9
#3  0x7fff8c681da5
#4  0x107b5f4d7
#5  0x107ab9c28
#6  0x107aae072
#7  0x107bd69c8
Floating point exception: 8
I also noticed that in the p-mode regime, the net order, n_pg is no longer n_p - n_g, but rather has an excess of one order.
Can you please comment on this?

Best,
Ehsan.
Thanks for the bug report; I'm looking into the issue today.

Regarding the radial order, this behavior is correct; in the Takata classification scheme for dipole modes, n_pg = n_p - n_g + 1 when n_p > n_g. This is why the radial order jumps from -1 to +1.

cheers,

Rich

Post Reply