diff options
Diffstat (limited to 'sysdeps/ieee754/ldbl-128/s_cbrtl.c')
-rw-r--r-- | sysdeps/ieee754/ldbl-128/s_cbrtl.c | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/sysdeps/ieee754/ldbl-128/s_cbrtl.c b/sysdeps/ieee754/ldbl-128/s_cbrtl.c index d93684e50e..eb88d29fc9 100644 --- a/sysdeps/ieee754/ldbl-128/s_cbrtl.c +++ b/sysdeps/ieee754/ldbl-128/s_cbrtl.c @@ -57,10 +57,10 @@ Adapted for glibc October, 2001. #include <math.h> #include <math_private.h> -static const _Float128 CBRT2 = 1.259921049894873164767210607278228350570251L; -static const _Float128 CBRT4 = 1.587401051968199474751705639272308260391493L; -static const _Float128 CBRT2I = 0.7937005259840997373758528196361541301957467L; -static const _Float128 CBRT4I = 0.6299605249474365823836053036391141752851257L; +static const _Float128 CBRT2 = L(1.259921049894873164767210607278228350570251); +static const _Float128 CBRT4 = L(1.587401051968199474751705639272308260391493); +static const _Float128 CBRT2I = L(0.7937005259840997373758528196361541301957467); +static const _Float128 CBRT4I = L(0.6299605249474365823836053036391141752851257); _Float128 @@ -89,11 +89,11 @@ __cbrtl (_Float128 x) /* Approximate cube root of number between .5 and 1, peak relative error = 1.2e-6 */ - x = ((((1.3584464340920900529734e-1L * x - - 6.3986917220457538402318e-1L) * x - + 1.2875551670318751538055e0L) * x - - 1.4897083391357284957891e0L) * x - + 1.3304961236013647092521e0L) * x + 3.7568280825958912391243e-1L; + x = ((((L(1.3584464340920900529734e-1) * x + - L(6.3986917220457538402318e-1)) * x + + L(1.2875551670318751538055e0)) * x + - L(1.4897083391357284957891e0)) * x + + L(1.3304961236013647092521e0)) * x + L(3.7568280825958912391243e-1); /* exponent divided by 3 */ if (e >= 0) @@ -123,9 +123,9 @@ __cbrtl (_Float128 x) x = __ldexpl (x, e); /* Newton iteration */ - x -= (x - (z / (x * x))) * 0.3333333333333333333333333333333333333333L; - x -= (x - (z / (x * x))) * 0.3333333333333333333333333333333333333333L; - x -= (x - (z / (x * x))) * 0.3333333333333333333333333333333333333333L; + x -= (x - (z / (x * x))) * L(0.3333333333333333333333333333333333333333); + x -= (x - (z / (x * x))) * L(0.3333333333333333333333333333333333333333); + x -= (x - (z / (x * x))) * L(0.3333333333333333333333333333333333333333); if (sign < 0) x = -x; |