e_pow.c (7525d42f705e49e74670dceacc09b3ee5388995c) e_pow.c (6813d08ff55ae587abd7e2297e051d491c218de0)
1/* @(#)e_pow.c 1.5 04/04/22 SMI */
2/*
3 * ====================================================
4 * Copyright (C) 2004 by Sun Microsystems, Inc. All rights reserved.
5 *
6 * Permission to use, copy, modify, and distribute this
7 * software is freely granted, provided that this notice
8 * is preserved.

--- 43 unchanged lines hidden (view full) ---

52 *
53 * Constants :
54 * The hexadecimal values are the intended ones for the following
55 * constants. The decimal values may be used, provided that the
56 * compiler will convert from decimal to binary accurately enough
57 * to produce the hexadecimal values shown.
58 */
59
1/* @(#)e_pow.c 1.5 04/04/22 SMI */
2/*
3 * ====================================================
4 * Copyright (C) 2004 by Sun Microsystems, Inc. All rights reserved.
5 *
6 * Permission to use, copy, modify, and distribute this
7 * software is freely granted, provided that this notice
8 * is preserved.

--- 43 unchanged lines hidden (view full) ---

52 *
53 * Constants :
54 * The hexadecimal values are the intended ones for the following
55 * constants. The decimal values may be used, provided that the
56 * compiler will convert from decimal to binary accurately enough
57 * to produce the hexadecimal values shown.
58 */
59
60#include <float.h>
60#include "math.h"
61#include "math_private.h"
62
63static const double
64bp[] = {1.0, 1.5,},
65dp_h[] = { 0.0, 5.84962487220764160156e-01,}, /* 0x3FE2B803, 0x40000000 */
66dp_l[] = { 0.0, 1.35003920212974897128e-08,}, /* 0x3E4CFDEB, 0x43CFD006 */
67zero = 0.0,

--- 234 unchanged lines hidden (view full) ---

302 r = (z*t1)/(t1-two)-(w+z*w);
303 z = one-(r-z);
304 GET_HIGH_WORD(j,z);
305 j += (n<<20);
306 if((j>>20)<=0) z = scalbn(z,n); /* subnormal output */
307 else SET_HIGH_WORD(z,j);
308 return s*z;
309}
61#include "math.h"
62#include "math_private.h"
63
64static const double
65bp[] = {1.0, 1.5,},
66dp_h[] = { 0.0, 5.84962487220764160156e-01,}, /* 0x3FE2B803, 0x40000000 */
67dp_l[] = { 0.0, 1.35003920212974897128e-08,}, /* 0x3E4CFDEB, 0x43CFD006 */
68zero = 0.0,

--- 234 unchanged lines hidden (view full) ---

303 r = (z*t1)/(t1-two)-(w+z*w);
304 z = one-(r-z);
305 GET_HIGH_WORD(j,z);
306 j += (n<<20);
307 if((j>>20)<=0) z = scalbn(z,n); /* subnormal output */
308 else SET_HIGH_WORD(z,j);
309 return s*z;
310}
311
312#if (LDBL_MANT_DIG == 53)
313__weak_reference(pow, powl);
314#endif