b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishra#pragma weak pow = __pow

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishra#include "libm.h"

b6917abefc343244b784f0cc34bc65b01469c3bfmishra#include "xpg6.h" /* __xpg6 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra#define _C99SUSv3_pow _C99SUSv3_pow_treats_Inf_as_an_even_int

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishrastatic const double zero = 0.0, one = 1.0, two = 2.0;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misraextern const double _TBL_log2_hi[], _TBL_log2_lo[];

b6917abefc343244b784f0cc34bc65b01469c3bfmishrastatic const double

b6917abefc343244b784f0cc34bc65b01469c3bfmishra two53 = 9007199254740992.0,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra A1_hi = 2.8853900432586669921875,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra A1_lo = 3.8519259825035041963606002e-8,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra A1 = 2.885390081777926817222541963606002026086e+0000,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra A2 = 9.617966939207270828380543979852286255862e-0001,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra A3 = 5.770807680887875964868853124873696201995e-0001,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra B0_hi = 2.8853900432586669921875,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra B0_lo = 3.8519259822532793056374320585e-8,

4c1c9391ff5b1b69f6fc02ce905b96759cc50fc0Joe Bonasera B0 = 2.885390081777926814720293056374320585689e+0000,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra B1 = 9.617966939259755138949202350396200257632e-0001,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra B2 = 5.770780163585687000782112776448797953382e-0001,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra B3 = 4.121985488948771523290174512461778354953e-0001,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra B4 = 3.207590534812432970433641789022666850193e-0001;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishrastatic double

b6917abefc343244b784f0cc34bc65b01469c3bfmishralog2_x(double x, double *w) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra double f, s, z, qn, h, t;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra int *px = (int *) &x;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra int *pz = (int *) &z;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra int i, j, ix, n;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishra n = 0;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra ix = px[HIWORD];

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (ix >= 0x3fef03f1 && ix < 0x3ff08208) { /* 65/63 > x > 63/65 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra double f1, v;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra f = x - one;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (((ix - 0x3ff00000) | px[LOWORD]) == 0) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra *w = zero;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (zero); /* log2(1)= +0 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra }

b6917abefc343244b784f0cc34bc65b01469c3bfmishra qn = one / (two + f);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra s = f * qn; /* |s|<2**-6 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra v = s * s;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra h = (double) ((float) s);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra f1 = (double) ((float) f);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra t = qn * (((f - two * h) - h * f1) - h * (f - f1));

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* s = h+t */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra f1 = h * B0_lo + s * (v * (B1 + v * (B2 + v * (B3 + v * B4))));

b6917abefc343244b784f0cc34bc65b01469c3bfmishra t = f1 + t * B0;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra h *= B0_hi;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra s = (double) ((float) (h + t));

b6917abefc343244b784f0cc34bc65b01469c3bfmishra *w = t - (s - h);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (s);

2ef50f010f7a3a07eb5a9f6001b9843fd868e26bJoe Bonasera }

2ef50f010f7a3a07eb5a9f6001b9843fd868e26bJoe Bonasera if (ix < 0x00100000) { /* subnormal x */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra x *= two53;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra n = -53;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra ix = px[HIWORD];

b6917abefc343244b784f0cc34bc65b01469c3bfmishra }

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* LARGE N */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra n += ((ix + 0x1000) >> 20) - 0x3ff;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra ix = (ix & 0x000fffff) | 0x3ff00000; /* scale x to [1,2] */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra px[HIWORD] = ix;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra i = ix + 0x1000;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra pz[HIWORD] = i & 0xffffe000;

2ef50f010f7a3a07eb5a9f6001b9843fd868e26bJoe Bonasera pz[LOWORD] = 0;

2ef50f010f7a3a07eb5a9f6001b9843fd868e26bJoe Bonasera qn = one / (x + z);

2ef50f010f7a3a07eb5a9f6001b9843fd868e26bJoe Bonasera f = x - z;

2ef50f010f7a3a07eb5a9f6001b9843fd868e26bJoe Bonasera s = f * qn;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra h = (double) ((float) s);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra t = qn * ((f - (h + h) * z) - h * f);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra j = (i >> 13) & 0x7f;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra f = s * s;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra t = t * A1 + h * A1_lo;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra t += (s * f) * (A2 + f * A3);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra qn = h * A1_hi;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra s = n + _TBL_log2_hi[j];

b6917abefc343244b784f0cc34bc65b01469c3bfmishra h = qn + s;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra t += _TBL_log2_lo[j] - ((h - s) - qn);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra f = (double) ((float) (h + t));

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra *w = t - (f - h);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (f);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra}

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishraextern const double _TBL_exp2_hi[], _TBL_exp2_lo[];

b6917abefc343244b784f0cc34bc65b01469c3bfmishrastatic const double /* poly app of 2^x-1 on [-1e-10,2^-7+1e-10] */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra E1 = 6.931471805599453100674958533810346197328e-0001,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra E2 = 2.402265069587779347846769151717493815979e-0001,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra E3 = 5.550410866475410512631124892773937864699e-0002,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra E4 = 9.618143209991026824853712740162451423355e-0003,

b6917abefc343244b784f0cc34bc65b01469c3bfmishra E5 = 1.333357676549940345096774122231849082991e-0003;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishrapow(double x, double y) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra double z, ax;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra double y1, y2, w1, w2;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra int sbx, sby, j, k, yisint;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra int hx, hy, ahx, ahy;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra unsigned lx, ly;

f9c480cd662f42d54475b87370bc2f4811ec0679Saurabh Misra int *pz = (int *) &z;

f9c480cd662f42d54475b87370bc2f4811ec0679Saurabh Misra

b6917abefc343244b784f0cc34bc65b01469c3bfmishra hx = ((int *) &x)[HIWORD];

b6917abefc343244b784f0cc34bc65b01469c3bfmishra lx = ((unsigned *) &x)[LOWORD];

b6917abefc343244b784f0cc34bc65b01469c3bfmishra hy = ((int *) &y)[HIWORD];

b6917abefc343244b784f0cc34bc65b01469c3bfmishra ly = ((unsigned *) &y)[LOWORD];

b6917abefc343244b784f0cc34bc65b01469c3bfmishra ahx = hx & ~0x80000000;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra ahy = hy & ~0x80000000;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if ((ahy | ly) == 0) { /* y==zero */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if ((ahx | lx) == 0)

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra z = _SVID_libm_err(x, y, 20); /* +-0**+-0 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else if ((ahx | (((lx | -lx) >> 31) & 1)) > 0x7ff00000)

b6917abefc343244b784f0cc34bc65b01469c3bfmishra z = _SVID_libm_err(x, y, 42); /* NaN**+-0 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else

b6917abefc343244b784f0cc34bc65b01469c3bfmishra z = one; /* x**+-0 = 1 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (z);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra } else if (hx == 0x3ff00000 && lx == 0 &&

b6917abefc343244b784f0cc34bc65b01469c3bfmishra (__xpg6 & _C99SUSv3_pow) != 0)

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (one); /* C99: 1**anything = 1 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else if (ahx > 0x7ff00000 || (ahx == 0x7ff00000 && lx != 0) ||

b6917abefc343244b784f0cc34bc65b01469c3bfmishra ahy > 0x7ff00000 || (ahy == 0x7ff00000 && ly != 0))

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (x * y); /* +-NaN return x*y; + -> * for Cheetah */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* includes Sun: 1**NaN = NaN */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra sbx = (unsigned) hx >> 31;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra sby = (unsigned) hy >> 31;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra ax = fabs(x);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /*

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett yisint = 0;

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett if (sbx) {

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett if (ahy >= 0x43400000)

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett yisint = 2; /* even integer y */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else if (ahy >= 0x3ff00000) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra k = (ahy >> 20) - 0x3ff; /* exponent */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (k > 20) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra j = ly >> (52 - k);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if ((j << (52 - k)) == ly)

b6917abefc343244b784f0cc34bc65b01469c3bfmishra yisint = 2 - (j & 1);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra } else if (ly == 0) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra j = ahy >> (20 - k);

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett if ((j << (20 - k)) == ahy)

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett yisint = 2 - (j & 1);

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett }

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett }

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett }

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett /* special value of y */

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett if (ly == 0) {

d23e508cd51e6b3ec3f10a427259d7dd2592fa94Edward Gillett if (ahy == 0x7ff00000) { /* y is +-inf */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (((ahx - 0x3ff00000) | lx) == 0) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if ((__xpg6 & _C99SUSv3_pow) != 0)

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (one);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* C99: (-1)**+-inf = 1 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (y - y);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* Sun: (+-1)**+-inf = NaN */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra } else if (ahx >= 0x3ff00000)

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* (|x|>1)**+,-inf = inf,0 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (sby == 0 ? y : zero);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else /* (|x|<1)**-,+inf = inf,0 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (sby != 0 ? -y : zero);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra }

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (ahy == 0x3ff00000) { /* y is +-1 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (sby != 0) { /* y is -1 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (x == zero) /* divided by zero */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (_SVID_libm_err(x, y, 23));

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else if (ahx < 0x40000 || ((ahx - 0x40000) |

b6917abefc343244b784f0cc34bc65b01469c3bfmishra lx) == 0) /* overflow */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (_SVID_libm_err(x, y, 21));

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (one / x);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra } else

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (x);

b6917abefc343244b784f0cc34bc65b01469c3bfmishra }

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (hy == 0x40000000) { /* y is 2 */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (ahx >= 0x5ff00000 && ahx < 0x7ff00000)

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (_SVID_libm_err(x, y, 21));

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* x*x overflow */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra else if ((ahx < 0x1e56a09e && (ahx | lx) != 0) ||

b6917abefc343244b784f0cc34bc65b01469c3bfmishra (ahx == 0x1e56a09e && lx < 0x667f3bcd))

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (_SVID_libm_err(x, y, 22));

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* x*x underflow */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra else

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra return (x * x);

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra }

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra if (hy == 0x3fe00000) {

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra if (!((ahx | lx) == 0 || ((ahx - 0x7ff00000) | lx) ==

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra 0 || sbx == 1))

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra return (sqrt(x)); /* y is 0.5 and x > 0 */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra }

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra }

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* special value of x */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (lx == 0) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (ahx == 0x7ff00000 || ahx == 0 || ahx == 0x3ff00000) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /* x is +-0,+-inf,-1 */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra z = ax;

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (sby == 1) {

b6917abefc343244b784f0cc34bc65b01469c3bfmishra z = one / z; /* z = |x|**y */

b6917abefc343244b784f0cc34bc65b01469c3bfmishra if (ahx == 0)

b6917abefc343244b784f0cc34bc65b01469c3bfmishra return (_SVID_libm_err(x, y, 23));

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra }

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra if (sbx == 1) {

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra if (ahx == 0x3ff00000 && yisint == 0)

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra z = _SVID_libm_err(x, y, 24);

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* neg**non-integral is NaN + invalid */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra else if (yisint == 1)

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra z = -z; /* (x<0)**odd = -(|x|**odd) */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra }

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra return (z);

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra }

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra }

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* (x<0)**(non-int) is NaN */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra if (sbx == 1 && yisint == 0)

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra return (_SVID_libm_err(x, y, 24));

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* Now ax is finite, y is finite */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* first compute log2(ax) = w1+w2, with 24 bits w1 */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra w1 = log2_x(ax, &w2);

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* split up y into y1+y2 and compute (y1+y2)*(w1+w2) */

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra if (((ly & 0x07ffffff) == 0) || ahy >= 0x47e00000 ||

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra ahy <= 0x38100000) {

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* no need to split if y is short or too large or too small */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra y1 = y * w1;

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra y2 = y * w2;

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra } else {

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra y1 = (double) ((float) y);

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra y2 = (y - y1) * w1 + y * w2;

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra y1 *= w1;

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra }

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra z = y1 + y2;

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra j = pz[HIWORD];

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra if (j >= 0x40900000) { /* z >= 1024 */

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra if (!(j == 0x40900000 && pz[LOWORD] == 0)) /* z > 1024 */

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra return (_SVID_libm_err(x, y, 21)); /* overflow */

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra else {

e511d54dfc1c7eb3aea1a9125b54791fc2f23d42Saurabh Misra w2 = y1 - z;

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra w2 += y2;

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* rounded to inf */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra if (w2 >= -8.008566259537296567160e-17)

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra return (_SVID_libm_err(x, y, 21));

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra /* overflow */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra }

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra } else if ((j & ~0x80000000) >= 0x4090cc00) { /* z <= -1075 */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra if (!(j == 0xc090cc00 && pz[LOWORD] == 0)) /* z < -1075 */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra return (_SVID_libm_err(x, y, 22)); /* underflow */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra else {

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra w2 = y1 - z;

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra w2 += y2;

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra if (w2 <= zero) /* underflow */

325e77f4fc68ea7eea21f5e705d3cf89fd3a7e3dSaurabh Misra return (_SVID_libm_err(x, y, 22));

b6917abefc343244b784f0cc34bc65b01469c3bfmishra }

b6917abefc343244b784f0cc34bc65b01469c3bfmishra }

b6917abefc343244b784f0cc34bc65b01469c3bfmishra /*

b6917abefc343244b784f0cc34bc65b01469c3bfmishra k = (int) (z * 64.0 + (((hy ^ (ahx - 0x3ff00000)) > 0) ? 0.5 : -0.5));

b6917abefc343244b784f0cc34bc65b01469c3bfmishra j = k & 63;

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra w1 = y2 - ((double) k * 0.015625 - y1);

87cc626978cedbea111888dba5e4042df5cf532cSaurabh Misra w2 = _TBL_exp2_hi[j];

87cc626978cedbea111888dba5e4042df5cf532cSaurabh Misra z = _TBL_exp2_lo[j] + (w2 * w1) * (E1 + w1 * (E2 + w1 * (E3 + w1 *

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra (E4 + w1 * E5))));

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra z += w2;

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra k >>= 6;

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra if (k < -1021)

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra z = scalbn(z, k);

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra else /* subnormal output */

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra pz[HIWORD] += k << 20;

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra if (sbx == 1 && yisint == 1)

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra z = -z; /* (-ve)**(odd int) */

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra return (z);

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra}

5d8efbbc545e36559256557618fcd256b7da72e4Saurabh Misra