------- Comment #20 from dominiq at lps dot ens dot fr 2010-04-14 13:06 -------
I have been able to get the following values before the crash:
DEBav =
0.59252398402327489 1.58848116996742547E-002 2.31157896165751706E-002
8.33002598145726886E-002 9.03564427292446321E-002 -596152.33325048804
-173.27129371256277 284192.62549692224 -226.34926169938268
-2145.2785283028134 6880.2385948782830 25825.135912393274
26997.892028831531 NaN 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
XA =
0.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 1.0000000000000000 1.0000000000000000
1.0000000000000000 0.99999999999999989 1.0000000000000000
1.0000000000000000 1.0000000000000000 0.49470835610659347
NaN NaN 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
XB =
0.0000000000000000 0.0000000000000000 4.95714580495132420E-018
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 7.34677397984116187E-002
NaN NaN 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 0.0000000000000000
ie = I - 1
debe = DEBav(ie)
debs = DEBav(I)
deve = XA(ie)*debe + XB(ie)
devs = XA(I)*debs + XB(I)
sec = SCApa(I)
debm = (deve+devs)/2.D+00
eps = 1.D+00
IF ( debm.LT.0. ) eps = -1.D+00
vis = (r+s*tf)/3600.D+00
Re = eps*debm*DIAhy/(vis*sec)
So if I= 13, 14, or 15, Re is a NaN.
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=43716
