Mercurial > octave-nkf
view liboctave/util/lo-ieee.h @ 18540:89bd70fae066 gui-release
fix initialization problem for Inf, NaN, and NA values (bug #41667)
* lo-cieee.c: Move function definitions to lo-ieee.cc and delete file.
* liboctave/util/module.mk (UTIL_C_SRC): Remove lo-cieee.c from the
list.
* lo-ieee.h (octave_Inf, octave_Float_Inf, octave_NA, octave_Float_NA,
octave_NaN, octave_Float_NaN): Define as macros that expand to
function calls instead of using global variables.
* lo-ieee.cc (lo_inf, lo_nan, lo_na, lo_float_inf, lo_float_nan,
lo_float_na): New static variables.
(octave_ieee_init): Set internal static variables here.
(__lo_ieee_isnan, __lo_ieee_finite, __lo_ieee_isinf,
__lo_ieee_signbit, __lo_ieee_float_isnan, __lo_ieee_float_finite,
__lo_ieee_float_isinf, __lo_ieee_float_signbit):
Use std:: functions if possible, otherwise rely on gnulib.
(lo_ieee_inf_value, lo_ieee_na_value, lo_ieee_nan_value,
lo_ieee_float_inf_value, lo_ieee_float_na_value,
lo_ieee_float_nan_value): Call init function.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Tue, 04 Mar 2014 10:49:31 -0500 |
parents | ebb3ef964372 |
children | 5b263e517c95 |
line wrap: on
line source
/* Copyright (C) 1996-2013 John W. Eaton This file is part of Octave. Octave is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version. Octave is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with Octave; see the file COPYING. If not, see <http://www.gnu.org/licenses/>. */ #if !defined (octave_lo_ieee_h) #define octave_lo_ieee_h 1 #ifdef __cplusplus extern "C" { #endif /* Octave's idea of infinity. */ #define octave_Inf (lo_ieee_inf_value ()) /* Octave's idea of a missing value. */ #define octave_NA (lo_ieee_na_value ()) /* Octave's idea of not a number. */ #define octave_NaN (lo_ieee_nan_value ()) /* Octave's idea of infinity. */ #define octave_Float_Inf (lo_ieee_float_inf_value ()) /* Octave's idea of a missing value. */ #define octave_Float_NA (lo_ieee_float_na_value ()) /* Octave's idea of not a number. */ #define octave_Float_NaN (lo_ieee_float_nan_value ()) /* FIXME -- this code assumes that a double has twice the number of bits as an int */ typedef union { double value; unsigned int word[2]; } lo_ieee_double; typedef union { float value; unsigned int word; } lo_ieee_float; #define LO_IEEE_NA_HW_OLD 0x7ff00000 #define LO_IEEE_NA_LW_OLD 1954 #define LO_IEEE_NA_HW 0x7FF840F4 #define LO_IEEE_NA_LW 0x40000000 #define LO_IEEE_NA_FLOAT 0x7FC207A2 extern OCTAVE_API void octave_ieee_init (void); extern OCTAVE_API int __lo_ieee_isnan (double x); extern OCTAVE_API int __lo_ieee_finite (double x); extern OCTAVE_API int __lo_ieee_isinf (double x); extern OCTAVE_API int __lo_ieee_is_NA (double); extern OCTAVE_API int __lo_ieee_is_old_NA (double); extern OCTAVE_API int __lo_ieee_is_NaN_or_NA (double) GCC_ATTR_DEPRECATED; extern OCTAVE_API double __lo_ieee_replace_old_NA (double); extern OCTAVE_API double lo_ieee_inf_value (void); extern OCTAVE_API double lo_ieee_na_value (void); extern OCTAVE_API double lo_ieee_nan_value (void); extern OCTAVE_API int __lo_ieee_signbit (double); extern OCTAVE_API int __lo_ieee_float_isnan (float x); extern OCTAVE_API int __lo_ieee_float_finite (float x); extern OCTAVE_API int __lo_ieee_float_isinf (float x); extern OCTAVE_API int __lo_ieee_float_is_NA (float); extern OCTAVE_API int __lo_ieee_float_is_NaN_or_NA (float) GCC_ATTR_DEPRECATED; extern OCTAVE_API float lo_ieee_float_inf_value (void); extern OCTAVE_API float lo_ieee_float_na_value (void); extern OCTAVE_API float lo_ieee_float_nan_value (void); extern OCTAVE_API int __lo_ieee_float_signbit (float); #ifdef __cplusplus } #endif #define lo_ieee_isnan(x) (sizeof (x) == sizeof (float) ? \ __lo_ieee_float_isnan (x) : __lo_ieee_isnan (x)) #define lo_ieee_finite(x) (sizeof (x) == sizeof (float) ? \ __lo_ieee_float_finite (x) : __lo_ieee_finite (x)) #define lo_ieee_isinf(x) (sizeof (x) == sizeof (float) ? \ __lo_ieee_float_isinf (x) : __lo_ieee_isinf (x)) #define lo_ieee_is_NA(x) (sizeof (x) == sizeof (float) ? \ __lo_ieee_float_is_NA (x) : __lo_ieee_is_NA (x)) #define lo_ieee_is_NaN_or_NA(x) (sizeof (x) == sizeof (float) ? \ __lo_ieee_float_is_NaN_or_NA (x) : __lo_ieee_is_NaN_or_NA (x)) #define lo_ieee_signbit(x) (sizeof (x) == sizeof (float) ? \ __lo_ieee_float_signbit (x) : __lo_ieee_signbit (x)) #ifdef __cplusplus template <typename T> struct octave_numeric_limits { static T NA (void) { return static_cast<T> (0); } }; template <> struct octave_numeric_limits<double> { static double NA (void) { return octave_NA; } }; template <> struct octave_numeric_limits<float> { static float NA (void) { return octave_Float_NA; } }; #endif #endif