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

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/hypotl.c

/* origin: FreeBSD /usr/src/lib/msun/src/e_hypotl.c */

* ====================================================

* Developed at SunSoft, a Sun Microsystems, Inc. business.

* Permission to use, copy, modify, and distribute this

* software is freely granted, provided that this notice

* is preserved.

* ====================================================

/* long double version of hypot(). See comments in hypot.c. */

#include "libm.h"

#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024

return hypot(x, y);

}

#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384

#define GET_LDBL_EXPSIGN(i, v) do { \

union IEEEl2bits uv; \

uv.e = v; \

i = uv.xbits.expsign; \

} while (0)

#define GET_LDBL_MAN(h, l, v) do { \

union IEEEl2bits uv; \

uv.e = v; \

h = uv.bits.manh; \

l = uv.bits.manl; \

} while (0)

#define SET_LDBL_EXPSIGN(v, i) do { \

union IEEEl2bits uv; \

uv.e = v; \

uv.xbits.expsign = i; \

v = uv.e; \

} while (0)

#undef GET_HIGH_WORD

#define GET_HIGH_WORD(i, v) GET_LDBL_EXPSIGN(i, v)

#undef SET_HIGH_WORD

#define SET_HIGH_WORD(v, i) SET_LDBL_EXPSIGN(v, i)

#define DESW(exp) (exp) /* delta expsign word */

#define ESW(exp) (MAX_EXP - 1 + (exp)) /* expsign word */

#define MANT_DIG LDBL_MANT_DIG

#define MAX_EXP LDBL_MAX_EXP

#if LDBL_MANL_SIZE > 32

typedef uint64_t man_t;

#else

typedef uint32_t man_t;

#if LDBL_MANT_DIG == 64

#define SPLIT (0x1p32L+1)

#elif LDBL_MANT_DIG == 113

#define SPLIT (0x1p57L+1)

#endif

static void sq(long double *hi, long double *lo, long double x)

{

long double xh, xl, xc;

xc = x*SPLIT;

xh = x - xc + xc;

xl = x - xh;

*hi = x*x;

*lo = xh*xh - *hi + 2*xh*xl + xl*xl;

}

long double hypotl(long double x, long double y)

{

long double a=x,b=y,t1,t2,y1,y2,w;

int32_t j,k,ha,hb;

GET_HIGH_WORD(ha, x);

ha &= 0x7fff;

GET_HIGH_WORD(hb, y);

hb &= 0x7fff;

if (hb > ha) {

a = y;

b = x;

j=ha; ha=hb; hb=j;

} else {

a = x;

b = y;

}

a = fabsl(a);

b = fabsl(b);

if (ha - hb > DESW(MANT_DIG+7)) /* x/y > 2**(MANT_DIG+7) */

return a+b;

k = 0;

if (ha > ESW(MAX_EXP/2-12)) { /* a>2**(MAX_EXP/2-12) */

if (ha >= ESW(MAX_EXP)) { /* Inf or NaN */

man_t manh, manl;

/* Use original arg order iff result is NaN; quieten sNaNs. */

w = fabsl(x+0.0)-fabsl(y+0.0);

GET_LDBL_MAN(manh,manl,a);

if (manh == LDBL_NBIT && manl == 0) w = a;

GET_LDBL_MAN(manh,manl,b);

if (hb >= ESW(MAX_EXP) && manh == LDBL_NBIT && manl == 0) w = b;

return w;

}

/* scale a and b by 2**-(MAX_EXP/2+88) */

ha -= DESW(MAX_EXP/2+88); hb -= DESW(MAX_EXP/2+88);

k += MAX_EXP/2+88;

SET_HIGH_WORD(a, ha);

SET_HIGH_WORD(b, hb);

100

}

101

if (hb < ESW(-(MAX_EXP/2-12))) { /* b < 2**-(MAX_EXP/2-12) */

102

if (hb <= 0) { /* subnormal b or 0 */

103

man_t manh, manl;

104

GET_LDBL_MAN(manh,manl,b);

105

if ((manh|manl) == 0)

106

return a;

107

t1 = 0;

108

SET_HIGH_WORD(t1, ESW(MAX_EXP-2)); /* t1 = 2^(MAX_EXP-2) */

109

b *= t1;

110

a *= t1;

111

k -= MAX_EXP-2;

112

} else { /* scale a and b by 2^(MAX_EXP/2+88) */

113

ha += DESW(MAX_EXP/2+88);

114

hb += DESW(MAX_EXP/2+88);

115

k -= MAX_EXP/2+88;

116

SET_HIGH_WORD(a, ha);

117

SET_HIGH_WORD(b, hb);

118

}

119

}

120

/* medium size a and b */

121

w = a - b;

122

if (w > b) {

123

t1 = a;

124

union IEEEl2bits uv;

125

uv.e = t1; uv.bits.manl = 0; t1 = uv.e;

126

t2 = a-t1;

127

w = sqrtl(t1*t1-(b*(-b)-t2*(a+t1)));

128

} else {

129

a = a+a;

130

y1 = b;

131

union IEEEl2bits uv;

132

uv.e = y1; uv.bits.manl = 0; y1 = uv.e;

133

y2 = b - y1;

134

t1 = a;

135

uv.e = t1; uv.bits.manl = 0; t1 = uv.e;

136

t2 = a - t1;

137

w = sqrtl(t1*y1-(w*(-w)-(t1*y2+t2*b)));

138

}

139

if(k!=0) {

140

uint32_t high;

141

t1 = 1.0;

142

GET_HIGH_WORD(high, t1);

143

SET_HIGH_WORD(t1, high+DESW(k));

144

return t1*w;

145

}

146

return w;

union ldshape ux = {x}, uy = {y};

int ex, ey;

long double hx, lx, hy, ly, z;

ux.i.se &= 0x7fff;

uy.i.se &= 0x7fff;

if (ux.i.se < uy.i.se) {

ex = uy.i.se;

ey = ux.i.se;

x = uy.f;

y = ux.f;

} else {

ex = ux.i.se;

ey = uy.i.se;

x = ux.f;

y = uy.f;

}

if (ex == 0x7fff && isinf(y))

return y;

if (ex == 0x7fff || y == 0)

return x;

if (ex - ey > LDBL_MANT_DIG)

return x + y;

z = 1;

if (ex > 0x3fff+8000) {

z = 0x1p10000L;

x *= 0x1p-10000L;

y *= 0x1p-10000L;

} else if (ey < 0x3fff-8000) {

z = 0x1p-10000L;

x *= 0x1p10000L;

y *= 0x1p10000L;

}

sq(&hx, &lx, x);

sq(&hy, &ly, y);

return z*sqrtl(ly+lx+hy+hx);

147

}

148

#endif

Older »