← Back to branch summary

~ubuntu-branches/ubuntu/wily/musl/wily

~ubuntu-branches/ubuntu/wily/musl/wily

« back to all changes in this revision

Viewing changes to src/math/exp2f.c

Committer: Package Import Robot
Author(s): Kevin Bortis
Date: 2013-09-27 23:47:18 UTC
mfrom: (1.1.1)
Revision ID: package-import@ubuntu.com-20130927234718-a96bgcnvzx5buf60

Tags: 0.9.14-1

http://bugs.debian.org/724277

http://bugs.debian.org/724247

* Import upstream version 0.9.14
* Only build on fully supported architectures
* Point to new homepage in control file (Closes: #724277)
* Revorked debian/rules
* Solved possible problem with postrm script (Closes: #724247)

files added:
include/wait.h

files removed:
debian/musl.postinst

debian/musl.postrm

src/internal/longdbl.h

files modified:
Makefile

WHATSNEW

arch/arm/atomic.h

arch/microblaze/bits/limits.h

arch/microblaze/reloc.h

arch/mips/bits/limits.h

arch/mips/bits/msg.h

arch/mips/bits/shm.h

arch/mips/reloc.h

arch/powerpc/bits/limits.h

arch/powerpc/reloc.h

debian/changelog

debian/control

debian/copyright

debian/musl-dev.install

debian/musl-dev.links

debian/musl-tools.install

debian/musl.install

debian/musl.links

debian/rules

debian/scripts/config.path.in

debian/scripts/ld-musl-config.in

debian/source/lintian-overrides

include/elf.h

include/limits.h

include/net/if_arp.h

include/netinet/in.h

include/sys/ptrace.h

include/sys/resource.h

include/sys/socket.h

src/aio/aio_readwrite.c

src/aio/lio_listio.c

src/conf/fpathconf.c

src/conf/sysconf.c

src/env/__libc_start_main.c

src/internal/libc.h

src/internal/libm.h

src/ldso/dynlink.c

src/legacy/getpagesize.c

src/legacy/valloc.c

src/malloc/malloc.c

src/math/__cosl.c

src/math/__fpclassify.c

src/math/__fpclassifyf.c

src/math/__fpclassifyl.c

src/math/__invtrigl.c

src/math/__invtrigl.h

src/math/__rem_pio2.c

src/math/__rem_pio2_large.c

src/math/__rem_pio2f.c

src/math/__rem_pio2l.c

src/math/__signbitl.c

src/math/__sinl.c

src/math/__tanl.c

src/math/acoshf.c

src/math/acoshl.c

src/math/acosl.c

src/math/asinhl.c

src/math/asinl.c

src/math/atan2l.c

src/math/atanh.c

src/math/atanhf.c

src/math/atanhl.c

src/math/atanl.c

src/math/cbrt.c

src/math/cbrtf.c

src/math/cbrtl.c

src/math/ceil.c

src/math/ceilf.c

src/math/ceill.c

src/math/copysign.c

src/math/copysignf.c

src/math/copysignl.c

src/math/coshl.c

src/math/cosl.c

src/math/erfl.c

src/math/exp.c

src/math/exp2.c

src/math/exp2f.c

src/math/exp2l.c

src/math/expf.c

src/math/expl.c

src/math/expm1.c

src/math/expm1f.c

src/math/fabs.c

src/math/fabsf.c

src/math/fabsl.c

src/math/floor.c

src/math/floorf.c

src/math/floorl.c

src/math/fma.c

src/math/fmal.c

src/math/fmod.c

src/math/fmodf.c

src/math/fmodl.c

src/math/frexpl.c

src/math/hypot.c

src/math/hypotf.c

src/math/hypotl.c

src/math/i386/exp.s

src/math/ilogb.c

src/math/ilogbf.c

src/math/ilogbl.c

src/math/lgamma_r.c

src/math/lgammaf_r.c

src/math/lgammal.c

src/math/log1p.c

src/math/log1pf.c

src/math/logb.c

src/math/logbf.c

src/math/logbl.c

src/math/modf.c

src/math/modff.c

src/math/modfl.c

src/math/nextafter.c

src/math/nextafterf.c

src/math/nextafterl.c

src/math/nexttoward.c

src/math/nexttowardf.c

src/math/remainder.c

src/math/remainderf.c

src/math/remquo.c

src/math/remquof.c

src/math/remquol.c

src/math/rint.c

src/math/rintf.c

src/math/rintl.c

src/math/round.c

src/math/roundf.c

src/math/roundl.c

src/math/scalbn.c

src/math/scalbnf.c

src/math/scalbnl.c

src/math/signgam.c

src/math/sincosl.c

src/math/sinhl.c

src/math/sinl.c

src/math/sqrtl.c

src/math/tanhl.c

src/math/tanl.c

src/math/trunc.c

src/math/truncf.c

src/math/truncl.c

src/math/x86_64/ceill.s

src/math/x86_64/exp2l.s

src/math/x86_64/lrintl.s

src/math/x86_64/truncl.s

src/misc/realpath.c

src/mman/mprotect.c

src/signal/sigemptyset.c

src/signal/sigfillset.c

src/signal/sigqueue.c

src/stdio/popen.c

src/stdio/vfscanf.c

src/stdio/vfwscanf.c

src/temp/__randname.c

src/thread/__lock.c

src/thread/mips/clone.s

src/thread/pthread_create.c

src/thread/pthread_getattr_np.c

src/thread/synccall.c

src/unistd/faccessat.c

Show diffs side-by-side

added added

removed removed

src/math/exp2f.c

63

63

* Method: (equally-spaced tables)

64

64

*

65

65

* Reduce x:

66

* x = 2**k + y, for integer k and |y| <= 1/2.

66

* x = k + y, for integer k and |y| <= 1/2.

67

67

* Thus we have exp2f(x) = 2**k * exp2(y).

68

68

*

69

69

* Reduce y:

83

83

*/

84

84

float exp2f(float x)

85

85

{

86

double tv, twopk, u, z;

87

float t;

88

uint32_t hx, ix, i0, k;

86

double_t t, r, z;

87

union {float f; uint32_t i;} u = {x};

88

union {double f; uint64_t i;} uk;

89

uint32_t ix, i0, k;

89

90

90

91

/* Filter out exceptional cases. */

91

GET_FLOAT_WORD(hx, x);

92

ix = hx & 0x7fffffff;

93

if (ix >= 0x43000000) { /* |x| >= 128 */

94

if (ix >= 0x7f800000) {

95

if (hx == 0xff800000) /* -inf */

92

ix = u.i & 0x7fffffff;

93

if (ix > 0x42fc0000) { /* |x| > 126 */

94

if (u.i >= 0x43000000 && u.i < 0x80000000) { /* x >= 128 */

95

x *= 0x1p127f;

96

return x;

97

}

98

if (u.i >= 0x80000000) { /* x < -126 */

99

if (u.i >= 0xc3160000 || (u.i & 0x0000ffff))

100

FORCE_EVAL(-0x1p-149f/x);

101

if (u.i >= 0xc3160000) /* x <= -150 */

96

102

return 0;

97

return x;

98

}

99

if (x >= 128) {

100

STRICT_ASSIGN(float, x, x * 0x1p127f);

101

return x;

102

}

103

if (x <= -150) {

104

STRICT_ASSIGN(float, x, 0x1p-100f*0x1p-100f);

105

return x;

106

103

}

107

104

} else if (ix <= 0x33000000) { /* |x| <= 0x1p-25 */

108

105

return 1.0f + x;

109

106

}

110

107

111

108

/* Reduce x, computing z, i0, and k. */

112

STRICT_ASSIGN(float, t, x + redux);

113

GET_FLOAT_WORD(i0, t);

109

u.f = x + redux;

110

i0 = u.i;

114

111

i0 += TBLSIZE / 2;

115

k = (i0 / TBLSIZE) << 20;

112

k = i0 / TBLSIZE;

113

uk.i = (uint64_t)(0x3ff + k)<<52;

116

114

i0 &= TBLSIZE - 1;

117

t -= redux;

118

z = x - t;

119

INSERT_WORDS(twopk, 0x3ff00000 + k, 0);

120

115

u.f -= redux;

116

z = x - u.f;

121

117

/* Compute r = exp2(y) = exp2ft[i0] * p(z). */

122

tv = exp2ft[i0];

123

u = tv * z;

124

tv = tv + u * (P1 + z * P2) + u * (z * z) * (P3 + z * P4);

118

r = exp2ft[i0];

119

t = r * z;

120

r = r + t * (P1 + z * P2) + t * (z * z) * (P3 + z * P4);

125

121

126

/* Scale by 2**(k>>20). */

127

return tv * twopk;

122

/* Scale by 2**k */

123

return r * uk.f;

128

124

}

Older »