Mercurial > octave-nkf
annotate 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 |
rev | line source |
---|---|
1967 | 1 /* |
2 | |
17744
d63878346099
maint: Update copyright notices for release.
John W. Eaton <jwe@octave.org>
parents:
15271
diff
changeset
|
3 Copyright (C) 1996-2013 John W. Eaton |
1967 | 4 |
5 This file is part of Octave. | |
6 | |
7 Octave is free software; you can redistribute it and/or modify it | |
8 under the terms of the GNU General Public License as published by the | |
7016 | 9 Free Software Foundation; either version 3 of the License, or (at your |
10 option) any later version. | |
1967 | 11 |
12 Octave is distributed in the hope that it will be useful, but WITHOUT | |
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
15 for more details. | |
16 | |
17 You should have received a copy of the GNU General Public License | |
7016 | 18 along with Octave; see the file COPYING. If not, see |
19 <http://www.gnu.org/licenses/>. | |
1967 | 20 |
21 */ | |
22 | |
17822
ebb3ef964372
maint: Use common #define syntax "octave_filename_h" in h_files.
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
23 #if !defined (octave_lo_ieee_h) |
ebb3ef964372
maint: Use common #define syntax "octave_filename_h" in h_files.
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
24 #define octave_lo_ieee_h 1 |
1967 | 25 |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
26 #ifdef __cplusplus |
4072 | 27 extern "C" { |
28 #endif | |
29 | |
4102 | 30 /* Octave's idea of infinity. */ |
18540
89bd70fae066
fix initialization problem for Inf, NaN, and NA values (bug #41667)
John W. Eaton <jwe@octave.org>
parents:
17822
diff
changeset
|
31 #define octave_Inf (lo_ieee_inf_value ()) |
4102 | 32 |
33 /* Octave's idea of a missing value. */ | |
18540
89bd70fae066
fix initialization problem for Inf, NaN, and NA values (bug #41667)
John W. Eaton <jwe@octave.org>
parents:
17822
diff
changeset
|
34 #define octave_NA (lo_ieee_na_value ()) |
4102 | 35 |
36 /* Octave's idea of not a number. */ | |
18540
89bd70fae066
fix initialization problem for Inf, NaN, and NA values (bug #41667)
John W. Eaton <jwe@octave.org>
parents:
17822
diff
changeset
|
37 #define octave_NaN (lo_ieee_nan_value ()) |
4102 | 38 |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
39 /* Octave's idea of infinity. */ |
18540
89bd70fae066
fix initialization problem for Inf, NaN, and NA values (bug #41667)
John W. Eaton <jwe@octave.org>
parents:
17822
diff
changeset
|
40 #define octave_Float_Inf (lo_ieee_float_inf_value ()) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
41 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
42 /* Octave's idea of a missing value. */ |
18540
89bd70fae066
fix initialization problem for Inf, NaN, and NA values (bug #41667)
John W. Eaton <jwe@octave.org>
parents:
17822
diff
changeset
|
43 #define octave_Float_NA (lo_ieee_float_na_value ()) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
44 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
45 /* Octave's idea of not a number. */ |
18540
89bd70fae066
fix initialization problem for Inf, NaN, and NA values (bug #41667)
John W. Eaton <jwe@octave.org>
parents:
17822
diff
changeset
|
46 #define octave_Float_NaN (lo_ieee_float_nan_value ()) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
47 |
5775 | 48 /* FIXME -- this code assumes that a double has twice the |
5388 | 49 number of bits as an int */ |
50 | |
4075 | 51 typedef union |
52 { | |
53 double value; | |
54 unsigned int word[2]; | |
55 } lo_ieee_double; | |
56 | |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
57 typedef union |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
58 { |
7814
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
59 float value; |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
60 unsigned int word; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
61 } lo_ieee_float; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
62 |
7991
139f47cf17ab
Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents:
7814
diff
changeset
|
63 #define LO_IEEE_NA_HW_OLD 0x7ff00000 |
139f47cf17ab
Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents:
7814
diff
changeset
|
64 #define LO_IEEE_NA_LW_OLD 1954 |
139f47cf17ab
Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents:
7814
diff
changeset
|
65 #define LO_IEEE_NA_HW 0x7FF840F4 |
139f47cf17ab
Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents:
7814
diff
changeset
|
66 #define LO_IEEE_NA_LW 0x40000000 |
139f47cf17ab
Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents:
7814
diff
changeset
|
67 #define LO_IEEE_NA_FLOAT 0x7FC207A2 |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
68 |
6108 | 69 extern OCTAVE_API void octave_ieee_init (void); |
1967 | 70 |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
71 extern OCTAVE_API int __lo_ieee_isnan (double x); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
72 extern OCTAVE_API int __lo_ieee_finite (double x); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
73 extern OCTAVE_API int __lo_ieee_isinf (double x); |
4072 | 74 |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
75 extern OCTAVE_API int __lo_ieee_is_NA (double); |
7991
139f47cf17ab
Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents:
7814
diff
changeset
|
76 extern OCTAVE_API int __lo_ieee_is_old_NA (double); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
77 extern OCTAVE_API int __lo_ieee_is_NaN_or_NA (double) GCC_ATTR_DEPRECATED; |
7991
139f47cf17ab
Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents:
7814
diff
changeset
|
78 extern OCTAVE_API double __lo_ieee_replace_old_NA (double); |
4075 | 79 |
6108 | 80 extern OCTAVE_API double lo_ieee_inf_value (void); |
81 extern OCTAVE_API double lo_ieee_na_value (void); | |
82 extern OCTAVE_API double lo_ieee_nan_value (void); | |
4102 | 83 |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
84 extern OCTAVE_API int __lo_ieee_signbit (double); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
85 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
86 extern OCTAVE_API int __lo_ieee_float_isnan (float x); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
87 extern OCTAVE_API int __lo_ieee_float_finite (float x); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
88 extern OCTAVE_API int __lo_ieee_float_isinf (float x); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
89 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
90 extern OCTAVE_API int __lo_ieee_float_is_NA (float); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
91 extern OCTAVE_API int __lo_ieee_float_is_NaN_or_NA (float) GCC_ATTR_DEPRECATED; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
92 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
93 extern OCTAVE_API float lo_ieee_float_inf_value (void); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
94 extern OCTAVE_API float lo_ieee_float_na_value (void); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
95 extern OCTAVE_API float lo_ieee_float_nan_value (void); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
96 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
97 extern OCTAVE_API int __lo_ieee_float_signbit (float); |
4349 | 98 |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
99 #ifdef __cplusplus |
4072 | 100 } |
2522 | 101 #endif |
102 | |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
103 #define lo_ieee_isnan(x) (sizeof (x) == sizeof (float) ? \ |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
104 __lo_ieee_float_isnan (x) : __lo_ieee_isnan (x)) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
105 #define lo_ieee_finite(x) (sizeof (x) == sizeof (float) ? \ |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
106 __lo_ieee_float_finite (x) : __lo_ieee_finite (x)) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
107 #define lo_ieee_isinf(x) (sizeof (x) == sizeof (float) ? \ |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
108 __lo_ieee_float_isinf (x) : __lo_ieee_isinf (x)) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
109 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
110 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
111 #define lo_ieee_is_NA(x) (sizeof (x) == sizeof (float) ? \ |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
112 __lo_ieee_float_is_NA (x) : __lo_ieee_is_NA (x)) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
113 #define lo_ieee_is_NaN_or_NA(x) (sizeof (x) == sizeof (float) ? \ |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
114 __lo_ieee_float_is_NaN_or_NA (x) : __lo_ieee_is_NaN_or_NA (x)) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
115 #define lo_ieee_signbit(x) (sizeof (x) == sizeof (float) ? \ |
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
116 __lo_ieee_float_signbit (x) : __lo_ieee_signbit (x)) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
117 |
14810
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
118 #ifdef __cplusplus |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
119 |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
120 template <typename T> |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
121 struct octave_numeric_limits |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
122 { |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
123 static T NA (void) { return static_cast<T> (0); } |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
124 }; |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
125 |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
126 template <> |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
127 struct octave_numeric_limits<double> |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
128 { |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
129 static double NA (void) { return octave_NA; } |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
130 }; |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
131 |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
132 template <> |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
133 struct octave_numeric_limits<float> |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
134 { |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
135 static float NA (void) { return octave_Float_NA; } |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
136 }; |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
137 |
1967 | 138 #endif |
14810
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
139 |
d40349493d87
provide template version of octave_NA
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
140 #endif |