diff options
| author | Szabolcs Nagy <nsz@port70.net> | 2012-11-13 10:55:35 +0100 |
|---|---|---|
| committer | Szabolcs Nagy <nsz@port70.net> | 2012-11-13 10:55:35 +0100 |
| commit | c4359e01303da2755fe7e8033826b132eb3659b1 (patch) | |
| tree | 93cdddff4278c4f0b092ed001f235051a1cd1a72 /src/math/modf.c | |
| parent | 666271c105e4137bdfa195e217799d74143370d4 (diff) | |
| download | musl-c4359e01303da2755fe7e8033826b132eb3659b1.tar.gz musl-c4359e01303da2755fe7e8033826b132eb3659b1.tar.xz musl-c4359e01303da2755fe7e8033826b132eb3659b1.zip | |
math: excess precision fix modf, modff, scalbn, scalbnf
old code was correct only if the result was stored (without the excess precision) or musl was compiled with -ffloat-store. now we use STRICT_ASSIGN to work around the issue. (see note 160 in c11 section 6.8.6.4)
Diffstat (limited to 'src/math/modf.c')
| -rw-r--r-- | src/math/modf.c | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/math/modf.c b/src/math/modf.c index de45069f..8f337ef0 100644 --- a/src/math/modf.c +++ b/src/math/modf.c @@ -1,5 +1,4 @@ -#include <math.h> -#include <stdint.h> +#include "libm.h" double modf(double x, double *iptr) { @@ -33,5 +32,6 @@ double modf(double x, double *iptr) } u.n &= ~mask; *iptr = u.x; - return x - *iptr; + STRICT_ASSIGN(double, x, x - *iptr); + return x; } |