← 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/exp2l.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/exp2l.c

33

33

return exp2(x);

34

34

}

35

35

#elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384

36

37

36

#define TBLBITS 7

38

37

#define TBLSIZE (1 << TBLBITS)

39

38

40

#define BIAS (LDBL_MAX_EXP - 1)

41

#define EXPMASK (BIAS + LDBL_MAX_EXP)

42

43

39

static const double

44

40

redux = 0x1.8p63 / TBLSIZE,

45

41

P1 = 0x1.62e42fefa39efp-1,

203

199

*/

204

200

long double exp2l(long double x)

205

201

{

206

union IEEEl2bits u, v;

202

union ldshape u = {x};

203

int e = u.i.se & 0x7fff;

207

204

long double r, z;

208

uint32_t hx, ix, i0;

205

uint32_t i0;

209

206

union {uint32_t u; int32_t i;} k;

210

207

211

208

/* Filter out exceptional cases. */

212

u.e = x;

213

hx = u.xbits.expsign;

214

ix = hx & EXPMASK;

215

if (ix >= BIAS + 14) { /* |x| >= 16384 or x is NaN */

216

if (ix == EXPMASK) {

217

if (u.xbits.man == 1ULL << 63 && hx == 0xffff) /* -inf */

209

if (e >= 0x3fff + 13) { /* |x| >= 8192 or x is NaN */

210

if (u.i.se >= 0x3fff + 14 && u.i.se >> 15 == 0)

211

/* overflow */

212

return x * 0x1p16383L;

213

if (e == 0x7fff) /* -inf or -nan */

214

return -1/x;

215

if (x < -16382) {

216

if (x <= -16446 || x - 0x1p63 + 0x1p63 != x)

217

/* underflow */

218

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

219

if (x <= -16446)

218

220

return 0;

219

return x;

220

}

221

if (x >= 16384) {

222

x *= 0x1p16383L;

223

return x;

224

}

225

if (x <= -16446)

226

return 0x1p-10000L*0x1p-10000L;

227

} else if (ix < BIAS - 64) /* |x| < 0x1p-64 */

221

}

222

} else if (e < 0x3fff - 64) {

228

223

return 1 + x;

224

}

229

225

230

226

/*

231

227

* Reduce x, computing z, i0, and k. The low bits of x + redux

238

234

* We split this into k = 0xabc and i0 = 0x12 (adjusted to

239

235

* index into the table), then we compute z = 0x0.003456p0.

240

236

*/

241

u.e = x + redux;

242

i0 = u.bits.manl + TBLSIZE / 2;

237

u.f = x + redux;

238

i0 = u.i.m + TBLSIZE / 2;

243

239

k.u = i0 / TBLSIZE * TBLSIZE;

244

240

k.i /= TBLSIZE;

245

241

i0 %= TBLSIZE;

246

u.e -= redux;

247

z = x - u.e;

242

u.f -= redux;

243

z = x - u.f;

248

244

249

245

/* Compute r = exp2l(y) = exp2lt[i0] * p(z). */

250

246

long double t_hi = tbl[2*i0];

Older »