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

*** 16,45 **** * fields enclosed by brackets "[]" replaced with your own identifying * information: Portions Copyright [yyyy] [name of copyright owner] * * CDDL HEADER END */ /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #pragma weak __exp2 = exp2 ! /* INDENT OFF */ /* * exp2(x) * Code by K.C. Ng for SUN 4.0 libm. * Method : * exp2(x) = 2**x = 2**((x-anint(x))+anint(x)) * = 2**anint(x)*2**(x-anint(x)) * = 2**anint(x)*exp((x-anint(x))*ln2) */ - /* INDENT ON */ #include "libm.h" static const double C[] = { 0.0, --- 16,46 ---- * fields enclosed by brackets "[]" replaced with your own identifying * information: Portions Copyright [yyyy] [name of copyright owner] * * CDDL HEADER END */ + /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ + /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #pragma weak __exp2 = exp2 ! /* * exp2(x) * Code by K.C. Ng for SUN 4.0 libm. * Method : * exp2(x) = 2**x = 2**((x-anint(x))+anint(x)) * = 2**anint(x)*2**(x-anint(x)) * = 2**anint(x)*exp((x-anint(x))*ln2) */ #include "libm.h" static const double C[] = { 0.0,
*** 56,88 **** #define ln2 C[3] #define huge C[4] #define tiny C[5] double ! exp2(double x) { int ix, hx, k; double t; ix = ((int *)&x)[HIWORD]; hx = ix & ~0x80000000; if (hx >= 0x4090e000) { /* |x| >= 1080 or x is nan */ if (hx >= 0x7ff00000) { /* x is inf or nan */ if (ix == 0xfff00000 && ((int *)&x)[LOWORD] == 0) return (zero); return (x * x); } ! t = (ix < 0)? tiny : huge; return (t * t); } if (hx < 0x3fe00000) { /* |x| < 0.5 */ if (hx < 0x3c000000) return (one + x); return (exp(ln2 * x)); } k = (int)x; if (x != (double)k) ! k = (int)((ix < 0)? x - half : x + half); return (scalbn(exp(ln2 * (x - (double)k)), k)); } --- 57,95 ---- #define ln2 C[3] #define huge C[4] #define tiny C[5] double ! exp2(double x) ! { int ix, hx, k; double t; ix = ((int *)&x)[HIWORD]; hx = ix & ~0x80000000; if (hx >= 0x4090e000) { /* |x| >= 1080 or x is nan */ if (hx >= 0x7ff00000) { /* x is inf or nan */ if (ix == 0xfff00000 && ((int *)&x)[LOWORD] == 0) return (zero); + return (x * x); } ! ! t = (ix < 0) ? tiny : huge; return (t * t); } if (hx < 0x3fe00000) { /* |x| < 0.5 */ if (hx < 0x3c000000) return (one + x); + return (exp(ln2 * x)); } k = (int)x; + if (x != (double)k) ! k = (int)((ix < 0) ? x - half : x + half); ! return (scalbn(exp(ln2 * (x - (double)k)), k)); }