1a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/************************************************************************* 2a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 3df143a5041f03a22808b59c76698770b74692815Selim Gurun * $Id$ 4a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 5a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Copyright (C) 2001 Bjorn Reese <breese@users.sourceforge.net> 6a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 7a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Permission to use, copy, modify, and distribute this software for any 8a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * purpose with or without fee is hereby granted, provided that the above 9a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * copyright notice and this permission notice appear in all copies. 10a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 11a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * THIS SOFTWARE IS PROVIDED ``AS IS'' AND WITHOUT ANY EXPRESS OR IMPLIED 12a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * WARRANTIES, INCLUDING, WITHOUT LIMITATION, THE IMPLIED WARRANTIES OF 13a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * MERCHANTIBILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE AUTHORS AND 14a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * CONTRIBUTORS ACCEPT NO RESPONSIBILITY IN ANY CONCEIVABLE MANNER. 15a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 16a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ************************************************************************ 17a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 18a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Functions to handle special quantities in floating-point numbers 19a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * (that is, NaNs and infinity). They provide the capability to detect 20a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * and fabricate special quantities. 21a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 22a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Although written to be as portable as possible, it can never be 23a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * guaranteed to work on all platforms, as not all hardware supports 24a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * special quantities. 25a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 26a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * The approach used here (approximately) is to: 27a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 28a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 1. Use C99 functionality when available. 29a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 2. Use IEEE 754 bit-patterns if possible. 30a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 3. Use platform-specific techniques. 31a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 32a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ************************************************************************/ 33a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 34a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* 35a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * TODO: 36a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * o Put all the magic into trio_fpclassify_and_signbit(), and use this from 37a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * trio_isnan() etc. 38a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 39a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 40a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/************************************************************************* 41a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Include files 42a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 43a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#include "triodef.h" 44a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#include "trionan.h" 45a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 46a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#include <math.h> 47a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#include <string.h> 48a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#include <limits.h> 49a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#include <float.h> 50a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(TRIO_PLATFORM_UNIX) 51a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# include <signal.h> 52a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 53a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(TRIO_COMPILER_DECC) 54a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(__linux__) 55a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# include <cpml.h> 56a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# else 57a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# include <fp_class.h> 58a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 59a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 60a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#include <assert.h> 61a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 62a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(TRIO_DOCUMENTATION) 63a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# include "doc/doc_nan.h" 64a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 65a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** @addtogroup SpecialQuantities 66a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @{ 67a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 68a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 69a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/************************************************************************* 70a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Definitions 71a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 72a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 73a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#define TRIO_TRUE (1 == 1) 74a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#define TRIO_FALSE (0 == 1) 75a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 76a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* 77a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * We must enable IEEE floating-point on Alpha 78a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 79a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(__alpha) && !defined(_IEEE_FP) 80a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_COMPILER_DECC) 81a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_VMS) 82a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# error "Must be compiled with option /IEEE_MODE=UNDERFLOW_TO_ZERO/FLOAT=IEEE" 83a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# else 84a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if !defined(_CFE) 85a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# error "Must be compiled with option -ieee" 86a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 87a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 88a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# elif defined(TRIO_COMPILER_GCC) && (defined(__osf__) || defined(__linux__)) 89a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# error "Must be compiled with option -mieee" 90a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 91a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif /* __alpha && ! _IEEE_FP */ 92a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 93a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* 94a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * In ANSI/IEEE 754-1985 64-bits double format numbers have the 95a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * following properties (amoungst others) 96a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 97a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * o FLT_RADIX == 2: binary encoding 98a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * o DBL_MAX_EXP == 1024: 11 bits exponent, where one bit is used 99a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * to indicate special numbers (e.g. NaN and Infinity), so the 100a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * maximum exponent is 10 bits wide (2^10 == 1024). 101a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * o DBL_MANT_DIG == 53: The mantissa is 52 bits wide, but because 102a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * numbers are normalized the initial binary 1 is represented 103a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * implicitly (the so-called "hidden bit"), which leaves us with 104a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * the ability to represent 53 bits wide mantissa. 105a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 106a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if (FLT_RADIX == 2) && (DBL_MAX_EXP == 1024) && (DBL_MANT_DIG == 53) 107a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define USE_IEEE_754 108a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 109a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 110a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 111a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/************************************************************************* 112a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Constants 113a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 114a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 115df143a5041f03a22808b59c76698770b74692815Selim Gurunstatic TRIO_CONST char rcsid[] = "@(#)$Id$"; 116a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 117a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(USE_IEEE_754) 118a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 119a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* 120a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Endian-agnostic indexing macro. 121a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 122a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * The value of internalEndianMagic, when converted into a 64-bit 123a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * integer, becomes 0x0706050403020100 (we could have used a 64-bit 124a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * integer value instead of a double, but not all platforms supports 125a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * that type). The value is automatically encoded with the correct 126a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * endianess by the compiler, which means that we can support any 127a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * kind of endianess. The individual bytes are then used as an index 128a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * for the IEEE 754 bit-patterns and masks. 129a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 130a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#define TRIO_DOUBLE_INDEX(x) (((unsigned char *)&internalEndianMagic)[7-(x)]) 131a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 132a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if (defined(__BORLANDC__) && __BORLANDC__ >= 0x0590) 133a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST double internalEndianMagic = 7.949928895127362e-275; 134a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#else 135a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST double internalEndianMagic = 7.949928895127363e-275; 136a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 137a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 138a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* Mask for the exponent */ 139a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST unsigned char ieee_754_exponent_mask[] = { 140a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 0x7F, 0xF0, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00 141a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang}; 142a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 143a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* Mask for the mantissa */ 144a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST unsigned char ieee_754_mantissa_mask[] = { 145a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 0x00, 0x0F, 0xFF, 0xFF, 0xFF, 0xFF, 0xFF, 0xFF 146a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang}; 147a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 148a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* Mask for the sign bit */ 149a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST unsigned char ieee_754_sign_mask[] = { 150a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 0x80, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00 151a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang}; 152a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 153a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* Bit-pattern for negative zero */ 154a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST unsigned char ieee_754_negzero_array[] = { 155a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 0x80, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00 156a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang}; 157a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 158a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* Bit-pattern for infinity */ 159a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST unsigned char ieee_754_infinity_array[] = { 160a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 0x7F, 0xF0, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00 161a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang}; 162a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 163a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* Bit-pattern for quiet NaN */ 164a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST unsigned char ieee_754_qnan_array[] = { 165a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 0x7F, 0xF8, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00 166a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang}; 167a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 168a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 169a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/************************************************************************* 170a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Functions 171a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 172a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 173a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* 174a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * trio_make_double 175a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 176a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PRIVATE double 177a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_make_double 178a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS1((values), 179a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang TRIO_CONST unsigned char *values) 180a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 181a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang TRIO_VOLATILE double result; 182a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int i; 183a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 184a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang for (i = 0; i < (int)sizeof(double); i++) { 185a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((TRIO_VOLATILE unsigned char *)&result)[TRIO_DOUBLE_INDEX(i)] = values[i]; 186a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 187a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return result; 188a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 189a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 190a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* 191a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * trio_is_special_quantity 192a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 193a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PRIVATE int 194a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_is_special_quantity 195a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS2((number, has_mantissa), 196a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number, 197a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int *has_mantissa) 198a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 199a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang unsigned int i; 200a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang unsigned char current; 201a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int is_special_quantity = TRIO_TRUE; 202a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 203a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *has_mantissa = 0; 204a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 205a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang for (i = 0; i < (unsigned int)sizeof(double); i++) { 206a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang current = ((unsigned char *)&number)[TRIO_DOUBLE_INDEX(i)]; 207a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang is_special_quantity 208a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang &= ((current & ieee_754_exponent_mask[i]) == ieee_754_exponent_mask[i]); 209a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *has_mantissa |= (current & ieee_754_mantissa_mask[i]); 210a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 211a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return is_special_quantity; 212a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 213a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 214a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* 215a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * trio_is_negative 216a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 217a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PRIVATE int 218a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_is_negative 219a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS1((number), 220a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number) 221a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 222a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang unsigned int i; 223a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int is_negative = TRIO_FALSE; 224a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 225a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang for (i = 0; i < (unsigned int)sizeof(double); i++) { 226a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang is_negative |= (((unsigned char *)&number)[TRIO_DOUBLE_INDEX(i)] 227a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang & ieee_754_sign_mask[i]); 228a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 229a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return is_negative; 230a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 231a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 232a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif /* USE_IEEE_754 */ 233a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 234a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 235a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 236a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Generate negative zero. 237a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 238a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return Floating-point representation of negative zero. 239a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 240a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC double 241a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_nzero(TRIO_NOARGS) 242a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 243a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(USE_IEEE_754) 244a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return trio_make_double(ieee_754_negzero_array); 245a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#else 246a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang TRIO_VOLATILE double zero = 0.0; 247a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 248a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return -zero; 249a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 250a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 251a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 252a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 253a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Generate positive infinity. 254a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 255a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return Floating-point representation of positive infinity. 256a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 257a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC double 258a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_pinf(TRIO_NOARGS) 259a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 260a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* Cache the result */ 261a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang static double result = 0.0; 262a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 263a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if (result == 0.0) { 26494442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 265a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(INFINITY) && defined(__STDC_IEC_559__) 266a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result = (double)INFINITY; 267a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 268a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(USE_IEEE_754) 269a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result = trio_make_double(ieee_754_infinity_array); 270a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 271a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#else 272a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 273a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * If HUGE_VAL is different from DBL_MAX, then HUGE_VAL is used 274a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * as infinity. Otherwise we have to resort to an overflow 275a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * operation to generate infinity. 276a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 277a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 278a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang void (*signal_handler)(int) = signal(SIGFPE, SIG_IGN); 279a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 280a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 281a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result = HUGE_VAL; 282a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if (HUGE_VAL == DBL_MAX) { 283a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* Force overflow */ 284a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result += HUGE_VAL; 285a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 28694442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 287a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 288a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang signal(SIGFPE, signal_handler); 289a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 290a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 291a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 292a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 293a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return result; 294a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 295a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 296a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 297a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Generate negative infinity. 298a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 299a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return Floating-point value of negative infinity. 300a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 301a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC double 302a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_ninf(TRIO_NOARGS) 303a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 304a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang static double result = 0.0; 305a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 306a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if (result == 0.0) { 307a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 308a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Negative infinity is calculated by negating positive infinity, 309a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * which can be done because it is legal to do calculations on 310a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * infinity (for example, 1 / infinity == 0). 311a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 312a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result = -trio_pinf(); 313a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 314a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return result; 315a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 316a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 317a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 318a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Generate NaN. 319a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 320a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return Floating-point representation of NaN. 321a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 322a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC double 323a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_nan(TRIO_NOARGS) 324a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 325a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* Cache the result */ 326a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang static double result = 0.0; 327a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 328a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if (result == 0.0) { 32994442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 330a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(TRIO_COMPILER_SUPPORTS_C99) 331a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result = nan(""); 332a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 333a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(NAN) && defined(__STDC_IEC_559__) 334a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result = (double)NAN; 33594442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 336a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(USE_IEEE_754) 337a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result = trio_make_double(ieee_754_qnan_array); 338a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 339a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#else 340a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 341a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * There are several ways to generate NaN. The one used here is 342a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * to divide infinity by infinity. I would have preferred to add 343a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * negative infinity to positive infinity, but that yields wrong 344a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * result (infinity) on FreeBSD. 345a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 346a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * This may fail if the hardware does not support NaN, or if 347a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * the Invalid Operation floating-point exception is unmasked. 348a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 349a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 350a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang void (*signal_handler)(int) = signal(SIGFPE, SIG_IGN); 351a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 35294442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 353a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang result = trio_pinf() / trio_pinf(); 35494442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 355a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 356a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang signal(SIGFPE, signal_handler); 357a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 35894442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 359a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 360a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 361a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return result; 362a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 363a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 364a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 365a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Check for NaN. 366a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 367a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @param number An arbitrary floating-point number. 368a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return Boolean value indicating whether or not the number is a NaN. 369a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 370a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC int 371a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_isnan 372a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS1((number), 373a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number) 374a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 375a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if (defined(TRIO_COMPILER_SUPPORTS_C99) && defined(isnan)) \ 376a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang || defined(TRIO_COMPILER_SUPPORTS_UNIX95) 377a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 378a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * C99 defines isnan() as a macro. UNIX95 defines isnan() as a 379a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * function. This function was already present in XPG4, but this 380a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * is a bit tricky to detect with compiler defines, so we choose 381a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * the conservative approach and only use it for UNIX95. 382a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 383a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return isnan(number); 38494442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 385a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(TRIO_COMPILER_MSVC) || defined(TRIO_COMPILER_BCB) 386a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 387a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Microsoft Visual C++ and Borland C++ Builder have an _isnan() 388a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * function. 389a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 390a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return _isnan(number) ? TRIO_TRUE : TRIO_FALSE; 391a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 392a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(USE_IEEE_754) 393a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 394a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Examine IEEE 754 bit-pattern. A NaN must have a special exponent 395a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * pattern, and a non-empty mantissa. 396a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 397a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int has_mantissa; 398a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int is_special_quantity; 399a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 400a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang is_special_quantity = trio_is_special_quantity(number, &has_mantissa); 40194442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 402a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return (is_special_quantity && has_mantissa); 40394442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 404a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#else 405a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 406a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Fallback solution 407a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 408a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int status; 409a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double integral, fraction; 41094442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 411a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 412a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang void (*signal_handler)(int) = signal(SIGFPE, SIG_IGN); 413a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 41494442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 415a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang status = (/* 416a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * NaN is the only number which does not compare to itself 417a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 418a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((TRIO_VOLATILE double)number != (TRIO_VOLATILE double)number) || 419a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 420a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Fallback solution if NaN compares to NaN 421a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 422a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((number != 0.0) && 423a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang (fraction = modf(number, &integral), 424a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang integral == fraction))); 42594442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 426a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 427a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang signal(SIGFPE, signal_handler); 428a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 42994442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 430a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return status; 43194442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 432a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 433a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 434a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 435a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 436a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Check for infinity. 437a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 438a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @param number An arbitrary floating-point number. 439a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return 1 if positive infinity, -1 if negative infinity, 0 otherwise. 440a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 441a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC int 442a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_isinf 443a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS1((number), 444a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number) 445a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 446a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(TRIO_COMPILER_DECC) && !defined(__linux__) 447a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 448a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * DECC has an isinf() macro, but it works differently than that 449a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * of C99, so we use the fp_class() function instead. 450a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 451a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return ((fp_class(number) == FP_POS_INF) 452a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ? 1 453a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang : ((fp_class(number) == FP_NEG_INF) ? -1 : 0)); 454a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 455a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(isinf) 456a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 457a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * C99 defines isinf() as a macro. 458a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 459a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return isinf(number) 460a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ? ((number > 0.0) ? 1 : -1) 461a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang : 0; 46294442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 463a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(TRIO_COMPILER_MSVC) || defined(TRIO_COMPILER_BCB) 464a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 465a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Microsoft Visual C++ and Borland C++ Builder have an _fpclass() 466a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * function that can be used to detect infinity. 467a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 468a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return ((_fpclass(number) == _FPCLASS_PINF) 469a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ? 1 470a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang : ((_fpclass(number) == _FPCLASS_NINF) ? -1 : 0)); 471a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 472a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(USE_IEEE_754) 473a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 474a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Examine IEEE 754 bit-pattern. Infinity must have a special exponent 475a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * pattern, and an empty mantissa. 476a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 477a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int has_mantissa; 478a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int is_special_quantity; 479a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 480a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang is_special_quantity = trio_is_special_quantity(number, &has_mantissa); 48194442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 482a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return (is_special_quantity && !has_mantissa) 483a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ? ((number < 0.0) ? -1 : 1) 484a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang : 0; 485a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 486a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#else 487a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 488a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Fallback solution. 489a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 490a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int status; 49194442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 492a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 493a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang void (*signal_handler)(int) = signal(SIGFPE, SIG_IGN); 494a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 49594442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 496a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double infinity = trio_pinf(); 49794442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 498a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang status = ((number == infinity) 499a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ? 1 500a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang : ((number == -infinity) ? -1 : 0)); 50194442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 502a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 503a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang signal(SIGFPE, signal_handler); 504a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 50594442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 506a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return status; 50794442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 508a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 509a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 510a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 511a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if 0 512a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* Temporary fix - this routine is not used anywhere */ 513a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 514a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Check for finity. 515a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 516a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @param number An arbitrary floating-point number. 517a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return Boolean value indicating whether or not the number is a finite. 518a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 519a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC int 520a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_isfinite 521a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS1((number), 522a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number) 523a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 524a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(TRIO_COMPILER_SUPPORTS_C99) && defined(isfinite) 525a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 526a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * C99 defines isfinite() as a macro. 527a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 528a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return isfinite(number); 52994442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 530a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(TRIO_COMPILER_MSVC) || defined(TRIO_COMPILER_BCB) 531a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 532a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Microsoft Visual C++ and Borland C++ Builder use _finite(). 533a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 534a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return _finite(number); 535a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 536a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#elif defined(USE_IEEE_754) 537a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 538a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Examine IEEE 754 bit-pattern. For finity we do not care about the 539a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * mantissa. 540a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 541a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int dummy; 542a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 543a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return (! trio_is_special_quantity(number, &dummy)); 544a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 545a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#else 546a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 547a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Fallback solution. 548a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 549a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return ((trio_isinf(number) == 0) && (trio_isnan(number) == 0)); 55094442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 551a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 552a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 553a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 554a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 555a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 556a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/* 557a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * The sign of NaN is always false 558a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 559a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC int 560a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_fpclassify_and_signbit 561a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS2((number, is_negative), 562a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number, 563a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int *is_negative) 564a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 565a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(fpclassify) && defined(signbit) 566a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 567a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * C99 defines fpclassify() and signbit() as a macros 568a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 569a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = signbit(number); 570a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang switch (fpclassify(number)) { 571a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case FP_NAN: 572a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_NAN; 573a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case FP_INFINITE: 574a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_INFINITE; 575a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case FP_SUBNORMAL: 576a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_SUBNORMAL; 577a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case FP_ZERO: 578a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_ZERO; 579a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang default: 580a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_NORMAL; 581a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 582a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 583a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#else 584a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_COMPILER_DECC) 585a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 586a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * DECC has an fp_class() function. 587a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 588a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_FPCLASSIFY(n) fp_class(n) 589a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_QUIET_NAN FP_QNAN 590a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_SIGNALLING_NAN FP_SNAN 591a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_INFINITY FP_POS_INF 592a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_INFINITY FP_NEG_INF 593a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_SUBNORMAL FP_POS_DENORM 594a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_SUBNORMAL FP_NEG_DENORM 595a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_ZERO FP_POS_ZERO 596a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_ZERO FP_NEG_ZERO 597a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_NORMAL FP_POS_NORM 598a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_NORMAL FP_NEG_NORM 59994442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 600a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# elif defined(TRIO_COMPILER_MSVC) || defined(TRIO_COMPILER_BCB) 601a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 602a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Microsoft Visual C++ and Borland C++ Builder have an _fpclass() 603a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * function. 604a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 605a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_FPCLASSIFY(n) _fpclass(n) 606a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_QUIET_NAN _FPCLASS_QNAN 607a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_SIGNALLING_NAN _FPCLASS_SNAN 608a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_INFINITY _FPCLASS_PINF 609a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_INFINITY _FPCLASS_NINF 610a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_SUBNORMAL _FPCLASS_PD 611a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_SUBNORMAL _FPCLASS_ND 612a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_ZERO _FPCLASS_PZ 613a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_ZERO _FPCLASS_NZ 614a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_NORMAL _FPCLASS_PN 615a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_NORMAL _FPCLASS_NN 61694442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 617a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# elif defined(FP_PLUS_NORM) 618a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 619a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * HP-UX 9.x and 10.x have an fpclassify() function, that is different 620a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * from the C99 fpclassify() macro supported on HP-UX 11.x. 621a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 622a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * AIX has class() for C, and _class() for C++, which returns the 623a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * same values as the HP-UX fpclassify() function. 624a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 625a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_AIX) 626a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(__cplusplus) 627a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_FPCLASSIFY(n) _class(n) 628a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# else 629a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_FPCLASSIFY(n) class(n) 630a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 631a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# else 632a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_FPCLASSIFY(n) fpclassify(n) 633a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 634a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_QUIET_NAN FP_QNAN 635a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_SIGNALLING_NAN FP_SNAN 636a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_INFINITY FP_PLUS_INF 637a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_INFINITY FP_MINUS_INF 638a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_SUBNORMAL FP_PLUS_DENORM 639a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_SUBNORMAL FP_MINUS_DENORM 640a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_ZERO FP_PLUS_ZERO 641a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_ZERO FP_MINUS_ZERO 642a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_POSITIVE_NORMAL FP_PLUS_NORM 643a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# define TRIO_NEGATIVE_NORMAL FP_MINUS_NORM 644a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 645a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 646a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_FPCLASSIFY) 647a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang switch (TRIO_FPCLASSIFY(number)) { 648a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_QUIET_NAN: 649a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_SIGNALLING_NAN: 650a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_FALSE; /* NaN has no sign */ 651a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_NAN; 652a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_POSITIVE_INFINITY: 653a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_FALSE; 654a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_INFINITE; 655a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_NEGATIVE_INFINITY: 656a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_TRUE; 657a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_INFINITE; 658a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_POSITIVE_SUBNORMAL: 659a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_FALSE; 660a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_SUBNORMAL; 661a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_NEGATIVE_SUBNORMAL: 662a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_TRUE; 663a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_SUBNORMAL; 664a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_POSITIVE_ZERO: 665a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_FALSE; 666a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_ZERO; 667a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_NEGATIVE_ZERO: 668a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_TRUE; 669a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_ZERO; 670a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_POSITIVE_NORMAL: 671a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_FALSE; 672a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_NORMAL; 673a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_NEGATIVE_NORMAL: 674a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_TRUE; 675a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_NORMAL; 676a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang default: 677a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* Just in case... */ 678a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = (number < 0.0); 679a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_NORMAL; 680a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 68194442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 682a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# else 683a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 684a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Fallback solution. 685a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 686a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int rc; 68794442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 688a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if (number == 0.0) { 689a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* 690a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * In IEEE 754 the sign of zero is ignored in comparisons, so we 691a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * have to handle this as a special case by examining the sign bit 692a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * directly. 693a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 694a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(USE_IEEE_754) 695a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = trio_is_negative(number); 696a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# else 697a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_FALSE; /* FIXME */ 698a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 699a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_ZERO; 700a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 701a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if (trio_isnan(number)) { 702a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_FALSE; 703a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_NAN; 704a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 705a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if ((rc = trio_isinf(number))) { 706a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = (rc == -1); 707a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_INFINITE; 708a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 709a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if ((number > 0.0) && (number < DBL_MIN)) { 710a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_FALSE; 711a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_SUBNORMAL; 712a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 713a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang if ((number < 0.0) && (number > -DBL_MIN)) { 714a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = TRIO_TRUE; 715a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_SUBNORMAL; 716a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 717a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang *is_negative = (number < 0.0); 718a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return TRIO_FP_NORMAL; 71994442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 720a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 721a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 722a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 723a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 724a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 725a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Examine the sign of a number. 726a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 727a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @param number An arbitrary floating-point number. 728a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return Boolean value indicating whether or not the number has the 729a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang sign bit set (i.e. is negative). 730a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 731a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC int 732a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_signbit 733a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS1((number), 734a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number) 735a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 736a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int is_negative; 73794442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 738a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang (void)trio_fpclassify_and_signbit(number, &is_negative); 739a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return is_negative; 740a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 741a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 742a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if 0 743a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang /* Temporary fix - this routine is not used in libxml */ 744a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** 745a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang Examine the class of a number. 746a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 747a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @param number An arbitrary floating-point number. 748a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang @return Enumerable value indicating the class of @p number 749a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang*/ 750a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_PUBLIC int 751a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangtrio_fpclassify 752a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS1((number), 753a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number) 754a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 755a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int dummy; 75694442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 757a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return trio_fpclassify_and_signbit(number, &dummy); 758a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 759a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 760a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 761a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 762a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/** @} SpecialQuantities */ 763a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 764a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang/************************************************************************* 765a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * For test purposes. 766a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 767a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Add the following compiler option to include this test code. 768a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * 769a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * Unix : -DSTANDALONE 770a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang * VMS : /DEFINE=(STANDALONE) 771a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang */ 772a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#if defined(STANDALONE) 773a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# include <stdio.h> 774a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 775a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic TRIO_CONST char * 776a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WanggetClassification 777a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS1((type), 778a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang int type) 779a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 780a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang switch (type) { 781a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_FP_INFINITE: 782a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return "FP_INFINITE"; 783a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_FP_NAN: 784a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return "FP_NAN"; 785a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_FP_NORMAL: 786a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return "FP_NORMAL"; 787a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_FP_SUBNORMAL: 788a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return "FP_SUBNORMAL"; 789a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang case TRIO_FP_ZERO: 790a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return "FP_ZERO"; 791a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang default: 792a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return "FP_UNKNOWN"; 793a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang } 794a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 795a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 796a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangstatic void 797a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangprint_class 798a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) WangTRIO_ARGS2((prefix, number), 799a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang TRIO_CONST char *prefix, 800a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double number) 801a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 802a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang printf("%-6s: %s %-15s %g\n", 803a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang prefix, 804a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang trio_signbit(number) ? "-" : "+", 805a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang getClassification(TRIO_FPCLASSIFY(number)), 806a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang number); 807a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 808a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 809a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wangint main(TRIO_NOARGS) 810a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang{ 811a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double my_nan; 812a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double my_pinf; 813a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang double my_ninf; 814a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 815a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang void (*signal_handler) TRIO_PROTO((int)); 816a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 817a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 818a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_nan = trio_nan(); 819a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_pinf = trio_pinf(); 820a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_ninf = trio_ninf(); 821a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 822a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("Nan", my_nan); 823a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("PInf", my_pinf); 824a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("NInf", my_ninf); 825a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("PZero", 0.0); 826a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("NZero", -0.0); 827a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("PNorm", 1.0); 828a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("NNorm", -1.0); 829a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("PSub", 1.01e-307 - 1.00e-307); 830a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang print_class("NSub", 1.00e-307 - 1.01e-307); 83194442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 832a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang printf("NaN : %4g 0x%02x%02x%02x%02x%02x%02x%02x%02x (%2d, %2d)\n", 833a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_nan, 834a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[0], 835a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[1], 836a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[2], 837a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[3], 838a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[4], 839a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[5], 840a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[6], 841a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[7], 842a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang trio_isnan(my_nan), trio_isinf(my_nan)); 843a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang printf("PInf: %4g 0x%02x%02x%02x%02x%02x%02x%02x%02x (%2d, %2d)\n", 844a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_pinf, 845a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[0], 846a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[1], 847a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[2], 848a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[3], 849a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[4], 850a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[5], 851a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[6], 852a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[7], 853a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang trio_isnan(my_pinf), trio_isinf(my_pinf)); 854a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang printf("NInf: %4g 0x%02x%02x%02x%02x%02x%02x%02x%02x (%2d, %2d)\n", 855a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_ninf, 856a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[0], 857a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[1], 858a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[2], 859a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[3], 860a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[4], 861a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[5], 862a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[6], 863a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[7], 864a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang trio_isnan(my_ninf), trio_isinf(my_ninf)); 86594442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 866a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 867a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang signal_handler = signal(SIGFPE, SIG_IGN); 868a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 86994442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 870a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_pinf = DBL_MAX + DBL_MAX; 871a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_ninf = -my_pinf; 872a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_nan = my_pinf / my_pinf; 873a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang 874a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# if defined(TRIO_PLATFORM_UNIX) 875a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang signal(SIGFPE, signal_handler); 876a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang# endif 87794442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 878a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang printf("NaN : %4g 0x%02x%02x%02x%02x%02x%02x%02x%02x (%2d, %2d)\n", 879a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_nan, 880a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[0], 881a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[1], 882a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[2], 883a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[3], 884a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[4], 885a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[5], 886a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[6], 887a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_nan)[7], 888a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang trio_isnan(my_nan), trio_isinf(my_nan)); 889a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang printf("PInf: %4g 0x%02x%02x%02x%02x%02x%02x%02x%02x (%2d, %2d)\n", 890a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_pinf, 891a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[0], 892a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[1], 893a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[2], 894a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[3], 895a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[4], 896a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[5], 897a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[6], 898a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_pinf)[7], 899a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang trio_isnan(my_pinf), trio_isinf(my_pinf)); 900a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang printf("NInf: %4g 0x%02x%02x%02x%02x%02x%02x%02x%02x (%2d, %2d)\n", 901a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang my_ninf, 902a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[0], 903a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[1], 904a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[2], 905a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[3], 906a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[4], 907a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[5], 908a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[6], 909a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang ((unsigned char *)&my_ninf)[7], 910a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang trio_isnan(my_ninf), trio_isinf(my_ninf)); 91194442ad4107000e6d49f9b85a46a591495a57632Selim Gurun 912a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang return 0; 913a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang} 914a66dd5f0b030e773d3af602682150eebdc9f2bddShimeng (Simon) Wang#endif 915