1*7c478bd9Sstevel@tonic-gate /* 2*7c478bd9Sstevel@tonic-gate * CDDL HEADER START 3*7c478bd9Sstevel@tonic-gate * 4*7c478bd9Sstevel@tonic-gate * The contents of this file are subject to the terms of the 5*7c478bd9Sstevel@tonic-gate * Common Development and Distribution License, Version 1.0 only 6*7c478bd9Sstevel@tonic-gate * (the "License"). You may not use this file except in compliance 7*7c478bd9Sstevel@tonic-gate * with the License. 8*7c478bd9Sstevel@tonic-gate * 9*7c478bd9Sstevel@tonic-gate * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE 10*7c478bd9Sstevel@tonic-gate * or http://www.opensolaris.org/os/licensing. 11*7c478bd9Sstevel@tonic-gate * See the License for the specific language governing permissions 12*7c478bd9Sstevel@tonic-gate * and limitations under the License. 13*7c478bd9Sstevel@tonic-gate * 14*7c478bd9Sstevel@tonic-gate * When distributing Covered Code, include this CDDL HEADER in each 15*7c478bd9Sstevel@tonic-gate * file and include the License file at usr/src/OPENSOLARIS.LICENSE. 16*7c478bd9Sstevel@tonic-gate * If applicable, add the following below this CDDL HEADER, with the 17*7c478bd9Sstevel@tonic-gate * fields enclosed by brackets "[]" replaced with your own identifying 18*7c478bd9Sstevel@tonic-gate * information: Portions Copyright [yyyy] [name of copyright owner] 19*7c478bd9Sstevel@tonic-gate * 20*7c478bd9Sstevel@tonic-gate * CDDL HEADER END 21*7c478bd9Sstevel@tonic-gate */ 22*7c478bd9Sstevel@tonic-gate /* 23*7c478bd9Sstevel@tonic-gate * Copyright 2004 Sun Microsystems, Inc. All rights reserved. 24*7c478bd9Sstevel@tonic-gate * Use is subject to license terms. 25*7c478bd9Sstevel@tonic-gate */ 26*7c478bd9Sstevel@tonic-gate 27*7c478bd9Sstevel@tonic-gate #pragma ident "%Z%%M% %I% %E% SMI" 28*7c478bd9Sstevel@tonic-gate 29*7c478bd9Sstevel@tonic-gate /* Copyright (c) 1988 AT&T */ 30*7c478bd9Sstevel@tonic-gate /* All Rights Reserved */ 31*7c478bd9Sstevel@tonic-gate 32*7c478bd9Sstevel@tonic-gate 33*7c478bd9Sstevel@tonic-gate /* IEEE recommended functions */ 34*7c478bd9Sstevel@tonic-gate 35*7c478bd9Sstevel@tonic-gate #pragma weak finite = _finite 36*7c478bd9Sstevel@tonic-gate #pragma weak fpclass = _fpclass 37*7c478bd9Sstevel@tonic-gate #pragma weak unordered = _unordered 38*7c478bd9Sstevel@tonic-gate 39*7c478bd9Sstevel@tonic-gate #include "synonyms.h" 40*7c478bd9Sstevel@tonic-gate #include <values.h> 41*7c478bd9Sstevel@tonic-gate #include "fpparts.h" 42*7c478bd9Sstevel@tonic-gate 43*7c478bd9Sstevel@tonic-gate #define P754_NOFAULT 1 /* avoid generating extra code */ 44*7c478bd9Sstevel@tonic-gate #include <ieeefp.h> 45*7c478bd9Sstevel@tonic-gate 46*7c478bd9Sstevel@tonic-gate /* 47*7c478bd9Sstevel@tonic-gate * FINITE(X) 48*7c478bd9Sstevel@tonic-gate * finite(x) returns 1 if x > -inf and x < +inf and 0 otherwise 49*7c478bd9Sstevel@tonic-gate * NaN returns 0 50*7c478bd9Sstevel@tonic-gate */ 51*7c478bd9Sstevel@tonic-gate 52*7c478bd9Sstevel@tonic-gate int 53*7c478bd9Sstevel@tonic-gate finite(double x) 54*7c478bd9Sstevel@tonic-gate { 55*7c478bd9Sstevel@tonic-gate return ((EXPONENT(x) != MAXEXP)); 56*7c478bd9Sstevel@tonic-gate } 57*7c478bd9Sstevel@tonic-gate 58*7c478bd9Sstevel@tonic-gate /* 59*7c478bd9Sstevel@tonic-gate * UNORDERED(x,y) 60*7c478bd9Sstevel@tonic-gate * unordered(x,y) returns 1 if x is unordered with y, otherwise 61*7c478bd9Sstevel@tonic-gate * it returns 0; x is unordered with y if either x or y is NAN 62*7c478bd9Sstevel@tonic-gate */ 63*7c478bd9Sstevel@tonic-gate 64*7c478bd9Sstevel@tonic-gate int 65*7c478bd9Sstevel@tonic-gate unordered(double x, double y) 66*7c478bd9Sstevel@tonic-gate { 67*7c478bd9Sstevel@tonic-gate if ((EXPONENT(x) == MAXEXP) && (HIFRACTION(x) || LOFRACTION(x))) 68*7c478bd9Sstevel@tonic-gate return (1); 69*7c478bd9Sstevel@tonic-gate if ((EXPONENT(y) == MAXEXP) && (HIFRACTION(y) || LOFRACTION(y))) 70*7c478bd9Sstevel@tonic-gate return (1); 71*7c478bd9Sstevel@tonic-gate return (0); 72*7c478bd9Sstevel@tonic-gate } 73*7c478bd9Sstevel@tonic-gate 74*7c478bd9Sstevel@tonic-gate /* 75*7c478bd9Sstevel@tonic-gate * FPCLASS(X) 76*7c478bd9Sstevel@tonic-gate * fpclass(x) returns the floating point class x belongs to 77*7c478bd9Sstevel@tonic-gate */ 78*7c478bd9Sstevel@tonic-gate 79*7c478bd9Sstevel@tonic-gate fpclass_t 80*7c478bd9Sstevel@tonic-gate fpclass(double x) 81*7c478bd9Sstevel@tonic-gate { 82*7c478bd9Sstevel@tonic-gate int sign, exp; 83*7c478bd9Sstevel@tonic-gate 84*7c478bd9Sstevel@tonic-gate exp = EXPONENT(x); 85*7c478bd9Sstevel@tonic-gate sign = SIGNBIT(x); 86*7c478bd9Sstevel@tonic-gate if (exp == 0) { /* de-normal or zero */ 87*7c478bd9Sstevel@tonic-gate if (HIFRACTION(x) || LOFRACTION(x)) /* de-normal */ 88*7c478bd9Sstevel@tonic-gate return (sign ? FP_NDENORM : FP_PDENORM); 89*7c478bd9Sstevel@tonic-gate else 90*7c478bd9Sstevel@tonic-gate return (sign ? FP_NZERO : FP_PZERO); 91*7c478bd9Sstevel@tonic-gate } 92*7c478bd9Sstevel@tonic-gate if (exp == MAXEXP) { /* infinity or NaN */ 93*7c478bd9Sstevel@tonic-gate if ((HIFRACTION(x) == 0) && (LOFRACTION(x) == 0)) /* infinity */ 94*7c478bd9Sstevel@tonic-gate return (sign ? FP_NINF : FP_PINF); 95*7c478bd9Sstevel@tonic-gate else 96*7c478bd9Sstevel@tonic-gate if (QNANBIT(x)) 97*7c478bd9Sstevel@tonic-gate /* hi-bit of mantissa set - quiet nan */ 98*7c478bd9Sstevel@tonic-gate return (FP_QNAN); 99*7c478bd9Sstevel@tonic-gate else return (FP_SNAN); 100*7c478bd9Sstevel@tonic-gate } 101*7c478bd9Sstevel@tonic-gate /* if we reach here we have non-zero normalized number */ 102*7c478bd9Sstevel@tonic-gate return (sign ? FP_NNORM : FP_PNORM); 103*7c478bd9Sstevel@tonic-gate } 104