1240dbabfSDavid Schultz /*- 2240dbabfSDavid Schultz * Copyright (c) 2004 David Schultz <das@FreeBSD.ORG> 3240dbabfSDavid Schultz * All rights reserved. 4240dbabfSDavid Schultz * 5240dbabfSDavid Schultz * Redistribution and use in source and binary forms, with or without 6240dbabfSDavid Schultz * modification, are permitted provided that the following conditions 7240dbabfSDavid Schultz * are met: 8240dbabfSDavid Schultz * 1. Redistributions of source code must retain the above copyright 9240dbabfSDavid Schultz * notice, this list of conditions and the following disclaimer. 10240dbabfSDavid Schultz * 2. Redistributions in binary form must reproduce the above copyright 11240dbabfSDavid Schultz * notice, this list of conditions and the following disclaimer in the 12240dbabfSDavid Schultz * documentation and/or other materials provided with the distribution. 13240dbabfSDavid Schultz * 14240dbabfSDavid Schultz * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND 15240dbabfSDavid Schultz * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 16240dbabfSDavid Schultz * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 17240dbabfSDavid Schultz * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE 18240dbabfSDavid Schultz * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 19240dbabfSDavid Schultz * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 20240dbabfSDavid Schultz * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 21240dbabfSDavid Schultz * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 22240dbabfSDavid Schultz * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 23240dbabfSDavid Schultz * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 24240dbabfSDavid Schultz * SUCH DAMAGE. 25240dbabfSDavid Schultz * 26240dbabfSDavid Schultz * $FreeBSD$ 27240dbabfSDavid Schultz */ 28240dbabfSDavid Schultz 29240dbabfSDavid Schultz #include <math.h> 30240dbabfSDavid Schultz 31240dbabfSDavid Schultz #include "fpmath.h" 32240dbabfSDavid Schultz 33240dbabfSDavid Schultz int 34240dbabfSDavid Schultz __isfinite(double d) 35240dbabfSDavid Schultz { 36240dbabfSDavid Schultz union IEEEd2bits u; 37240dbabfSDavid Schultz 38240dbabfSDavid Schultz u.d = d; 39240dbabfSDavid Schultz return (u.bits.exp != 2047); 40240dbabfSDavid Schultz } 41240dbabfSDavid Schultz 42240dbabfSDavid Schultz int 43240dbabfSDavid Schultz __isfinitef(float f) 44240dbabfSDavid Schultz { 45240dbabfSDavid Schultz union IEEEf2bits u; 46240dbabfSDavid Schultz 47240dbabfSDavid Schultz u.f = f; 48240dbabfSDavid Schultz return (u.bits.exp != 255); 49240dbabfSDavid Schultz } 50240dbabfSDavid Schultz 51240dbabfSDavid Schultz int 52240dbabfSDavid Schultz __isfinitel(long double e) 53240dbabfSDavid Schultz { 54240dbabfSDavid Schultz union IEEEl2bits u; 55240dbabfSDavid Schultz 56240dbabfSDavid Schultz u.e = e; 57240dbabfSDavid Schultz return (u.bits.exp != 32767); 58240dbabfSDavid Schultz } 59