| /programs/develop/libraries/menuetlibc/src/libc/ansi/math/ldexp.c |
|---|
| 0,0 → 1,32 |
| /* Copyright (C) 1996 DJ Delorie, see COPYING.DJ for details */ |
| /* Copyright (C) 1995 DJ Delorie, see COPYING.DJ for details */ |
| #include <math.h> |
| double |
| ldexp(double v, int e) |
| { |
| double two = 2.0; |
| if (e < 0) |
| { |
| e = -e; /* This just might overflow on two-complement machines. */ |
| if (e < 0) return 0.0; |
| while (e > 0) |
| { |
| if (e & 1) v /= two; |
| two *= two; |
| e >>= 1; |
| } |
| } |
| else if (e > 0) |
| { |
| while (e > 0) |
| { |
| if (e & 1) v *= two; |
| two *= two; |
| e >>= 1; |
| } |
| } |
| return v; |
| } |