jyapayne/musl
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

math: raise exception flags in log1pl.c on <= -1 arguments

Browse source
This commit is contained in:
Szabolcs Nagy 2012-11-13 00:21:09 +01:00
commit 04ccbdca6d
1 changed files with 2 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

9
src/math/log1pl.c
View file

@ -46,11 +46,6 @@
* Relative error: * Relative error:
* arithmetic domain # trials peak rms * arithmetic domain # trials peak rms
* IEEE -1.0, 9.0 100000 8.2e-20 2.5e-20 * IEEE -1.0, 9.0 100000 8.2e-20 2.5e-20
*
* ERROR MESSAGES:
*
* log singularity: x-1 = 0; returns -INFINITY
* log domain: x-1 < 0; returns NAN
*/ */
#include "libm.h" #include "libm.h"
@ -123,8 +118,8 @@ long double log1pl(long double xm1)
/* Test for domain errors. */ /* Test for domain errors. */
if (x <= 0.0) { if (x <= 0.0) {
if (x == 0.0) if (x == 0.0)
return -INFINITY; return -1/x; /* -inf with divbyzero */
return NAN; return 0/0.0f; /* nan with invalid */
} }
/* Separate mantissa from exponent. /* Separate mantissa from exponent.