00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
#if defined(LIBM_SCCS) && !defined(lint)
00014
static char rcsid[] =
"$NetBSD: w_log.c,v 1.6 1995/05/10 20:49:33 jtc Exp $";
00015
#endif
00016
00017
00018
00019
00020
00021
#include "math.h"
00022
#include "mathP.h"
00023
00024
00025
#ifdef __STDC__
00026
double log(
double x)
00027 #
else
00028 double log(x)
00029 double x;
00030 #endif
00031 {
00032
#ifdef _IEEE_LIBM
00033
return __ieee754_log(x);
00034
#else
00035
double z;
00036 z =
__ieee754_log(x);
00037
if(
_LIB_VERSION == _IEEE_ || isnan(x) || x > 0.0)
return z;
00038
if(x==0.0)
00039
return __kernel_standard(x,x,16);
00040
else
00041
return __kernel_standard(x,x,17);
00042
#endif
00043
}