Last change
on this file since 445 was 444, checked in by coas-nagasima, 4 years ago |
muslのソースコードを追加
|
-
Property svn:eol-style
set to
native
-
Property svn:mime-type
set to
text/x-csrc;charset=UTF-8
|
File size:
374 bytes
|
Line | |
---|
1 | #include <math.h>
|
---|
2 | #include <stdint.h>
|
---|
3 |
|
---|
4 | double frexp(double x, int *e)
|
---|
5 | {
|
---|
6 | union { double d; uint64_t i; } y = { x };
|
---|
7 | int ee = y.i>>52 & 0x7ff;
|
---|
8 |
|
---|
9 | if (!ee) {
|
---|
10 | if (x) {
|
---|
11 | x = frexp(x*0x1p64, e);
|
---|
12 | *e -= 64;
|
---|
13 | } else *e = 0;
|
---|
14 | return x;
|
---|
15 | } else if (ee == 0x7ff) {
|
---|
16 | return x;
|
---|
17 | }
|
---|
18 |
|
---|
19 | *e = ee - 0x3fe;
|
---|
20 | y.i &= 0x800fffffffffffffull;
|
---|
21 | y.i |= 0x3fe0000000000000ull;
|
---|
22 | return y.d;
|
---|
23 | }
|
---|
Note:
See
TracBrowser
for help on using the repository browser.