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

34

34

}

35

35

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

36

36

#include <fenv.h>

37

#if LDBL_MANT_DIG == 64

38

#define LASTBIT(u) (u.i.m & 1)

39

#define SPLIT (0x1p32L + 1)

40

#elif LDBL_MANT_DIG == 113

41

#define LASTBIT(u) (u.i.lo & 1)

42

#define SPLIT (0x1p57L + 1)

43

#endif

37

44

38

45

/*

39

46

* A struct dd represents a floating-point number with twice the precision

75

82

static inline long double add_adjusted(long double a, long double b)

76

83

{

77

84

struct dd sum;

78

union IEEEl2bits u;

85

union ldshape u;

79

86

80

87

sum = dd_add(a, b);

81

88

if (sum.lo != 0) {

82

u.e = sum.hi;

83

if ((u.bits.manl & 1) == 0)

89

u.f = sum.hi;

90

if (!LASTBIT(u))

84

91

sum.hi = nextafterl(sum.hi, INFINITY * sum.lo);

85

92

}

86

93

return (sum.hi);

95

102

{

96

103

struct dd sum;

97

104

int bits_lost;

98

union IEEEl2bits u;

105

union ldshape u;

99

106

100

107

sum = dd_add(a, b);

101

108

110

117

* break the ties manually.

111

118

*/

112

119

if (sum.lo != 0) {

113

u.e = sum.hi;

114

bits_lost = -u.bits.exp - scale + 1;

115

if (bits_lost != 1 ^ (int)(u.bits.manl & 1))

120

u.f = sum.hi;

121

bits_lost = -u.i.se - scale + 1;

122

if ((bits_lost != 1) ^ LASTBIT(u))

116

123

sum.hi = nextafterl(sum.hi, INFINITY * sum.lo);

117

124

}

118

125

return scalbnl(sum.hi, scale);

125

132

*/

126

133

static inline struct dd dd_mul(long double a, long double b)

127

134

{

128

#if LDBL_MANT_DIG == 64

129

static const long double split = 0x1p32L + 1.0;

130

#elif LDBL_MANT_DIG == 113

131

static const long double split = 0x1p57L + 1.0;

132

#endif

133

135

struct dd ret;

134

136

long double ha, hb, la, lb, p, q;

135

137

136

p = a * split;

138

p = a * SPLIT;

137

139

ha = a - p;

138

140

ha += p;

139

141

la = a - ha;

140

142

141

p = b * split;

143

p = b * SPLIT;

142

144

hb = b - p;

143

145

hb += p;

144

146

lb = b - hb;

Older »