/*
* ====================================================
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
*
* Developed at SunPro, a Sun Microsystems, Inc. business.
* Permission to use, copy, modify, and distribute this
* software is freely granted, provided that this notice
* is preserved.
* ====================================================
*/

#include "softmath_private.h"
#include <inttypes.h>

static const float
    bp[]   = {1.0, 1.5,},
    dp_h[] = { 0.0, 5.84960938e-01,}, /* 0x3f15c000 */
    dp_l[] = { 0.0, 1.56322085e-06,}, /* 0x35d1cfdc */
    Zero[] = {0.0, -0.0,},
    zero   =  0.0,
    one    =  1.0,
    two    =  2.0,
    two24  =  16777216.0, /* 0x4b800000 */
    huge   =  1.0e30,
    tiny   =  1.0e-30,
    /* poly coefs for (3/2)*(log(x)-2s-2/3*s**3 */
    L1     =  6.0000002384e-01, /* 0x3f19999a */
    L2     =  4.2857143283e-01, /* 0x3edb6db7 */
    L3     =  3.3333334327e-01, /* 0x3eaaaaab */
    L4     =  2.7272811532e-01, /* 0x3e8ba305 */
    L5     =  2.3066075146e-01, /* 0x3e6c3255 */
    L6     =  2.0697501302e-01, /* 0x3e53f142 */
    P1     =  1.6666667163e-01, /* 0x3e2aaaab */
    P2     = -2.7777778450e-03, /* 0xbb360b61 */
    P3     =  6.6137559770e-05, /* 0x388ab355 */
    P4     = -1.6533901999e-06, /* 0xb5ddea0e */
    P5     =  4.1381369442e-08, /* 0x3331bb4c */
    lg2    =  6.9314718246e-01, /* 0x3f317218 */
    lg2_h  =  6.93145752e-01, /* 0x3f317200 */
    lg2_l  =  1.42860654e-06, /* 0x35bfbe8c */
    ovt    =  4.2995665694e-08, /* -(128-log2(ovfl+.5ulp)) */
    cp     =  9.6179670095e-01, /* 0x3f76384f =2/(3ln2) */
    cp_h   =  9.6191406250e-01, /* 0x3f764000 =12b cp */
    cp_l   = -1.1736857402e-04, /* 0xb8f623c6 =tail of cp_h */
    ivln2  =  1.4426950216e+00, /* 0x3fb8aa3b =1/ln2 */
    ivln2_h=  1.4426879883e+00, /* 0x3fb8aa00 =16b 1/ln2*/
    ivln2_l=  7.0526075433e-06; /* 0x36eca570 =1/ln2 tail*/

typedef unsigned int u_int32_t;

typedef union {
    float value;
    u_int32_t word;
} ieee_float_shape_type;

#define GET_FLOAT_WORD(i,d) do \
{ \
    ieee_float_shape_type gf_u; \
    gf_u.value = (d); \
    (i) = gf_u.word; \
} while(0)

#define SET_FLOAT_WORD(d,i) do \
{ \
    ieee_float_shape_type gf_u; \
    gf_u.word = (i); \
    (d) = gf_u.value; \
} while(0)