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

*** 20,37 **** */ /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #include "libm.h" ! /* INDENT OFF */ /* * float __k_tan(double x); * kernel (float) tan function on [-pi/4, pi/4], pi/4 ~ 0.785398164 * Input x is in double and assumed to be bounded by ~pi/4 in magnitude. * --- 20,38 ---- */ /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ + /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #include "libm.h" ! /* * float __k_tan(double x); * kernel (float) tan function on [-pi/4, pi/4], pi/4 ~ 0.785398164 * Input x is in double and assumed to be bounded by ~pi/4 in magnitude. *
*** 41,51 **** * from decimal to binary accurately enough to produce the hexadecimal values * shown. */ static const double q[] = { ! /* one */ 1.0, /* P0 */ 4.46066928428959230679140546271810308098793029785e-0003, /* P1 */ 4.92165316309189027066395283327437937259674072266e+0000, /* P2 */ -7.11410648161473480044492134766187518835067749023e-0001, /* P3 */ 4.08549808374053391446523164631798863410949707031e+0000, /* P4 */ 2.50411070398050927821032018982805311679840087891e+0000, --- 42,53 ---- * from decimal to binary accurately enough to produce the hexadecimal values * shown. */ static const double q[] = { ! /* one */ ! 1.0, /* P0 */ 4.46066928428959230679140546271810308098793029785e-0003, /* P1 */ 4.92165316309189027066395283327437937259674072266e+0000, /* P2 */ -7.11410648161473480044492134766187518835067749023e-0001, /* P3 */ 4.08549808374053391446523164631798863410949707031e+0000, /* P4 */ 2.50411070398050927821032018982805311679840087891e+0000,
*** 53,63 **** /* P6 */ -1.50565540968422650891511693771462887525558471680e+0000, /* P7 */ -1.81484378878349295050043110677506774663925170898e+0000, /* T0 */ 3.333335997532835641297409611782510896641e-0001, /* T1 */ 2.999997598248363761541668282006867229939e+00, }; ! /* INDENT ON */ #define one q[0] #define P0 q[1] #define P1 q[2] #define P2 q[3] --- 55,65 ---- /* P6 */ -1.50565540968422650891511693771462887525558471680e+0000, /* P7 */ -1.81484378878349295050043110677506774663925170898e+0000, /* T0 */ 3.333335997532835641297409611782510896641e-0001, /* T1 */ 2.999997598248363761541668282006867229939e+00, }; ! #define one q[0] #define P0 q[1] #define P1 q[2] #define P2 q[3]
*** 68,96 **** #define P7 q[8] #define T0 q[9] #define T1 q[10] float ! __k_tanf(double x, int n) { float ft = 0.0; double z, w; int ix; ! ix = ((int *) &x)[HIWORD] & ~0x80000000; /* ix = leading |x| */ /* small argument */ if (ix < 0x3f800000) { /* if |x| < 0.0078125 = 2**-7 */ if (ix < 0x3f100000) { /* if |x| < 2**-14 */ ! if ((int) x == 0) { /* raise inexact if x != 0 */ ! ft = n == 0 ? (float) x : (float) (-one / x); ! } return (ft); } z = (x * T0) * (T1 + x * x); ! ft = n == 0 ? (float) z : (float) (-one / z); return (ft); } z = x * x; ! w = ((P0 * x) * (P1 + z * (P2 + z)) * (P3 + z * (P4 + z))) ! * (P5 + z * (P6 + z * (P7 + z))); ! ft = n == 0 ? (float) w : (float) (-one / w); return (ft); } --- 70,102 ---- #define P7 q[8] #define T0 q[9] #define T1 q[10] float ! __k_tanf(double x, int n) ! { float ft = 0.0; double z, w; int ix; ! ix = ((int *)&x)[HIWORD] & ~0x80000000; /* ix = leading |x| */ ! /* small argument */ if (ix < 0x3f800000) { /* if |x| < 0.0078125 = 2**-7 */ if (ix < 0x3f100000) { /* if |x| < 2**-14 */ ! if ((int)x == 0) /* raise inexact if x != 0 */ ! ft = n == 0 ? (float)x : (float)(-one / x); ! return (ft); } + z = (x * T0) * (T1 + x * x); ! ft = n == 0 ? (float)z : (float)(-one / z); return (ft); } + z = x * x; ! w = ((P0 * x) * (P1 + z * (P2 + z)) * (P3 + z * (P4 + z))) * (P5 + z * ! (P6 + z * (P7 + z))); ! ft = n == 0 ? (float)w : (float)(-one / w); return (ft); }