summaryrefslogtreecommitdiffstats
path: root/lib
diff options
context:
space:
mode:
authorbde <bde@FreeBSD.org>1997-02-24 17:44:08 +0000
committerbde <bde@FreeBSD.org>1997-02-24 17:44:08 +0000
commit139936b34a895520e9ee8b9a6d3f774754cdcc48 (patch)
tree1b9f54c259d2b4a658d83e51242c6a2516d396f0 /lib
parent9e266561b379c627a821de0a2e93f47aa96cff01 (diff)
downloadFreeBSD-src-139936b34a895520e9ee8b9a6d3f774754cdcc48.zip
FreeBSD-src-139936b34a895520e9ee8b9a6d3f774754cdcc48.tar.gz
Use the C library version of log10() instead of the inaccurate formula
log10(x) = log10e * log(x). This fixes some small (one or two ULP) inaccuracies. Found by: ucbtest
Diffstat (limited to 'lib')
-rw-r--r--lib/libF77/d_lg10.c6
1 files changed, 2 insertions, 4 deletions
diff --git a/lib/libF77/d_lg10.c b/lib/libF77/d_lg10.c
index f03ff00..3135881 100644
--- a/lib/libF77/d_lg10.c
+++ b/lib/libF77/d_lg10.c
@@ -1,9 +1,7 @@
#include "f2c.h"
-#define log10e 0.43429448190325182765
-
#ifdef KR_headers
-double log();
+double log10();
double d_lg10(x) doublereal *x;
#else
#undef abs
@@ -11,5 +9,5 @@ double d_lg10(x) doublereal *x;
double d_lg10(doublereal *x)
#endif
{
-return( log10e * log(*x) );
+return( log10(*x) );
}
OpenPOWER on IntegriCloud