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

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) );
 }