summaryrefslogtreecommitdiff
path: root/libquadmath/math/cbrtq.c
diff options
context:
space:
mode:
authorjakub <jakub@138bc75d-0d04-0410-961f-82ee72b054a4>2012-12-13 12:35:56 +0000
committerjakub <jakub@138bc75d-0d04-0410-961f-82ee72b054a4>2012-12-13 12:35:56 +0000
commit3f11fdf1afaf467ae428801fb96b4f23a9a2c1c0 (patch)
treef6549255d38c82758564c06a82765b9e0d30bb92 /libquadmath/math/cbrtq.c
parentcc4b03ac450a259ceaba67e684b22c5ef621d1d0 (diff)
downloadgcc-3f11fdf1afaf467ae428801fb96b4f23a9a2c1c0.tar.gz
* math/cbrtq.c (cbrtq): Use Q suffixed floating point constants
instead of L suffixed ones. * math/fmaq.c (fmaq): Likewise. * math/rintq.c (TWO112): Likewise. git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@194475 138bc75d-0d04-0410-961f-82ee72b054a4
Diffstat (limited to 'libquadmath/math/cbrtq.c')
-rw-r--r--libquadmath/math/cbrtq.c16
1 files changed, 8 insertions, 8 deletions
diff --git a/libquadmath/math/cbrtq.c b/libquadmath/math/cbrtq.c
index f1f05cac789..2567d4d5bd3 100644
--- a/libquadmath/math/cbrtq.c
+++ b/libquadmath/math/cbrtq.c
@@ -88,11 +88,11 @@ cbrtq ( __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 = ((((1.3584464340920900529734e-1Q * x
+ - 6.3986917220457538402318e-1Q) * x
+ + 1.2875551670318751538055e0Q) * x
+ - 1.4897083391357284957891e0Q) * x
+ + 1.3304961236013647092521e0Q) * x + 3.7568280825958912391243e-1Q;
/* exponent divided by 3 */
if (e >= 0)
@@ -122,9 +122,9 @@ cbrtq ( __float128 x)
x = ldexpq (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))) * 0.3333333333333333333333333333333333333333Q;
+ x -= (x - (z / (x * x))) * 0.3333333333333333333333333333333333333333Q;
+ x -= (x - (z / (x * x))) * 0.3333333333333333333333333333333333333333Q;
if (sign < 0)
x = -x;