Print this page
11210 libm should be cstyle(1ONBLD) clean
@@ -20,10 +20,11 @@
*/
/*
* Copyright 2011 Nexenta Systems, Inc. All rights reserved.
*/
+
/*
* Copyright 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
@@ -31,39 +32,44 @@
#include "libm.h"
#include "longdouble.h"
static const double d_one = 1.0;
-
long double
-cbrtl(long double x) {
+cbrtl(long double x)
+{
long double s, t, r, w, y;
double dx, dy;
- int *py = (int *) &dy;
+ int *py = (int *)&dy;
int n, m, m3, n0, sx;
if (!finitel(x))
return (x + x);
+
if (iszerol(x))
return (x);
+
n0 = 0;
- if (*((int *) &d_one) == 0)
+
+ if (*((int *)&d_one) == 0)
n0 = 1;
+
sx = signbitl(x);
x = fabsl(x);
n = ilogbl(x);
m = n / 3;
m3 = m + m + m;
y = scalbnl(x, -m3);
- dx = (double) y;
+ dx = (double)y;
dy = cbrt(dx);
py[1 - n0] += 2;
+
if (py[1 - n0] == 0)
py[n0] += 1;
/* one step newton iteration to 113 bits with error < 0.667ulps */
- t = (long double) dy;
+ t = (long double)dy;
t = scalbnl(t, m);
s = t * t;
r = x / s;
w = t + t;
r = (r - t) / (w + r);