Print this page
11210 libm should be cstyle(1ONBLD) clean

Split Close
Expand all
Collapse all
          --- old/usr/src/lib/libm/common/LD/cbrtl.c
          +++ new/usr/src/lib/libm/common/LD/cbrtl.c
↓ open down ↓ 14 lines elided ↑ open up ↑
  15   15   * If applicable, add the following below this CDDL HEADER, with the
  16   16   * fields enclosed by brackets "[]" replaced with your own identifying
  17   17   * information: Portions Copyright [yyyy] [name of copyright owner]
  18   18   *
  19   19   * CDDL HEADER END
  20   20   */
  21   21  
  22   22  /*
  23   23   * Copyright 2011 Nexenta Systems, Inc.  All rights reserved.
  24   24   */
       25 +
  25   26  /*
  26   27   * Copyright 2006 Sun Microsystems, Inc.  All rights reserved.
  27   28   * Use is subject to license terms.
  28   29   */
  29   30  
  30   31  #pragma weak __cbrtl = cbrtl
  31   32  
  32   33  #include "libm.h"
  33   34  #include "longdouble.h"
  34   35  
  35   36  static const double d_one = 1.0;
  36      -
  37   37  long double
  38      -cbrtl(long double x) {
       38 +cbrtl(long double x)
       39 +{
  39   40          long double s, t, r, w, y;
  40   41          double dx, dy;
  41      -        int *py = (int *) &dy;
       42 +        int *py = (int *)&dy;
  42   43          int n, m, m3, n0, sx;
  43   44  
  44   45          if (!finitel(x))
  45   46                  return (x + x);
       47 +
  46   48          if (iszerol(x))
  47   49                  return (x);
       50 +
  48   51          n0 = 0;
  49      -        if (*((int *) &d_one) == 0)
       52 +
       53 +        if (*((int *)&d_one) == 0)
  50   54                  n0 = 1;
       55 +
  51   56          sx = signbitl(x);
  52   57          x = fabsl(x);
  53   58          n = ilogbl(x);
  54   59          m = n / 3;
  55   60          m3 = m + m + m;
  56   61          y = scalbnl(x, -m3);
  57      -        dx = (double) y;
       62 +        dx = (double)y;
  58   63          dy = cbrt(dx);
  59   64          py[1 - n0] += 2;
       65 +
  60   66          if (py[1 - n0] == 0)
  61   67                  py[n0] += 1;
  62   68  
  63   69          /* one step newton iteration to 113 bits with error < 0.667ulps */
  64      -        t = (long double) dy;
       70 +        t = (long double)dy;
  65   71          t = scalbnl(t, m);
  66   72          s = t * t;
  67   73          r = x / s;
  68   74          w = t + t;
  69   75          r = (r - t) / (w + r);
  70   76          t += t * r;
  71   77  
  72   78          return (sx == 0 ? t : -t);
  73   79  }
    
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX