octave-nkf: liboctave/lo-cieee.c annotate

annotate liboctave/lo-cieee.c @ 11542:695141f1c05c ss-3-3-55

snapshot 3.3.55

author	John W. Eaton <jwe@octave.org>
date	Sat, 15 Jan 2011 04:53:04 -0500
parents	fd0a3ac60b0e
children	12df7854fa7c

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
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 11235 diff changeset	3 Copyright (C) 2002-2011 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 #include "lo-ieee.h"
7231 2eb392d058bb [project @ 2007-11-30 18:53:29 by jwe] jwe parents: 7017 diff changeset	38 #include "lo-math.h"
4072 3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe] jwe parents: diff changeset	39
6061 c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	40 #if ! defined (HAVE_ISNAN) && defined (HAVE__ISNAN)
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	41 #define isnan _isnan
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	42 #define HAVE_ISNAN 1
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	43 #endif
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	44
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	45 #if ! defined (HAVE_FINITE) && defined (HAVE__FINITE)
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	46 #define finite _finite
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	47 #define HAVE_FINITE 1
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	48 #endif
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	49
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	50 #if ! defined (HAVE_COPYSIGN) && defined (HAVE__COPYSIGN)
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	51 #define copysign _copysign
c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	52 #define HAVE_COPYSIGN 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
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	55 #if defined (_AIX) && defined (__GNUG__)
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	56 #undef finite
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	57 #define finite(x) ((x) < DBL_MAX && (x) > -DBL_MAX)
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	58 #endif
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	59
4102 cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	60 /* Octave's idea of infinity. */
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	61 double octave_Inf;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	62 float octave_Float_Inf;
4102 cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	63
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	64 /* Octave's idea of a missing value. */
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	65 double octave_NA;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	66 float octave_Float_NA;
4102 cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	67
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	68 /* Octave's idea of not a number. */
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	69 double octave_NaN;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	70 float octave_Float_NaN;
4102 cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	71
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	72 int lo_ieee_hw;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	73 int lo_ieee_lw;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	74
4072 3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe] jwe parents: diff changeset	75 int
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	76 __lo_ieee_isnan (double x)
4072 3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe] jwe parents: diff changeset	77 {
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	78 #if defined (HAVE_ISNAN)
5389 25c8956d2204 [project @ 2005-06-15 03:45:46 by jwe] jwe parents: 5307 diff changeset	79 return isnan (x);
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	80 #else
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	81 return 0;
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	82 #endif
4072 3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe] jwe parents: diff changeset	83 }
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe] jwe parents: diff changeset	84
3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe] jwe parents: diff changeset	85 int
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	86 __lo_ieee_finite (double x)
4072 3cc39e3b8fa5 [project @ 2002-09-27 17:54:53 by jwe] jwe parents: diff changeset	87 {
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	88 #if defined (HAVE_FINITE)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	89 return finite (x) != 0 && ! __lo_ieee_isnan (x);
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	90 #elif defined (HAVE_ISINF)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	91 return (! isinf (x) && ! __lo_ieee_isnan (x));
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	92 #else
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	93 return ! __lo_ieee_isnan (x);
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	94 #endif
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	95 }
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	96
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	97 int
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	98 __lo_ieee_isinf (double x)
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	99 {
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	100 #if defined (HAVE_ISINF)
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	101 return isinf (x);
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	102 #elif defined (HAVE_FINITE)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	103 return (! (finite (x) \|\| __lo_ieee_isnan (x)));
4074 a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	104 #else
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	105 return 0;
a67f76924703 [project @ 2002-09-27 21:06:02 by jwe] jwe parents: 4072 diff changeset	106 #endif
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
4075 9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	109 int
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	110 __lo_ieee_is_NA (double x)
4075 9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	111 {
6061 c968f4198067 [project @ 2006-10-18 20:57:04 by jwe] jwe parents: 5389 diff changeset	112 #if defined (HAVE_ISNAN)
4075 9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	113 lo_ieee_double t;
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	114 t.value = x;
7991 139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	115 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	116 && 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	117 #else
53ee020af847 [project @ 2003-07-26 03:45:10 by jwe] jwe parents: 4102 diff changeset	118 return 0;
53ee020af847 [project @ 2003-07-26 03:45:10 by jwe] jwe parents: 4102 diff changeset	119 #endif
4075 9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	120 }
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	121
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	122 int
7991 139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	123 __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	124 {
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	125 #if defined (HAVE_ISNAN)
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	126 lo_ieee_double t;
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	127 t.value = x;
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	128 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	129 && 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	130 #else
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	131 return 0;
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	132 #endif
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	133 }
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	134
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	135 double
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	136 __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	137 {
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	138 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	139 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	140 else
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	141 return x;
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	142 }
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	143
139f47cf17ab Change NA value to support single to double precision conversion David Bateman <dbateman@free.fr> parents: 7814 diff changeset	144 int
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	145 __lo_ieee_is_NaN_or_NA (double x)
4075 9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	146 {
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	147 return __lo_ieee_isnan (x);
4075 9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	148 }
9a2102864eb0 [project @ 2002-09-27 21:37:56 by jwe] jwe parents: 4074 diff changeset	149
4102 cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	150 double
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	151 lo_ieee_inf_value (void)
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	152 {
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	153 return octave_Inf;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	154 }
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	155
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	156 double
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	157 lo_ieee_na_value (void)
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	158 {
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	159 return octave_NA;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	160 }
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	161
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	162 double
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	163 lo_ieee_nan_value (void)
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	164 {
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	165 return octave_NaN;
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	166 }
cbac6756967e [project @ 2002-10-14 20:16:44 by jwe] jwe parents: 4075 diff changeset	167
5098 ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	168 #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	169 extern int signbit (double);
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	170 #endif
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	171
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	172 int
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	173 __lo_ieee_signbit (double x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	174 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	175 /* 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	176 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	177 functions. */
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	178
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	179 #if defined (signbit)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	180 return signbit (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	181 #elif defined (HAVE_SIGNBIT)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	182 return (x < 0 \|\| signbit (x));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	183 #elif defined (copysign)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	184 return (copysign (1.0, x) < 0);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	185 #elif defined (HAVE_COPYSIGN)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	186 return (x < 0 \|\| copysign (1.0, x) < 0);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	187 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	188 return x < 0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	189 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	190 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	191
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	192 int
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	193 __lo_ieee_float_isnan (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	194 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	195 #if defined (HAVE_ISNAN)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	196 return isnan (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	197 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	198 return 0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	199 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	200 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	201
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	202 int
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	203 __lo_ieee_float_finite (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	204 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	205 #if defined (HAVE_FINITE)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	206 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	207 #elif defined (HAVE_ISINF)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	208 return (! isinf (x) && ! __lo_ieee_float_isnan (x));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	209 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	210 return ! __lo_ieee_float_isnan (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	211 #endif
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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	214 int
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	215 __lo_ieee_float_isinf (float x)
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 (HAVE_ISINF)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	218 return isinf (x);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	219 #elif defined (HAVE_FINITE)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	220 return (! (finite (x) \|\| __lo_ieee_float_isnan (x)));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	221 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	222 return 0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	223 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	224 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	225
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	226 int
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	227 __lo_ieee_float_is_NA (float x)
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 #if defined (HAVE_ISNAN)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	230 lo_ieee_float t;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	231 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	232 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	233 #else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	234 return 0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	235 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	236 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	237
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	238 int
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	239 __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	240 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	241 return __lo_ieee_float_isnan (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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	244 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	245 lo_ieee_float_inf_value (void)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	246 {
7814 87865ed7405f Second set of single precision test code and fix of resulting bugs David Bateman <dbateman@free.fr> parents: 7789 diff changeset	247 return octave_Float_Inf;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	248 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	249
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	250 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	251 lo_ieee_float_na_value (void)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	252 {
7814 87865ed7405f Second set of single precision test code and fix of resulting bugs David Bateman <dbateman@free.fr> parents: 7789 diff changeset	253 return octave_Float_NA;
7789 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
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	256 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	257 lo_ieee_float_nan_value (void)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	258 {
7814 87865ed7405f Second set of single precision test code and fix of resulting bugs David Bateman <dbateman@free.fr> parents: 7789 diff changeset	259 return octave_Float_NaN;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	260 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	261
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	262 #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	263 extern int signbit (float);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	264 #endif
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	265
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	266 int
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7231 diff changeset	267 __lo_ieee_float_signbit (float x)
5098 ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	268 {
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	269 /* 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	270 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	271 functions. */
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	272
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	273 #if defined (signbit)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	274 return signbit (x);
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	275 #elif defined (HAVE_SIGNBIT)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	276 return (x < 0 \|\| signbit (x));
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	277 #elif defined (copysign)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	278 return (copysign (1.0, x) < 0);
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	279 #elif defined (HAVE_COPYSIGN)
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	280 return (x < 0 \|\| copysign (1.0, x) < 0);
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	281 #else
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	282 return x < 0;
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	283 #endif
ab4e64f92526 [project @ 2004-12-18 15:04:20 by jwe] jwe parents: 4498 diff changeset	284 }

Mercurial > octave-nkf

annotate liboctave/lo-cieee.c @ 11542:695141f1c05c ss-3-3-55