fix potentially wrong-sign zero in cproj functions at infinity

[musl] / src / math / __polevll.c

diff --git a/src/math/__polevll.c b/src/math/__polevll.c
index 08e68d4..ce1a840 100644 (file)
--- a/src/math/__polevll.c
+++ b/src/math/__polevll.c
@@ -56,11 +56,13 @@
 
 #include "libm.h"
 
+#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
+#else
 /*
  * Polynomial evaluator:
  *  P[0] x^n  +  P[1] x^(n-1)  +  ...  +  P[n]
  */
-long double __polevll(long double x, long double *P, int n)
+long double __polevll(long double x, const long double *P, int n)
 {
        long double y;
 
@@ -76,7 +78,7 @@ long double __polevll(long double x, long double *P, int n)
  * Polynomial evaluator:
  *  x^n  +  P[0] x^(n-1)  +  P[1] x^(n-2)  +  ...  +  P[n]
  */
-long double __p1evll(long double x, long double *P, int n)
+long double __p1evll(long double x, const long double *P, int n)
 {
        long double y;
 
@@ -88,3 +90,4 @@ long double __p1evll(long double x, long double *P, int n)
 
        return y;
 }
+#endif