← 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/exp.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

80

80

81

81

double exp(double x)

82

82

{

83

double hi, lo, c, xx;

83

double_t hi, lo, c, xx, y;

84

84

int k, sign;

85

85

uint32_t hx;

86

86

89

89

hx &= 0x7fffffff; /* high word of |x| */

90

90

91

91

/* special cases */

92

if (hx >= 0x40862e42) { /* if |x| >= 709.78... */

92

if (hx >= 0x4086232b) { /* if |x| >= 708.39... */

93

93

if (isnan(x))

94

94

return x;

95

if (hx == 0x7ff00000 && sign) /* -inf */

96

return 0;

97

95

if (x > 709.782712893383973096) {

98

96

/* overflow if x!=inf */

99

STRICT_ASSIGN(double, x, 0x1p1023 * x);

97

x *= 0x1p1023;

100

98

return x;

101

99

}

102

if (x < -745.13321910194110842) {

103

/* underflow */

104

STRICT_ASSIGN(double, x, 0x1p-1000 * 0x1p-1000);

105

return x;

100

if (x < -708.39641853226410622) {

101

/* underflow if x!=-inf */

102

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

103

if (x < -745.13321910194110842)

104

return 0;

106

105

}

107

106

}

108

107

114

113

k = 1 - sign - sign;

115

114

hi = x - k*ln2hi; /* k*ln2hi is exact here */

116

115

lo = k*ln2lo;

117

STRICT_ASSIGN(double, x, hi - lo);

116

x = hi - lo;

118

117

} else if (hx > 0x3e300000) { /* if |x| > 2**-28 */

119

118

k = 0;

120

119

hi = x;

128

127

/* x is now in primary range */

129

128

xx = x*x;

130

129

c = x - xx*(P1+xx*(P2+xx*(P3+xx*(P4+xx*P5))));

131

x = 1 + (x*c/(2-c) - lo + hi);

130

y = 1 + (x*c/(2-c) - lo + hi);

132

131

if (k == 0)

133

return x;

134

return scalbn(x, k);

132

return y;

133

return scalbn(y, k);

135

134

}

Older »