From 5d8d8dd3b68ba4176368c890f38c03b1e6e2be45 Mon Sep 17 00:00:00 2001
From: Bruce Evans <bde@FreeBSD.org>
Date: Mon, 24 Feb 1997 17:44:08 +0000
Subject: [PATCH] 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
---
 lib/libF77/d_lg10.c | 6 ++----
 1 file changed, 2 insertions(+), 4 deletions(-)

diff --git a/lib/libF77/d_lg10.c b/lib/libF77/d_lg10.c
index f03ff0043f8..3135881c56a 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) );
 }