POK
/home/jaouen/pok_official/pok/trunk/libpok/libm/log10f.c
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