mirrors/phobos
1
0
Fork 0
mirror of https://github.com/dlang/phobos.git synced 2025-04-28 14:10:30 +03:00
Code Issues Projects Releases Packages Wiki Activity Actions

Merge pull request #4975 from somzzz/issue_16634

Browse source 
fix issue 16634 -  std.math exposes yl2x and yl2xp1 publicly
This commit is contained in:
Andrei Alexandrescu 2016-12-21 21:18:05 -05:00 • committed by GitHub
commit 6acf4a289d
1 changed files with 6 additions and 23 deletions
Download patch file Download diff file Expand all files Collapse all files

29
std/math.d
View file

@ -3020,7 +3020,7 @@ typed_allocator.d
real log(real x) @safe pure nothrow @nogc real log(real x) @safe pure nothrow @nogc
{ {
version (INLINE_YL2X) version (INLINE_YL2X)
return yl2x(x, LN2); return core.math.yl2x(x, LN2);
else else
{ {
// Coefficients for log(1 + x) // Coefficients for log(1 + x)
@ -3147,7 +3147,7 @@ real log(real x) @safe pure nothrow @nogc
real log10(real x) @safe pure nothrow @nogc real log10(real x) @safe pure nothrow @nogc
{ {
version (INLINE_YL2X) version (INLINE_YL2X)
return yl2x(x, LOG2); return core.math.yl2x(x, LOG2);
else else
{ {
// Coefficients for log(1 + x) // Coefficients for log(1 + x)
@ -3285,7 +3285,7 @@ real log1p(real x) @safe pure nothrow @nogc
{ {
// On x87, yl2xp1 is valid if and only if -0.5 <= lg(x) <= 0.5, // On x87, yl2xp1 is valid if and only if -0.5 <= lg(x) <= 0.5,
// ie if -0.29<=x<=0.414 // ie if -0.29<=x<=0.414
return (fabs(x) <= 0.25) ? yl2xp1(x, LN2) : yl2x(x+1, LN2); return (fabs(x) <= 0.25) ? core.math.yl2xp1(x, LN2) : core.math.yl2x(x+1, LN2);
} }
else else
{ {
@ -3317,7 +3317,7 @@ real log1p(real x) @safe pure nothrow @nogc
real log2(real x) @safe pure nothrow @nogc real log2(real x) @safe pure nothrow @nogc
{ {
version (INLINE_YL2X) version (INLINE_YL2X)
return yl2x(x, 1); return core.math.yl2x(x, 1);
else else
{ {
// Coefficients for log(1 + x) // Coefficients for log(1 + x)
@ -3579,7 +3579,7 @@ real cbrt(real x) @trusted nothrow @nogc
version (CRuntime_Microsoft) version (CRuntime_Microsoft)
{ {
version (INLINE_YL2X) version (INLINE_YL2X)
return copysign(exp2(yl2x(fabs(x), 1.0L/3.0L)), x); return copysign(exp2(core.math.yl2x(fabs(x), 1.0L/3.0L)), x);
else else
return core.stdc.math.cbrtl(x); return core.stdc.math.cbrtl(x);
} }
@ -6378,7 +6378,7 @@ Unqual!(Largest!(F, G)) pow(F, G)(F x, G y) @nogc @trusted pure nothrow
// (though complicated) method is described in: // (though complicated) method is described in:
// "An efficient rounding boundary test for pow(x, y) // "An efficient rounding boundary test for pow(x, y)
// in double precision", C.Q. Lauter and V. Lefèvre, INRIA (2007). // in double precision", C.Q. Lauter and V. Lefèvre, INRIA (2007).
return sign * exp2( yl2x(x, y) ); return sign * exp2( core.math.yl2x(x, y) );
} }
else else
{ {
@ -7094,23 +7094,6 @@ deprecated("Phobos1 math functions are deprecated, use isNormal ") alias isnorma
deprecated("Phobos1 math functions are deprecated, use isSubnormal ") alias issubnormal = isSubnormal; deprecated("Phobos1 math functions are deprecated, use isSubnormal ") alias issubnormal = isSubnormal;
deprecated("Phobos1 math functions are deprecated, use isInfinity ") alias isinf = isInfinity; deprecated("Phobos1 math functions are deprecated, use isInfinity ") alias isinf = isInfinity;
/* **********************************
* Building block functions, they
* translate to a single x87 instruction.
*/
real yl2x(real x, real y) @nogc @safe pure nothrow; // y * log2(x)
real yl2xp1(real x, real y) @nogc @safe pure nothrow; // y * log2(x + 1)
@safe pure nothrow @nogc unittest
{
version (INLINE_YL2X)
{
assert(yl2x(1024, 1) == 10);
assert(yl2xp1(1023, 1) == 10);
}
}
@safe pure nothrow @nogc unittest @safe pure nothrow @nogc unittest
{ {
real num = real.infinity; real num = real.infinity;