| OLD | NEW |
| 1 #include "libm.h" | 1 #include "libm.h" |
| 2 | 2 |
| 3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 | 3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 |
| 4 long double nextafterl(long double x, long double y) | 4 long double nextafterl(long double x, long double y) { |
| 5 { | 5 return nextafter(x, y); |
| 6 » return nextafter(x, y); | |
| 7 } | 6 } |
| 8 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 | 7 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 |
| 9 long double nextafterl(long double x, long double y) | 8 long double nextafterl(long double x, long double y) { |
| 10 { | 9 union ldshape ux, uy; |
| 11 » union ldshape ux, uy; | |
| 12 | 10 |
| 13 » if (isnan(x) || isnan(y)) | 11 if (isnan(x) || isnan(y)) |
| 14 » » return x + y; | 12 return x + y; |
| 15 » if (x == y) | 13 if (x == y) |
| 16 » » return y; | 14 return y; |
| 17 » ux.f = x; | 15 ux.f = x; |
| 18 » if (x == 0) { | 16 if (x == 0) { |
| 19 » » uy.f = y; | 17 uy.f = y; |
| 20 » » ux.i.m = 1; | 18 ux.i.m = 1; |
| 21 » » ux.i.se = uy.i.se & 0x8000; | 19 ux.i.se = uy.i.se & 0x8000; |
| 22 » } else if ((x < y) == !(ux.i.se & 0x8000)) { | 20 } else if ((x < y) == !(ux.i.se & 0x8000)) { |
| 23 » » ux.i.m++; | 21 ux.i.m++; |
| 24 » » if (ux.i.m << 1 == 0) { | 22 if (ux.i.m << 1 == 0) { |
| 25 » » » ux.i.m = 1ULL << 63; | 23 ux.i.m = 1ULL << 63; |
| 26 » » » ux.i.se++; | 24 ux.i.se++; |
| 27 » » } | 25 } |
| 28 » } else { | 26 } else { |
| 29 » » if (ux.i.m << 1 == 0) { | 27 if (ux.i.m << 1 == 0) { |
| 30 » » » ux.i.se--; | 28 ux.i.se--; |
| 31 » » » if (ux.i.se) | 29 if (ux.i.se) |
| 32 » » » » ux.i.m = 0; | 30 ux.i.m = 0; |
| 33 » » } | 31 } |
| 34 » » ux.i.m--; | 32 ux.i.m--; |
| 35 » } | 33 } |
| 36 » /* raise overflow if ux is infinite and x is finite */ | 34 /* raise overflow if ux is infinite and x is finite */ |
| 37 » if ((ux.i.se & 0x7fff) == 0x7fff) | 35 if ((ux.i.se & 0x7fff) == 0x7fff) |
| 38 » » return x + x; | 36 return x + x; |
| 39 » /* raise underflow if ux is subnormal or zero */ | 37 /* raise underflow if ux is subnormal or zero */ |
| 40 » if ((ux.i.se & 0x7fff) == 0) | 38 if ((ux.i.se & 0x7fff) == 0) |
| 41 » » FORCE_EVAL(x*x + ux.f*ux.f); | 39 FORCE_EVAL(x * x + ux.f * ux.f); |
| 42 » return ux.f; | 40 return ux.f; |
| 43 } | 41 } |
| 44 #elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384 | 42 #elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384 |
| 45 long double nextafterl(long double x, long double y) | 43 long double nextafterl(long double x, long double y) { |
| 46 { | 44 union ldshape ux, uy; |
| 47 » union ldshape ux, uy; | |
| 48 | 45 |
| 49 » if (isnan(x) || isnan(y)) | 46 if (isnan(x) || isnan(y)) |
| 50 » » return x + y; | 47 return x + y; |
| 51 » if (x == y) | 48 if (x == y) |
| 52 » » return y; | 49 return y; |
| 53 » ux.f = x; | 50 ux.f = x; |
| 54 » if (x == 0) { | 51 if (x == 0) { |
| 55 » » uy.f = y; | 52 uy.f = y; |
| 56 » » ux.i.lo = 1; | 53 ux.i.lo = 1; |
| 57 » » ux.i.se = uy.i.se & 0x8000; | 54 ux.i.se = uy.i.se & 0x8000; |
| 58 » } else if ((x < y) == !(ux.i.se & 0x8000)) { | 55 } else if ((x < y) == !(ux.i.se & 0x8000)) { |
| 59 » » ux.i2.lo++; | 56 ux.i2.lo++; |
| 60 » » if (ux.i2.lo == 0) | 57 if (ux.i2.lo == 0) |
| 61 » » » ux.i2.hi++; | 58 ux.i2.hi++; |
| 62 » } else { | 59 } else { |
| 63 » » if (ux.i2.lo == 0) | 60 if (ux.i2.lo == 0) |
| 64 » » » ux.i2.hi--; | 61 ux.i2.hi--; |
| 65 » » ux.i2.lo--; | 62 ux.i2.lo--; |
| 66 » } | 63 } |
| 67 » /* raise overflow if ux is infinite and x is finite */ | 64 /* raise overflow if ux is infinite and x is finite */ |
| 68 » if ((ux.i.se & 0x7fff) == 0x7fff) | 65 if ((ux.i.se & 0x7fff) == 0x7fff) |
| 69 » » return x + x; | 66 return x + x; |
| 70 » /* raise underflow if ux is subnormal or zero */ | 67 /* raise underflow if ux is subnormal or zero */ |
| 71 » if ((ux.i.se & 0x7fff) == 0) | 68 if ((ux.i.se & 0x7fff) == 0) |
| 72 » » FORCE_EVAL(x*x + ux.f*ux.f); | 69 FORCE_EVAL(x * x + ux.f * ux.f); |
| 73 » return ux.f; | 70 return ux.f; |
| 74 } | 71 } |
| 75 #endif | 72 #endif |
| OLD | NEW |