Last change
on this file since 367 was 337, checked in by coas-nagasima, 6 years ago |
ASP3版ECNLを追加
|
-
Property svn:eol-style
set to
native
-
Property svn:mime-type
set to
text/x-csrc;charset=UTF-8
|
File size:
374 bytes
|
Rev | Line | |
---|
[337] | 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.