/* * CDDL HEADER START * * The contents of this file are subject to the terms of the * Common Development and Distribution License, Version 1.0 only * (the "License"). You may not use this file except in compliance * with the License. * * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE * or http://www.opensolaris.org/os/licensing. * See the License for the specific language governing permissions * and limitations under the License. * * When distributing Covered Code, include this CDDL HEADER in each * file and include the License file at usr/src/OPENSOLARIS.LICENSE. * If applicable, add the following below this CDDL HEADER, with the * fields enclosed by brackets "[]" replaced with your own identifying * information: Portions Copyright [yyyy] [name of copyright owner] * * CDDL HEADER END */ /* Copyright (c) 1988 AT&T */ /* All Rights Reserved */ /* * Copyright 2014 Garrett D'Amore <garrett@damore.org> * * Copyright (c) 1996, by Sun Microsystems, Inc. * All Rights Reserved */ #ifndef _NAN_H #define _NAN_H /* * Handling of Not_a_Number's (only in IEEE floating-point standard) */ #include <sys/isa_defs.h> #include <values.h> #ifdef __cplusplus extern "C" { #endif #if defined(_IEEE_754) /* * Structure order is endian dependent. Only the common variants of * big and little endian are supported. */ #if defined(_BIG_ENDIAN) typedef union { struct { unsigned sign : 1; unsigned exponent :11; unsigned bits :20; unsigned fraction_low :32; } inf_parts; struct { unsigned sign : 1; unsigned exponent :11; unsigned qnan_bit : 1; unsigned bits :19; unsigned fraction_low :32; } nan_parts; double d; } dnan; #else /* Must be _LITTLE_ENDIAN */ typedef union { struct { unsigned fraction_low :32; unsigned bits :20; unsigned exponent :11; unsigned sign : 1; } inf_parts; struct { unsigned fraction_low :32; unsigned bits :19; unsigned qnan_bit : 1; unsigned exponent :11; unsigned sign : 1; } nan_parts; double d; } dnan; #endif /* Endian based selection */ /* * IsNANorINF checks that exponent of double == 2047 * i.e. that number is a NaN or an infinity */ #define IsNANorINF(X) (((dnan *)&(X))->nan_parts.exponent == 0x7ff) /* * IsINF must be used after IsNANorINF has checked the exponent */ #define IsINF(X) (((dnan *)&(X))->inf_parts.bits == 0 && \ ((dnan *)&(X))->inf_parts.fraction_low == 0) /* * IsPosNAN and IsNegNAN can be used to check the sign of infinities too */ #define IsPosNAN(X) (((dnan *)&(X))->nan_parts.sign == 0) #define IsNegNAN(X) (((dnan *)&(X))->nan_parts.sign == 1) /* * GETNaNPC gets the leftmost 32 bits of the fraction part */ #define GETNaNPC(dval) (((dnan *)&(dval))->inf_parts.bits << 12 | \ ((dnan *)&(dval))->nan_parts.fraction_low >> 20) #define KILLFPE() (void) _kill(_getpid(), 8) #define NaN(X) (((dnan *)&(X))->nan_parts.exponent == 0x7ff) #define KILLNaN(X) if (NaN(X)) KILLFPE() #else /* defined(_IEEE_754) */ /* #error is strictly ansi-C, but works as well as anything for K&R systems. */ #error ISA not supported #endif /* defined(_IEEE_754) */ #ifdef __cplusplus } #endif #endif /* _NAN_H */