blob: 588f90cee25b67af69040f3a0e3d552530380873 (
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
|
#include <math.h>
#include <errno.h>
#define _RESEARCH_SOURCE
#include <float.h>
#define MASK 0x7ffL
#define SHIFT 20
#define BIAS 1022L
typedef union
{
double d;
struct
{
#ifdef IEEE_8087
long ls;
long ms;
#else
long ms;
long ls;
#endif
};
} Cheat;
double
frexp(double d, int *ep)
{
Cheat x;
if(d == 0) {
*ep = 0;
return 0;
}
x.d = d;
*ep = ((x.ms >> SHIFT) & MASK) - BIAS;
x.ms &= ~(MASK << SHIFT);
x.ms |= BIAS << SHIFT;
return x.d;
}
double
ldexp(double d, int e)
{
Cheat x;
if(d == 0)
return 0;
x.d = d;
e += (x.ms >> SHIFT) & MASK;
if(e <= 0)
return 0;
if(e >= MASK){
errno = ERANGE;
if(d < 0)
return -HUGE_VAL;
return HUGE_VAL;
}
x.ms &= ~(MASK << SHIFT);
x.ms |= (long)e << SHIFT;
return x.d;
}
double
modf(double d, double *ip)
{
double f;
Cheat x;
int e;
if(d < 1) {
if(d < 0) {
f = modf(-d, ip);
*ip = -*ip;
return -f;
}
*ip = 0;
return d;
}
x.d = d;
e = ((x.ms >> SHIFT) & MASK) - BIAS;
if(e <= SHIFT+1) {
x.ms &= ~(0x1fffffL >> e);
x.ls = 0;
} else
if(e <= SHIFT+33)
x.ls &= ~(0x7fffffffL >> (e-SHIFT-2));
*ip = x.d;
return d - x.d;
}
|
