From 139936b34a895520e9ee8b9a6d3f774754cdcc48 Mon Sep 17 00:00:00 2001
From: bde <bde@FreeBSD.org>
Date: Mon, 24 Feb 1997 17:44:08 +0000
Subject: 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(-)

(limited to 'lib/libF77')

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) );
 }
-- 
cgit v1.1