mirrors/phobos
1
0
Fork 0
mirror of https://github.com/dlang/phobos.git synced 2025-05-12 23:29:01 +03:00
Code Issues Projects Releases Packages Wiki Activity Actions 1

std.math: Remove special cases that attempt to raise FE overflow/underflow flag.

Browse source
This commit is contained in:
Iain Buclaw 2018-03-24 10:42:23 +01:00
parent f55a6a797d
commit d79b27eab3
1 changed files with 8 additions and 36 deletions
Download patch file Download diff file Expand all files Collapse all files

44
std/math.d
View file

@ -1834,24 +1834,13 @@ real exp(real x) @trusted pure nothrow @nogc
else
static assert(0, "Not implemented for this architecture");
// Special cases. Raises an overflow or underflow flag accordingly,
// except in the case for CTFE, where there are no hardware controls.
// Special cases.
if (isNaN(x))
return x;
if (x > OF)
{
if (__ctfe)
return real.infinity;
else
return real.max * copysign(real.max, real.infinity);
}
return real.infinity;
if (x < UF)
{
if (__ctfe)
return 0.0;
else
return real.min_normal * copysign(real.min_normal, 0.0);
}
return 0.0;
// Express: e^^x = e^^g * 2^^n
// = e^^g * e^^(n * LOG2E)
@ -2121,15 +2110,9 @@ L_largenegative:
enum real C1 = 6.9314575195312500000000E-1L;
enum real C2 = 1.428606820309417232121458176568075500134E-6L;
// Special cases. Raises an overflow flag, except in the case
// for CTFE, where there are no hardware controls.
// Special cases.
if (x > OF)
{
if (__ctfe)
return real.infinity;
else
return real.max * copysign(real.max, real.infinity);
}
return real.infinity;
if (x == 0.0)
return x;
if (x < UF)
@ -2419,24 +2402,13 @@ private real exp2Impl(real x) @nogc @trusted pure nothrow
enum real OF = 16_384.0L;
enum real UF = -16_382.0L;
// Special cases. Raises an overflow or underflow flag accordingly,
// except in the case for CTFE, where there are no hardware controls.
// Special cases.
if (isNaN(x))
return x;
if (x > OF)
{
if (__ctfe)
return real.infinity;
else
return real.max * copysign(real.max, real.infinity);
}
return real.infinity;
if (x < UF)
{
if (__ctfe)
return 0.0;
else
return real.min_normal * copysign(real.min_normal, 0.0);
}
return 0.0;
// Separate into integer and fractional parts.
int n = cast(int) floor(x + 0.5);