math: use sqrtl if FLT_EVAL_METHOD==2 in acosh and acoshf

this makes acosh slightly more precise around 1.0 on i386

diff --git a/src/math/acosh.c b/src/math/acosh.c
index 4ce9b3d1c252d88ed165357fa4b24d13430a4cc9..badbf9081e644ff9df80bd2ae86cfc32e32a888b 100644 (file)
--- a/src/math/acosh.c
+++ b/src/math/acosh.c
@@ -1,5 +1,10 @@
 #include "libm.h"
 
+#if FLT_EVAL_METHOD==2
+#undef sqrt
+#define sqrt sqrtl
+#endif
+
 /* acosh(x) = log(x + sqrt(x*x-1)) */
 double acosh(double x)
 {

diff --git a/src/math/acoshf.c b/src/math/acoshf.c
index 16550f19d397a5c4b267b00d946cecf340a3d21e..8a4ec4d57ead3cbcd473d2ddfdd43d910956f1b9 100644 (file)
--- a/src/math/acoshf.c
+++ b/src/math/acoshf.c
@@ -1,5 +1,13 @@
 #include "libm.h"
 
+#if FLT_EVAL_METHOD==2
+#undef sqrtf
+#define sqrtf sqrtl
+#elif FLT_EVAL_METHOD==1
+#undef sqrtf
+#define sqrtf sqrt
+#endif
+
 /* acosh(x) = log(x + sqrt(x*x-1)) */
 float acoshf(float x)
 {