From e2014c1dcb07666d51c07fb90d794554e384d85e Mon Sep 17 00:00:00 2001 From: "Jordan K. Hubbard" Date: Fri, 19 Aug 1994 23:52:29 +0000 Subject: [PATCH] Latest fix from jtc: The fyl2xp1 instruction has such a limited range: -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1 it's not worth trying to use it. Also, I'm not sure fyl2xp1's extra precision will matter once the result is converted from extended real (80 bits) back to double real (64 bits). Reviewed by: jkh Submitted by: jtc --- lib/msun/i387/s_log1p.S | 15 +++++++++++++-- 1 file changed, 13 insertions(+), 2 deletions(-) diff --git a/lib/msun/i387/s_log1p.S b/lib/msun/i387/s_log1p.S index 4f144513ed8..ad511bee358 100644 --- a/lib/msun/i387/s_log1p.S +++ b/lib/msun/i387/s_log1p.S @@ -35,10 +35,21 @@ #include -RCSID("$Id: s_log1p.S,v 1.1.1.1 1994/08/19 09:40:01 jkh Exp $") +RCSID("$Id: s_log1p.S,v 1.4 1994/08/19 03:55:57 jtc Exp $") +/* + * The fyl2xp1 instruction has such a limited range: + * -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1 + * it's not worth trying to use it. + * + * Also, I'm not sure fyl2xp1's extra precision will + * matter once the result is converted from extended + * real (80 bits) back to double real (64 bits). + */ ENTRY(log1p) fldln2 fldl 4(%esp) - fyl2xp1 + fld1 + faddp + fyl2x ret