blob: bd40ef1a491aca36b64d6ed646b0d89095eca72d (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
|
/*
* CDDL HEADER START
*
* The contents of this file are subject to the terms of the
* Common Development and Distribution License (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 2008 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
/* Copyright (c) 1988 AT&T */
/* All Rights Reserved */
#pragma ident "%Z%%M% %I% %E% SMI"
/* IEEE recommended functions */
#pragma weak _finite = finite
#pragma weak _fpclass = fpclass
#pragma weak _unordered = unordered
#include "lint.h"
#include <values.h>
#include "fpparts.h"
#define P754_NOFAULT 1 /* avoid generating extra code */
#include <ieeefp.h>
/*
* FINITE(X)
* finite(x) returns 1 if x > -inf and x < +inf and 0 otherwise
* NaN returns 0
*/
int
finite(double x)
{
return ((EXPONENT(x) != MAXEXP));
}
/*
* UNORDERED(x,y)
* unordered(x,y) returns 1 if x is unordered with y, otherwise
* it returns 0; x is unordered with y if either x or y is NAN
*/
int
unordered(double x, double y)
{
if ((EXPONENT(x) == MAXEXP) && (HIFRACTION(x) || LOFRACTION(x)))
return (1);
if ((EXPONENT(y) == MAXEXP) && (HIFRACTION(y) || LOFRACTION(y)))
return (1);
return (0);
}
/*
* FPCLASS(X)
* fpclass(x) returns the floating point class x belongs to
*/
fpclass_t
fpclass(double x)
{
int sign, exp;
exp = EXPONENT(x);
sign = SIGNBIT(x);
if (exp == 0) { /* de-normal or zero */
if (HIFRACTION(x) || LOFRACTION(x)) /* de-normal */
return (sign ? FP_NDENORM : FP_PDENORM);
else
return (sign ? FP_NZERO : FP_PZERO);
}
if (exp == MAXEXP) { /* infinity or NaN */
if ((HIFRACTION(x) == 0) && (LOFRACTION(x) == 0)) /* infinity */
return (sign ? FP_NINF : FP_PINF);
else
if (QNANBIT(x))
/* hi-bit of mantissa set - quiet nan */
return (FP_QNAN);
else return (FP_SNAN);
}
/* if we reach here we have non-zero normalized number */
return (sign ? FP_NNORM : FP_PNORM);
}
|
