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

annotate liboctave/lo-cieee.c @ 16088:b29b10fbb744 stable release-3-6-4

Version 3.6.4 released. * configure.ac (AC_INIT): Version is now 3.6.4. (OCTAVE_RELEASE_DATE): Now 2013-02-21.

author	John W. Eaton <jwe@octave.org>
date	Thu, 21 Feb 2013 15:17:54 -0500
parents	72c96de7a403
children

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

Mercurial > octave-nkf

annotate liboctave/lo-cieee.c @ 16088:b29b10fbb744 stable release-3-6-4