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 }