blob: f1d042ff5fa6dedc75814f3a8ae87e01a2cc20cc (
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
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
/*
* 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 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
#pragma ident "@(#)nextafterl.c 1.3 06/01/31 SMI"
#if defined(ELFOBJ)
#pragma weak nextafterl = __nextafterl
#endif
#include "libm.h"
#include <float.h> /* LDBL_MAX, LDBL_MIN */
#if defined(__sparc)
#define n0 0
#define n1 1
#define n2 2
#define n3 3
#define X86PDNRM1(x)
#define INC(px) { \
if (++px[n3] == 0) \
if (++px[n2] == 0) \
if (++px[n1] == 0) \
++px[n0]; \
}
#define DEC(px) { \
if (--px[n3] == 0xffffffff) \
if (--px[n2] == 0xffffffff) \
if (--px[n1] == 0xffffffff) \
--px[n0]; \
}
#elif defined(__i386)
#define n0 2
#define n1 1
#define n2 0
#define n3 0
/*
* if pseudo-denormal, replace by the equivalent normal
*/
#define X86PDNRM1(x) if (XBIASED_EXP(x) == 0 && (((int *) &x)[1] & \
0x80000000) != 0) \
((int *) &x)[2] |= 1
#define INC(px) { \
if (++px[n2] == 0) \
if ((++px[n1] & ~0x80000000) == 0) \
px[n1] = 0x80000000, ++px[n0]; \
}
#define DEC(px) { \
if (--px[n2] == 0xffffffff) \
if (--px[n1] == 0x7fffffff) \
if ((--px[n0] & 0x7fff) != 0) \
px[n1] |= 0x80000000; \
}
#endif
long double
nextafterl(long double x, long double y) {
int *px = (int *) &x;
int *py = (int *) &y;
if (x == y)
return (y); /* C99 requirement */
if (x != x || y != y)
return (x * y);
if (ISZEROL(x)) { /* x == 0.0 */
px[n0] = py[n0] & XSGNMSK;
px[n1] = px[n2] = 0;
px[n3] = 1;
} else {
X86PDNRM1(x);
if ((px[n0] & XSGNMSK) == 0) { /* x > 0.0 */
if (x > y) /* x > y */
DEC(px)
else
INC(px)
} else {
if (x < y) /* x < y */
DEC(px)
else
INC(px)
}
}
#ifndef lint
{
volatile long double dummy;
int k = XBIASED_EXP(x);
if (k == 0)
dummy = LDBL_MIN * copysignl(LDBL_MIN, x);
else if (k == 0x7fff)
dummy = LDBL_MAX * copysignl(LDBL_MAX, x);
}
#endif
return (x);
}
|
