Update a comment to reflect that we are using an endpoint of

an interval instead of a midpoint.
This commit is contained in:
Steve Kargl 2013-06-03 18:14:18 +00:00
parent ad36b00fcb
commit f3049ab5f3
2 changed files with 2 additions and 2 deletions

View File

@ -226,7 +226,7 @@ expl(long double x)
return (1.0L + x);
}
/* Reduce x to (k*ln2 + midpoint[n2] + r1 + r2). */
/* Reduce x to (k*ln2 + endpoint[n2] + r1 + r2). */
fn = x * INV_L + 0x1.8p112 - 0x1.8p112;
n = (int)fn;
n2 = (unsigned)n % INTERVALS;

View File

@ -261,7 +261,7 @@ expl(long double x)
ENTERI();
/* Reduce x to (k*ln2 + midpoint[n2] + r1 + r2). */
/* Reduce x to (k*ln2 + endpoint[n2] + r1 + r2). */
/* Use a specialized rint() to get fn. Assume round-to-nearest. */
fn = x * INV_L + 0x1.8p63 - 0x1.8p63;
r = x - fn * L1 - fn * L2; /* r = r1 + r2 done independently. */