Use __ieee754_sqrt() instead of sqrt() internally. Similarly for the

float versions.  Using sqrt() was inefficient.

Obtained from:	NetBSD

1 files changed, 3 insertions, 3 deletions

diff --git a/lib/msun/src/e_acoshf.c b/lib/msun/src/e_acoshf.c
index ba4006b..1d234d2 100644
--- a/lib/msun/src/e_acoshf.c
+++ b/lib/msun/src/e_acoshf.c
@@ -14,7 +14,7 @@
  */
 
 #ifndef lint
-static char rcsid[] = "$Id$";
+static char rcsid[] = "$Id: e_acoshf.c,v 1.4 1997/02/22 15:09:56 peter Exp $";
 #endif
 
 #include "math.h"
@@ -49,9 +49,9 @@ ln2	= 6.9314718246e-01;  /* 0x3f317218 */
 	    return 0.0;			/* acosh(1) = 0 */
 	} else if (hx > 0x40000000) {	/* 2**28 > x > 2 */
 	    t=x*x;
-	    return __ieee754_logf((float)2.0*x-one/(x+sqrtf(t-one)));
+	    return __ieee754_logf((float)2.0*x-one/(x+__ieee754_sqrtf(t-one)));
 	} else {			/* 1<x<2 */
 	    t = x-one;
-	    return log1pf(t+sqrtf((float)2.0*t+t*t));
+	    return log1pf(t+__ieee754_sqrtf((float)2.0*t+t*t));
 	}
 }