POK
/home/jaouen/pok_official/pok/trunk/libpok/libm/atanhf.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 Fri Jan 30 14:41:34 2009 
00015  */
00016 
00017 /* w_atanhf.c -- float version of w_atanh.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 atanhf(x)
00036  */
00037 
00038 #include <libm.h>
00039 #include "math_private.h"
00040 
00041 
00042 float
00043 atanhf(float x)         /* wrapper atanhf */
00044 {
00045 #ifdef _IEEE_LIBM
00046         return __ieee754_atanhf(x);
00047 #else
00048         float z,y;
00049         z = __ieee754_atanhf(x);
00050         if(_LIB_VERSION == _IEEE_ || isnanf(x)) return z;
00051         y = fabsf(x);
00052         if(y>=(float)1.0) {
00053             if(y>(float)1.0)
00054                 /* atanhf(|x|>1) */
00055                 return (float)__kernel_standard((double)x,(double)x,130);
00056             else
00057                 /* atanhf(|x|==1) */
00058                 return (float)__kernel_standard((double)x,(double)x,131);
00059         } else
00060             return z;
00061 #endif
00062 }
00063 
00064 #endif