Print this page
11210 libm should be cstyle(1ONBLD) clean
*** 20,29 ****
--- 20,30 ----
*/
/*
* Copyright 2011 Nexenta Systems, Inc. All rights reserved.
*/
+
/*
* Copyright 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
*** 31,69 ****
#include "libm.h"
#include "longdouble.h"
static const double d_one = 1.0;
-
long double
! cbrtl(long double x) {
long double s, t, r, w, y;
double dx, 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)
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;
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 = scalbnl(t, m);
s = t * t;
r = x / s;
w = t + t;
r = (r - t) / (w + r);
--- 32,75 ----
#include "libm.h"
#include "longdouble.h"
static const double d_one = 1.0;
long double
! cbrtl(long double x)
! {
long double s, t, r, w, y;
double dx, 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)
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;
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 = scalbnl(t, m);
s = t * t;
r = x / s;
w = t + t;
r = (r - t) / (w + r);