Whitespace.

* Rename the polynomial coefficients from P2, P3, ... to A2, A3, ....
The names now coincide with the name used in PTP Tang's paper.

* Rename the variable from s to tbl to better reflect that
th

* Rename the polynomial coefficients from P2, P3, ... to A2, A3, ....
The names now coincide with the name used in PTP Tang's paper.

* Rename the variable from s to tbl to better reflect that
this is a table, and to be consistent with the naming scheme
in s_exp2l.c

Reviewed by: bde (as part of larger diff)

show more ...

* Style(9). Start non-Copyright fancy formatted comments with /**.

Reviewed by: bde (as part of larger diff)

ld80/s_expl.c:

* Update Copyright years to include 2013.

ld128/s_expl.c:

* Correct and update Copyright years. This code originated from
the ld80 version, so it should reflect the same time per

ld80/s_expl.c:

* Update Copyright years to include 2013.

ld128/s_expl.c:

* Correct and update Copyright years. This code originated from
the ld80 version, so it should reflect the same time period.

Reviewed by: bde (as part of larger diff)

show more ...

MFC

MFC @ r241285

IFC @ r242684

Sync from head

* Update the comment that explains the choice of values in the
table and the requirement on trailing zero bits.

* Remove the __aligned() compiler directives as these were found
to have a negativ

* Update the comment that explains the choice of values in the
table and the requirement on trailing zero bits.

* Remove the __aligned() compiler directives as these were found
to have a negative effect on the produced code.

Submitted by: bde
Approved by: das (mentor)

show more ...

* src/math_private.h:
. Change the API for the LD80C by removing the explicit passing
of the sign bit. The sign can be determined from the last
parameter of the macro.
. On i386, load lo

* src/math_private.h:
. Change the API for the LD80C by removing the explicit passing
of the sign bit. The sign can be determined from the last
parameter of the macro.
. On i386, load long double by bit manipulations to work around
at least a gcc compiler issue. On non-i386 ld80 architectures,
use a simple assignment.

* ld80/s_expl.c:
. Update the only consumer of LD80C.

Submitted by: bde
Approved by: das (mentor)

show more ...

* ld80/s_expl.c:
. Fix the threshold for expl(x) where |x| is small.
. Also update the previously incorrect comment to match the
new threshold.

* ld128/s_expl.c:
. Re-order logic in except

* ld80/s_expl.c:
. Fix the threshold for expl(x) where |x| is small.
. Also update the previously incorrect comment to match the
new threshold.

* ld128/s_expl.c:
. Re-order logic in exceptional cases to match the logic used in
other long double functions.
. Fix the threshold for expl(x) where is |x| is small.
. Also update the previously incorrect comment to match the
new threshold.

Submitted by: bde
Approved by: das (mentor)

show more ...

Fix whitespace issue.

Approved by: das (mentor, implicit)

* ld80/s_expl.c:
. Guard a comment from reformatting by indent(1).
. Re-order variables in declarations to alphabetical order.
. Remove a banal comment.

* ld128/s_expl.c:
. Add a comment to

* ld80/s_expl.c:
. Guard a comment from reformatting by indent(1).
. Re-order variables in declarations to alphabetical order.
. Remove a banal comment.

* ld128/s_expl.c:
. Add a comment to point to ld80/s_expl.c for implementation details.
. Move the #define of INTERVAL to reduce the diff with ld80/s_expl.c.
. twom10000 does not need to be volatile, so move its declaration.
. Re-order variables in declarations to alphabetical order.
. Add a comment that describes the argument reduction.
. Remove the same banal comment found in ld80/s_expl.c.

Reviewed by: bde
Approved by: das (mentor)

show more ...

* Update the lookup table to use 53-bit high and low values.
Also, update the comment to describe the choice of using
a high and low decomposition of 2^(i/INTERNVAL) for
0 <= i <= INTERVAL in p

* Update the lookup table to use 53-bit high and low values.
Also, update the comment to describe the choice of using
a high and low decomposition of 2^(i/INTERNVAL) for
0 <= i <= INTERVAL in preparation for an implementation of
expm1l.

* Move the #define of INTERVAL above the comment, because the
comment refers to INTERVAL.

Reviewed by: bde
Approved by: das (mentor)

show more ...

Merge head r233826 through r239010.

Whitespace.

Submitted by: bde
Approved by: das (pre-approved)

Replace the macro name NUM with INTERVALS. This change provides
compatibility with the INTERVALS macro used in the soon-to-be-commmitted
expm1l() and someday-to-be-committed log*l() functions.

Add

Replace the macro name NUM with INTERVALS. This change provides
compatibility with the INTERVALS macro used in the soon-to-be-commmitted
expm1l() and someday-to-be-committed log*l() functions.

Add a comment into ld128/s_expl.c noting at gcc issue that was
deleted when rewriting ld80/e_expl.c as ld128/s_expl.c.

Requested by: bde
Approved by: das (mentor)

show more ...

* ld80/expl.c:
. Remove a few #ifdefs that should have been removed in the initial
commit.
. Sort fpmath.h to its rightful place.

* ld128/s_expl.c:
. Replace EXPMASK with its actual value.

* ld80/expl.c:
. Remove a few #ifdefs that should have been removed in the initial
commit.
. Sort fpmath.h to its rightful place.

* ld128/s_expl.c:
. Replace EXPMASK with its actual value.
. Sort fpmath.h to its rightful place.

Requested by: bde
Approved by: das (mentor)

show more ...

Compute the exponential of x for Intel 80-bit format and IEEE 128-bit
format. These implementations are based on

PTP Tang, "Table-driven implementation of the exponential function
in IEEE floating-

Compute the exponential of x for Intel 80-bit format and IEEE 128-bit
format. These implementations are based on

PTP Tang, "Table-driven implementation of the exponential function
in IEEE floating-point arithmetic," ACM Trans. Math. Soft., 15,
144-157 (1989).

PR: standards/152415
Submitted by: kargl
Reviewed by: bde, das
Approved by: das (mentor)

show more ...