# HG changeset patch # User jwe # Date 1068661829 0 # Node ID 1ed6cfbc8ea4054f3ba7b87b93292ebc2f3ec359 # Parent 3680bb30de0b17a87f5d36152ba58e1fa1d64f1f [project @ 2003-11-12 18:30:29 by jwe] diff -r 3680bb30de0b -r 1ed6cfbc8ea4 liboctave/lo-ieee.cc --- a/liboctave/lo-ieee.cc Wed Nov 12 18:25:53 2003 +0000 +++ b/liboctave/lo-ieee.cc Wed Nov 12 18:30:29 2003 +0000 @@ -51,24 +51,34 @@ void octave_ieee_init (void) { + oct_mach_info::float_format ff = oct_mach_info::native_float_format (); + + if (ff == octave_mach_info::flt_fmt_vax_d + || ff == octave_mach_info::flt_fmt_vax_g + || ff == octave_mach_info::flt_fmt_cray) + { + octave_Inf = octave_NaN = octaveNA = DBL_MAX; + } + else + { #if defined (HAVE_ISINF) || defined (HAVE_FINITE) #if defined (SCO) - double tmp = 1.0; - octave_Inf = 1.0 / (tmp - tmp); + double tmp = 1.0; + octave_Inf = 1.0 / (tmp - tmp); #elif defined (__alpha__) && defined (__osf__) - extern unsigned int DINFINITY[2]; - octave_Inf = (*(X_CAST(double *, DINFINITY))); + extern unsigned int DINFINITY[2]; + octave_Inf = (*(X_CAST(double *, DINFINITY))); #else - double tmp = 1e+10; - octave_Inf = tmp; - for (;;) - { - octave_Inf *= 1e+10; - if (octave_Inf == tmp) - break; - tmp = octave_Inf; - } + double tmp = 1e+10; + octave_Inf = tmp; + for (;;) + { + octave_Inf *= 1e+10; + if (octave_Inf == tmp) + break; + tmp = octave_Inf; + } #endif #endif @@ -76,34 +86,33 @@ #if defined (HAVE_ISNAN) #if defined (__alpha__) && defined (__osf__) - extern unsigned int DQNAN[2]; - octave_NaN = (*(X_CAST(double *, DQNAN))); + extern unsigned int DQNAN[2]; + octave_NaN = (*(X_CAST(double *, DQNAN))); #else - octave_NaN = octave_Inf / octave_Inf; + octave_NaN = octave_Inf / octave_Inf; #endif - // This is patterned after code in R. - - oct_mach_info::float_format ff = oct_mach_info::native_float_format (); + // This is patterned after code in R. - if (ff == oct_mach_info::flt_fmt_ieee_big_endian) - { - lo_ieee_hw = 0; - lo_ieee_lw = 1; - } - else - { - lo_ieee_hw = 1; - lo_ieee_lw = 0; - } + if (ff == oct_mach_info::flt_fmt_ieee_big_endian) + { + lo_ieee_hw = 0; + lo_ieee_lw = 1; + } + else + { + lo_ieee_hw = 1; + lo_ieee_lw = 0; + } - lo_ieee_double t; - t.word[lo_ieee_hw] = LO_IEEE_NA_HW; - t.word[lo_ieee_lw] = LO_IEEE_NA_LW; + lo_ieee_double t; + t.word[lo_ieee_hw] = LO_IEEE_NA_HW; + t.word[lo_ieee_lw] = LO_IEEE_NA_LW; - octave_NA = t.value; + octave_NA = t.value; #endif + } } /*