Mirrors_Framework/nuttx
1
0
Fork 0
mirror of https://github.com/apache/nuttx.git synced 2026-09-27 18:58:27 +08:00
Code Issues Packages Projects Releases Wiki Activity

libm: Fix/improve range checks in lib_expi() and lib_expif().

Browse Source
This commit is contained in:
Gregory Nutt
2016-07-12 09:59:10 -06:00
parent de3b8d85a4
commit 96edfdb4c7
2 changed files with 10 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files
libc/math/lib_libexpi.c
+5 -5
View File
@@ -52,13 +52,12 @@
#define M_E128 (M_E64 * M_E64) #define M_E128 (M_E64 * M_E64)
#define M_E256 (M_E128 * M_E128) #define M_E256 (M_E128 * M_E128)
#define M_E512 (M_E256 * M_E256) #define M_E512 (M_E256 * M_E256)
#define M_E1024 (M_E512 * M_E512)
/**************************************************************************** /****************************************************************************
* Private Data * Private Data
****************************************************************************/ ****************************************************************************/
static const double g_expi_square_tbl[11] = static const double g_expi_square_tbl[] =
{ {
M_E, /* e^1 */ M_E, /* e^1 */
M_E2, /* e^2 */ M_E2, /* e^2 */
@@ -69,8 +68,7 @@ static const double g_expi_square_tbl[11] =
M_E64, /* e^64 */ M_E64, /* e^64 */
M_E128, /* e^128 */ M_E128, /* e^128 */
M_E256, /* e^256 */ M_E256, /* e^256 */
M_E512, /* e^512 */ M_E512 /* e^512 */
M_E1024, /* e^1024 */
}; };
/**************************************************************************** /****************************************************************************
@@ -82,7 +80,9 @@ double lib_expi(size_t n)
size_t i; size_t i;
double val; double val;
if (n > 1024) /* The largest calculable value for n is floor(ln(DBL_MAX)) */
if (n > 709)
{ {
return INFINITY; return INFINITY;
} }
libc/math/lib_libexpif.c
+5 -3
View File
@@ -58,7 +58,7 @@
*(3.402823e+38). *(3.402823e+38).
*/ */
static const float g_expif_square_tbl[7] = static const float g_expif_square_tbl[] =
{ {
(float)M_E, /* e^1 */ (float)M_E, /* e^1 */
(float)M_E2, /* e^2 */ (float)M_E2, /* e^2 */
@@ -66,7 +66,7 @@ static const float g_expif_square_tbl[7] =
(float)M_E8, /* e^8 */ (float)M_E8, /* e^8 */
(float)M_E16, /* e^16 */ (float)M_E16, /* e^16 */
(float)M_E32, /* e^32 */ (float)M_E32, /* e^32 */
(float)M_E64, /* e^64 */ (float)M_E64 /* e^64 */
}; };
/**************************************************************************** /****************************************************************************
@@ -78,7 +78,9 @@ float lib_expif(size_t n)
size_t i; size_t i;
float val; float val;
if (n >= 128) /* The largest calculable value for n is floor(ln(FLT_MAX)) */
if (n > 88)
{ {
return INFINITY_F; return INFINITY_F;
} }