JehanneOS
/
newlib
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

2007-10-17 Jeff Johnston <jjohnstn@redhat.com> 

* libm/mathfp/s_logarithm.c: Fix error introduced by previous
        fix when handling negative input values.  Make function
        consistent with math directory and glibc version such that
        inf and nan values return inf and nan respectively with no
        errno setting.
        * libm/mathfp/sf_logarithm.c: Ditto.
        * libm/math/w_log.c: Set errno to ERANGE when input is 0.0.
        * libm/math/wf_log.c: Ditto.
        * libm/math/w_log10.c: Ditto.
        * libm/math/wf_log10.c: Ditto.
This commit is contained in:
Jeff Johnston 2007-10-18 00:03:32 +00:00
parent 923f9573a3
commit 70317d8506
7 changed files with 45 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
newlib/ChangeLog
View File

@ -1,3 +1,16 @@
2007-10-17 Jeff Johnston <jjohnstn@redhat.com>
* libm/mathfp/s_logarithm.c: Fix error introduced by previous
fix when handling negative input values. Make function
consistent with math directory and glibc version such that
inf and nan values return inf and nan respectively with no
errno setting.
* libm/mathfp/sf_logarithm.c: Ditto.
* libm/math/w_log.c: Set errno to ERANGE when input is 0.0.
* libm/math/wf_log.c: Ditto.
* libm/math/w_log10.c: Ditto.
* libm/math/wf_log10.c: Ditto.
2007-10-17 Jeff Johnston <jjohnstn@redhat.com>
* libm/mathfp/s_logarithm.c: Fix case where input is 0 to

2
newlib/libm/math/w_log.c
View File

@ -95,7 +95,7 @@ PORTABILITY
if (_LIB_VERSION == _POSIX_)
errno = ERANGE;
else if (!matherr(&exc)) {
errno = EDOM;
errno = ERANGE;
}
} else {
/* log(x<0) */

2
newlib/libm/math/w_log10.c
View File

@ -93,7 +93,7 @@ PORTABILITY
if (_LIB_VERSION == _POSIX_)
errno = ERANGE;
else if (!matherr(&exc)) {
errno = EDOM;
errno = ERANGE;
}
} else {
/* log10(x<0) */

2
newlib/libm/math/wf_log.c
View File

@ -53,7 +53,7 @@
if (_LIB_VERSION == _POSIX_)
errno = ERANGE;
else if (!matherr(&exc)) {
errno = EDOM;
errno = ERANGE;
}
} else {
/* logf(x<0) */

2
newlib/libm/math/wf_log10.c
View File

@ -54,7 +54,7 @@
if (_LIB_VERSION == _POSIX_)
errno = ERANGE;
else if (!matherr(&exc)) {
errno = EDOM;
errno = ERANGE;
}
} else {
/* log10f(x<0) */

16
newlib/libm/mathfp/s_logarithm.c
View File

@ -100,12 +100,24 @@ _DEFUN (logarithm, (double, int),
int N;
double f, w, z;
/* Check for domain error here. */
if (x <= 0.0)
/* Check for range and domain errors here. */
if (x == 0.0)
{
errno = ERANGE;
return (-z_infinity.d);
}
else if (x < 0.0)
{
errno = EDOM;
return (z_notanum.d);
}
else if (!isfinite(x))
{
if (isnan(x))
return (z_notanum.d);
else
return (z_infinity.d);
}
/* Get the exponent and mantissa where x = f * 2^N. */
f = frexp (x, &N);

16
newlib/libm/mathfp/sf_logarithm.c
View File

@ -38,12 +38,24 @@ _DEFUN (logarithmf, (float, int),
int N;
float f, w, z;
/* Check for domain error here. */
if (x <= 0.0)
/* Check for domain/range errors here. */
if (x == 0.0)
{
errno = ERANGE;
return (-z_infinity_f.f);
}
else if (x < 0.0)
{
errno = EDOM;
return (z_notanum_f.f);
}
else if (!isfinitef(x))
{
if (isnanf(x))
return (z_notanum_f.f);
else
return (z_infinity_f.f);
}
/* Get the exponent and mantissa where x = f * 2^N. */
f = frexpf (x, &N);