Mercurial > hg > octave-lyh

annotate liboctave/lo-cieee.c @ 10396:a0b51ac0f88a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
optimize accumdim with summation
author Jaroslav Hajek <highegg@gmail.com>
date Fri, 05 Mar 2010 12:31:30 +0100
parents 42d098307c30
children 9900d375e585
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
1 /*
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
2
8920
eb63fbe60fab update copyright notices
John W. Eaton <jwe@octave.org>
parents: 7991
diff changeset
3 Copyright (C) 2002, 2003, 2004, 2005, 2006, 2007, 2008 John W. Eaton
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
4
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
5 This file is part of Octave.
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
6
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
7 Octave is free software; you can redistribute it and/or modify it
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
8 under the terms of the GNU General Public License as published by the
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6061
diff changeset
9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6061
diff changeset
10 option) any later version.
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
11
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
12 Octave is distributed in the hope that it will be useful, but WITHOUT
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
15 for more details.
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
16
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
17 You should have received a copy of the GNU General Public License
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6061
diff changeset
18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6061
diff changeset
19 <http://www.gnu.org/licenses/>.
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
20
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
21 */
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
22
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
23 #ifdef HAVE_CONFIG_H
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
24 #include <config.h>
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
25 #endif
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
26
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
27 #include <float.h>
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
28
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
29 #if defined (HAVE_FLOATINGPOINT_H)
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
30 #include <floatingpoint.h>
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
31 #endif
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
32
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
33 #if defined (HAVE_IEEEFP_H)
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
34 #include <ieeefp.h>
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
35 #endif
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
36
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
37 #if defined (HAVE_NAN_H)
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
38 #if defined (SCO)
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
39 #define _IEEE 1
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
40 #endif
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
41 #include <nan.h>
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
42 #if defined (SCO)
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
43 #undef _IEEE
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
44 #endif
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
45 #endif
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
46
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
47 #include "lo-ieee.h"
7231
2eb392d058bb [project @ 2007-11-30 18:53:29 by jwe]
jwe
parents: 7017
diff changeset
48 #include "lo-math.h"
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
49
6061
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
50 #if ! defined (HAVE_ISNAN) && defined (HAVE__ISNAN)
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
51 #define isnan _isnan
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
52 #define HAVE_ISNAN 1
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
53 #endif
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
54
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
55 #if ! defined (HAVE_FINITE) && defined (HAVE__FINITE)
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
56 #define finite _finite
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
57 #define HAVE_FINITE 1
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
58 #endif
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
59
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
60 #if ! defined (HAVE_COPYSIGN) && defined (HAVE__COPYSIGN)
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
61 #define copysign _copysign
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
62 #define HAVE_COPYSIGN 1
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
63 #endif
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
64
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
65 #if defined (_AIX) && defined (__GNUG__)
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
66 #undef finite
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
67 #define finite(x) ((x) < DBL_MAX && (x) > -DBL_MAX)
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
68 #endif
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
69
4102
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
70 /* Octave's idea of infinity. */
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
71 double octave_Inf;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
72 float octave_Float_Inf;
4102
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
73
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
74 /* Octave's idea of a missing value. */
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
75 double octave_NA;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
76 float octave_Float_NA;
4102
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
77
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
78 /* Octave's idea of not a number. */
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
79 double octave_NaN;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
80 float octave_Float_NaN;
4102
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
81
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
82 int lo_ieee_hw;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
83 int lo_ieee_lw;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
84
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
85 #if defined (SCO)
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
86
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
87 int
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
88 __isnan (double x)
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
89 {
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
90 return (IsNANorINF (x) && NaN (x) && ! IsINF (x)) ? 1 : 0;
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
91 }
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
92
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
93 int
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
94 __isinf (double x)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
95 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
96 return (IsNANorINF (x) && IsINF (x)) ? 1 : 0;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
97 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
98
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
99 int
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
100 __isnanf (float x)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
101 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
102 return (IsNANorINF (x) && NaN (x) && ! IsINF (x)) ? 1 : 0;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
103 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
104
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
105 int
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
106 __isinff (float x)
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
107 {
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
108 return (IsNANorINF (x) && IsINF (x)) ? 1 : 0;
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
109 }
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
110
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
111 #endif
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
112
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
113 int
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
114 __lo_ieee_isnan (double x)
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
115 {
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
116 #if defined (HAVE_ISNAN)
5389
25c8956d2204 [project @ 2005-06-15 03:45:46 by jwe]
jwe
parents: 5307
diff changeset
117 return isnan (x);
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
118 #else
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
119 return 0;
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
120 #endif
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
121 }
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
122
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
123 int
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
124 __lo_ieee_finite (double x)
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
125 {
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
126 #if defined (HAVE_FINITE)
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
127 return finite (x) != 0 && ! __lo_ieee_isnan (x);
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
128 #elif defined (HAVE_ISINF)
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
129 return (! isinf (x) && ! __lo_ieee_isnan (x));
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
130 #else
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
131 return ! __lo_ieee_isnan (x);
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
132 #endif
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
133 }
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
134
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
135 int
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
136 __lo_ieee_isinf (double x)
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
137 {
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
138 #if defined (HAVE_ISINF)
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
139 return isinf (x);
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
140 #elif defined (HAVE_FINITE)
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
141 return (! (finite (x) || __lo_ieee_isnan (x)));
4074
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
142 #else
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
143 return 0;
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe]
jwe
parents: 4072
diff changeset
144 #endif
4072
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
145 }
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe]
jwe
parents:
diff changeset
146
4075
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
147 int
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
148 __lo_ieee_is_NA (double x)
4075
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
149 {
6061
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe]
jwe
parents: 5389
diff changeset
150 #if defined (HAVE_ISNAN)
4075
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
151 lo_ieee_double t;
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
152 t.value = x;
7991
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
153 return (isnan (x) && t.word[lo_ieee_hw] == LO_IEEE_NA_HW
10317
42d098307c30 untabify additional source files
John W. Eaton <jwe@octave.org>
parents: 10158
diff changeset
154 && t.word[lo_ieee_lw] == LO_IEEE_NA_LW) ? 1 : 0;
4469
53ee020af847 [project @ 2003-07-26 03:45:10 by jwe]
jwe
parents: 4102
diff changeset
155 #else
53ee020af847 [project @ 2003-07-26 03:45:10 by jwe]
jwe
parents: 4102
diff changeset
156 return 0;
53ee020af847 [project @ 2003-07-26 03:45:10 by jwe]
jwe
parents: 4102
diff changeset
157 #endif
4075
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
158 }
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
159
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
160 int
7991
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
161 __lo_ieee_is_old_NA (double x)
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
162 {
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
163 #if defined (HAVE_ISNAN)
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
164 lo_ieee_double t;
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
165 t.value = x;
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
166 return (isnan (x) && t.word[lo_ieee_lw] == LO_IEEE_NA_LW_OLD
10317
42d098307c30 untabify additional source files
John W. Eaton <jwe@octave.org>
parents: 10158
diff changeset
167 && t.word[lo_ieee_hw] == LO_IEEE_NA_HW_OLD) ? 1 : 0;
7991
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
168 #else
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
169 return 0;
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
170 #endif
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
171 }
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
172
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
173 double
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
174 __lo_ieee_replace_old_NA (double x)
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
175 {
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
176 if (__lo_ieee_is_old_NA (x))
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
177 return lo_ieee_na_value ();
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
178 else
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
179 return x;
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
180 }
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
181
139f47cf17ab Change NA value to support single to double precision conversion
David Bateman <dbateman@free.fr>
parents: 7814
diff changeset
182 int
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
183 __lo_ieee_is_NaN_or_NA (double x)
4075
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
184 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
185 return __lo_ieee_isnan (x);
4075
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
186 }
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe]
jwe
parents: 4074
diff changeset
187
4102
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
188 double
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
189 lo_ieee_inf_value (void)
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
190 {
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
191 return octave_Inf;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
192 }
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
193
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
194 double
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
195 lo_ieee_na_value (void)
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
196 {
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
197 return octave_NA;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
198 }
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
199
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
200 double
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
201 lo_ieee_nan_value (void)
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
202 {
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
203 return octave_NaN;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
204 }
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe]
jwe
parents: 4075
diff changeset
205
5098
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
206 #if ! (defined (signbit) || defined (HAVE_DECL_SIGNBIT)) && defined (HAVE_SIGNBIT)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
207 extern int signbit (double);
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
208 #endif
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
209
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
210 int
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
211 __lo_ieee_signbit (double x)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
212 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
213 /* In the following definitions, only check x < 0 explicitly to avoid
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
214 a function call when it looks like signbit or copysign are actually
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
215 functions. */
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
216
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
217 #if defined (signbit)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
218 return signbit (x);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
219 #elif defined (HAVE_SIGNBIT)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
220 return (x < 0 || signbit (x));
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
221 #elif defined (copysign)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
222 return (copysign (1.0, x) < 0);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
223 #elif defined (HAVE_COPYSIGN)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
224 return (x < 0 || copysign (1.0, x) < 0);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
225 #else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
226 return x < 0;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
227 #endif
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
228 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
229
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
230 int
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
231 __lo_ieee_float_isnan (float x)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
232 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
233 #if defined (HAVE_ISNAN)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
234 return isnan (x);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
235 #else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
236 return 0;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
237 #endif
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
238 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
239
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
240 int
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
241 __lo_ieee_float_finite (float x)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
242 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
243 #if defined (HAVE_FINITE)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
244 return finite (x) != 0 && ! __lo_ieee_float_isnan (x);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
245 #elif defined (HAVE_ISINF)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
246 return (! isinf (x) && ! __lo_ieee_float_isnan (x));
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
247 #else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
248 return ! __lo_ieee_float_isnan (x);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
249 #endif
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
250 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
251
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
252 int
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
253 __lo_ieee_float_isinf (float x)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
254 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
255 #if defined (HAVE_ISINF)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
256 return isinf (x);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
257 #elif defined (HAVE_FINITE)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
258 return (! (finite (x) || __lo_ieee_float_isnan (x)));
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
259 #else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
260 return 0;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
261 #endif
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
262 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
263
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
264 int
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
265 __lo_ieee_float_is_NA (float x)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
266 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
267 #if defined (HAVE_ISNAN)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
268 lo_ieee_float t;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
269 t.value = x;
7814
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
270 return (isnan (x) && (t.word == LO_IEEE_NA_FLOAT)) ? 1 : 0;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
271 #else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
272 return 0;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
273 #endif
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
274 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
275
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
276 int
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
277 __lo_ieee_float_is_NaN_or_NA (float x)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
278 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
279 return __lo_ieee_float_isnan (x);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
280 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
281
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
282 float
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
283 lo_ieee_float_inf_value (void)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
284 {
7814
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
285 return octave_Float_Inf;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
286 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
287
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
288 float
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
289 lo_ieee_float_na_value (void)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
290 {
7814
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
291 return octave_Float_NA;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
292 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
293
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
294 float
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
295 lo_ieee_float_nan_value (void)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
296 {
7814
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
297 return octave_Float_NaN;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
298 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
299
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
300 #if ! (defined (signbit) || defined (HAVE_DECL_SIGNBIT)) && defined (HAVE_SIGNBIT)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
301 extern int signbit (float);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
302 #endif
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
303
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
304 int
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7231
diff changeset
305 __lo_ieee_float_signbit (float x)
5098
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
306 {
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
307 /* In the following definitions, only check x < 0 explicitly to avoid
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
308 a function call when it looks like signbit or copysign are actually
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
309 functions. */
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
310
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
311 #if defined (signbit)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
312 return signbit (x);
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
313 #elif defined (HAVE_SIGNBIT)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
314 return (x < 0 || signbit (x));
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
315 #elif defined (copysign)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
316 return (copysign (1.0, x) < 0);
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
317 #elif defined (HAVE_COPYSIGN)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
318 return (x < 0 || copysign (1.0, x) < 0);
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
319 #else
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
320 return x < 0;
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
321 #endif
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe]
jwe
parents: 4498
diff changeset
322 }