1 /* $OpenBSD: isnan.c,v 1.9 2013/11/12 06:09:48 deraadt Exp $ */ 2 /* 3 * Copyright (c) 2008 Martynas Venckus <martynas@openbsd.org> 4 * 5 * Permission to use, copy, modify, and distribute this software for any 6 * purpose with or without fee is hereby granted, provided that the above 7 * copyright notice and this permission notice appear in all copies. 8 * 9 * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES 10 * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF 11 * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR 12 * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES 13 * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN 14 * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF 15 * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. 16 */ 17 18 #include <sys/types.h> 19 #include <machine/ieee.h> 20 #include <float.h> 21 #include <math.h> 22 23 int 24 __isnan(double d) 25 { 26 struct ieee_double *p = (struct ieee_double *)&d; 27 28 return (p->dbl_exp == DBL_EXP_INFNAN && 29 (p->dbl_frach != 0 || p->dbl_fracl != 0)); 30 } 31 32 int 33 __isnanf(float f) 34 { 35 struct ieee_single *p = (struct ieee_single *)&f; 36 37 return (p->sng_exp == SNG_EXP_INFNAN && p->sng_frac != 0); 38 } 39 40 #if LDBL_MANT_DIG == DBL_MANT_DIG 41 __strong_alias(__isnanl, __isnan); 42 #endif /* LDBL_MANT_DIG == DBL_MANT_DIG */ 43 44 /* 45 * 3BSD compatibility aliases. 46 */ 47 __strong_alias(isnan, __isnan); 48 __strong_alias(isnanf, __isnanf); 49