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

Split Close
Expand all
Collapse all
          --- old/usr/src/lib/libm/common/R/__sinf.c
          +++ new/usr/src/lib/libm/common/R/__sinf.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  #include "libm.h"
  31   32  
  32      -/* INDENT OFF */
       33 +
  33   34  /*
  34   35   * float __k_sin(double x);
  35   36   * kernel (float) sin function on [-pi/4, pi/4], pi/4 ~ 0.785398164
  36   37   * Input x is in double and assumed to be bounded by ~pi/4 in magnitude.
  37   38   *
  38   39   * Method: Let z = x * x, then
  39   40   *      S(x) = x(S0 + S1*z)(S2 + S3*z + z*z)
  40   41   * where
  41   42   *      S0 =   1.85735322054308378716204874632872525989806770558e-0003,
  42   43   *      S1 =  -1.95035094218403635082921458859320791358115801259e-0004,
↓ open down ↓ 1 lines elided ↑ open up ↑
  44   45   *      S3 =  -3.31975110777873728964197739157371509422022905947e+0001,
  45   46   *
  46   47   * The remez error is bound by  |(sin(x) - S(x))/x| < 2**(-28.2)
  47   48   *
  48   49   * Constants:
  49   50   * The hexadecimal values are the intended ones for the following constants.
  50   51   * The decimal values may be used, provided that the compiler will convert
  51   52   * from decimal to binary accurately enough to produce the hexadecimal values
  52   53   * shown.
  53   54   */
  54      -/* INDENT ON */
  55   55  
  56   56  static const double q[] = {
  57      -/* S0 = */  1.85735322054308378716204874632872525989806770558e-0003,
       57 +/* S0 = */
       58 +        1.85735322054308378716204874632872525989806770558e-0003,
  58   59  /* S1 = */ -1.95035094218403635082921458859320791358115801259e-0004,
  59      -/* S2 = */  5.38400550766074785970952495168558701485841707252e+0002,
  60      -/* S3 = */ -3.31975110777873728964197739157371509422022905947e+0001,
  61      -};
  62      -
  63      -#define S0  q[0]
  64      -#define S1  q[1]
  65      -#define S2  q[2]
  66      -#define S3  q[3]
       60 +/* S2 = */ 5.38400550766074785970952495168558701485841707252e+0002,
       61 +/* S3 = */ -3.31975110777873728964197739157371509422022905947e+0001, };
       62 +
       63 +#define S0              q[0]
       64 +#define S1              q[1]
       65 +#define S2              q[2]
       66 +#define S3              q[3]
  67   67  
  68   68  float
  69      -__k_sinf(double x) {
       69 +__k_sinf(double x)
       70 +{
  70   71          float ft;
  71   72          double z;
  72   73          int hx;
  73   74  
  74      -        hx = ((int *) &x)[HIWORD];      /* hx = leading x */
       75 +        hx = ((int *)&x)[HIWORD];               /* hx = leading x */
       76 +
  75   77          if ((hx & ~0x80000000) < 0x3f100000) {  /* if |x| < 2**-14 */
  76      -                ft = (float) x;
  77      -                if ((int) x == 0)       /* raise inexact if x != 0 */
       78 +                ft = (float)x;
       79 +
       80 +                if ((int)x == 0)                /* raise inexact if x != 0 */
  78   81                          return (ft);
  79   82          }
       83 +
  80   84          z = x * x;
  81      -        ft = (float) ((x * (S0 + z * S1)) * (S2 + z * (S3 + z)));
       85 +        ft = (float)((x * (S0 + z * S1)) * (S2 + z * (S3 + z)));
  82   86          return (ft);
  83   87  }
    
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX