freebsd-dev/lib/libF77/d_lg10.c
Bruce Evans 5d8d8dd3b6 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
1997-02-24 17:44:08 +00:00

14 lines
180 B
C

#include "f2c.h"
#ifdef KR_headers
double log10();
double d_lg10(x) doublereal *x;
#else
#undef abs
#include "math.h"
double d_lg10(doublereal *x)
#endif
{
return( log10(*x) );
}