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

this makes acosh slightly more precise around 1.0 on i386

math: use sqrtl if FLT_EVAL_METHOD==2 in acosh and acoshf
this makes acosh slightly more precise around 1.0 on i386
4b539a82 · Szabolcs Nagy · 551da2c1 · 4b539a82 · 4b539a82
隐藏空白更改
内联并排

Showing with 13 addition and 0 deletion

src/math/acosh.c src/math/acosh.c +5 -0

src/math/acoshf.c src/math/acoshf.c +8 -0

未找到文件。
--- a/src/math/acosh.c
+++ b/src/math/acosh.c
 #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)
 {

--- a/src/math/acoshf.c
+++ b/src/math/acoshf.c
 #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)
 {