X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Ffma.c;h=87d450c731ac55d95bd95cbb032c222ec8139f2e;hb=dbdec9722ee604969dddf51e81f8c9e4a9ce064a;hp=a9de2cb37e0f70d73312edeadec8fbd6a21c9992;hpb=b1cbd70743f0e0e0295e92c3e38d7599114db8c6;p=musl diff --git a/src/math/fma.c b/src/math/fma.c index a9de2cb3..87d450c7 100644 --- a/src/math/fma.c +++ b/src/math/fma.c @@ -23,10 +23,6 @@ static void add(long double *hi, long double *lo, long double x, long double y) *lo = y - r; } -/* -TODO(nsz): probably simpler mul is enough if we assume x and y are doubles -so last 11bits are all zeros, no subnormals etc -*/ /* exact mul, assumes no over/underflow */ static void mul(long double *hi, long double *lo, long double x, long double y) { @@ -115,7 +111,7 @@ double fma(double x, double y, double z) add(&hi, &lo2, xy, z); if (hi == 0) { fesetround(round); - /* TODO: verify that the sign of 0 is always correct */ + /* make sure that the sign of 0 is correct */ return (xy + z) + lo1; } } else {