POK
|
00001 /* 00002 * POK header 00003 * 00004 * The following file is a part of the POK project. Any modification should 00005 * made according to the POK licence. You CANNOT use this file or a part of 00006 * this file is this part of a file for your own project 00007 * 00008 * For more information on the POK licence, please see our LICENCE FILE 00009 * 00010 * Please follow the coding guidelines described in doc/CODING_GUIDELINES 00011 * 00012 * Copyright (c) 2007-2009 POK team 00013 * 00014 * Created by julien on Sat Jan 31 15:19:57 2009 00015 */ 00016 00017 /* w_log10f.c -- float version of w_log10.c. 00018 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. 00019 */ 00020 00021 /* 00022 * ==================================================== 00023 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 00024 * 00025 * Developed at SunPro, a Sun Microsystems, Inc. business. 00026 * Permission to use, copy, modify, and distribute this 00027 * software is freely granted, provided that this notice 00028 * is preserved. 00029 * ==================================================== 00030 */ 00031 00032 #ifdef POK_NEEDS_LIBMATH 00033 00034 /* 00035 * wrapper log10f(X) 00036 */ 00037 00038 #include <libm.h> 00039 #include "math_private.h" 00040 00041 00042 float 00043 log10f(float x) /* wrapper log10f */ 00044 { 00045 #ifdef _IEEE_LIBM 00046 return __ieee754_log10f(x); 00047 #else 00048 float z; 00049 z = __ieee754_log10f(x); 00050 if(_LIB_VERSION == _IEEE_ || isnanf(x)) return z; 00051 if(x<=(float)0.0) { 00052 if(x==(float)0.0) 00053 /* log10(0) */ 00054 return (float)__kernel_standard((double)x,(double)x,118); 00055 else 00056 /* log10(x<0) */ 00057 return (float)__kernel_standard((double)x,(double)x,119); 00058 } else 00059 return z; 00060 #endif 00061 } 00062 00063 #endif