/* vi: set sw=4 ts=4: */
/*
 * Wrapper functions implementing all the float math functions
 * defined by SuSv3 by actually calling the double version of
 * each function and then casting the result back to a float
 * to return to the user.
 *
 * Copyright (C) 2005 by Erik Andersen <andersen@uclibc.org>
 *
 * This program is free software; you can redistribute it and/or modify it
 * under the terms of the GNU Library General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or (at your
 * option) any later version.
 *
 * This program is distributed in the hope that it will be useful, but WITHOUT
 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
 * FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public License
 * for more details.
 *
 * You should have received a copy of the GNU Library General Public License
 * along with this program; if not, write to the Free Software Foundation,
 * Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
 */

#include "math.h"

/* For the time being, do _NOT_ implement these functions
 * that are defined by SuSv3 */
#if 0
float       exp2f(float);
float       fmaf(float, float, float);
float       fmaxf(float, float);
float       fminf(float, float);
float       fdimf(float, float);
long long   llrintf(float);
long long   llroundf(float);
long        lroundf(float);
float       log2f(float);
long        lrintf(float);
float       nexttowardf(float, long double);
float       remquof(float, float, int *);
float       roundf(float);
float       scalblnf(float, long);
float       truncf(float);
float       tgammaf(float);
#endif

/* Implement the following, as defined by SuSv3 */
#if 0
float       acosf(float);
float       acoshf(float);
float       asinf(float);
float       asinhf(float);
float       atan2f(float, float);
float       atanf(float);
float       atanhf(float);
float       cbrtf(float);
float       ceilf(float);
float       copysignf(float, float);
float       cosf(float);
float       coshf(float);
float       erfcf(float);
float       erff(float);
float       expf(float);
float       expm1f(float);
float       fabsf(float);
float       floorf(float);
float       fmodf(float, float);
float       frexpf(float value, int *);
float       hypotf(float, float);
int         ilogbf(float);
float       ldexpf(float, int);
float       lgammaf(float);
float       log10f(float);
float       log1pf(float);
float       logbf(float);
float       logf(float);
float       modff(float, float *);
float       nearbyintf(float);
float       nextafterf(float, float);
float       powf(float, float);
float       remainderf(float, float);
float       rintf(float);
float       scalbnf(float, int);
float       sinf(float);
float       sinhf(float);
float       sqrtf(float);
float       tanf(float);
float       tanhf(float);
#endif

#ifdef L_acosf
float acosf (float x)
{
	return (float) acos( (double)x );
}
#endif


#ifdef L_acoshf
float acoshf (float x)
{
	return (float) acosh( (double)x );
}
#endif


#ifdef L_asinf
float asinf (float x)
{
	return (float) asin( (double)x );
}
#endif


#ifdef L_asinhf
float asinhf (float x)
{
	return (float) asinh( (double)x );
}
#endif


#ifdef L_atan2f
float atan2f (float x, float y)
{
	return (float) atan2( (double)x, (double)y );
}
#endif


#ifdef L_atanf
float atanf (float x)
{
	return (float) atan( (double)x );
}
#endif


#ifdef L_atanhf
float atanhf (float x)
{
	return (float) atanh( (double)x );
}
#endif


#ifdef L_cbrtf
float cbrtf (float x)
{
	return (float) cbrt( (double)x );
}
#endif


#ifdef L_ceilf
float ceilf (float x)
{
	return (float) ceil( (double)x );
}
#endif


#ifdef L_copysignf
float copysignf (float x, float y)
{
	return (float) copysign( (double)x, (double)y );
}
#endif


#ifdef L_cosf
float cosf (float x)
{
	return (float) cos( (double)x );
}
#endif


#ifdef L_coshf
float coshf (float x)
{
	return (float) cosh( (double)x );
}
#endif


#ifdef L_erfcf
float erfcf (float x)
{
	return (float) erfc( (double)x );
}
#endif


#ifdef L_erff
float erff (float x)
{
	return (float) erf( (double)x );
}
#endif


#if 0
#ifdef L_exp2f
float exp2f (float x)
{
	return (float) exp2( (double)x );
}
#endif
#endif


#ifdef L_expf
float expf (float x)
{
	return (float) exp( (double)x );
}
#endif


#ifdef L_expm1f
float expm1f (float x)
{
	return (float) expm1( (double)x );
}
#endif


#ifdef L_fabsf
float fabsf (float x)
{
	return (float) fabs( (double)x );
}
#endif


#if 0
#ifdef L_fdimf
float fdimf (float x, float y)
{
	return (float) fdim( (double)x, (double)y );
}
#endif
#endif


#ifdef L_floorf
float floorf (float x)
{
	return (float) floor( (double)x );
}
#endif


#if 0
#ifdef L_fmaf
float fmaf (float x, float y, float z)
{
	return (float) fma( (double)x, (double)y, (double)z );
}
#endif


#ifdef L_fmaxf
float fmaxf (float x, float y)
{
	return (float) fmax( (double)x, (double)y );
}
#endif


#ifdef L_fminf
float fminf (float x, float y)
{
	return (float) fmin( (double)x, (double)y );
}
#endif
#endif


#ifdef L_fmodf
float fmodf (float x, float y)
{
	return (float) fmod( (double)x, (double)y );
}
#endif


#ifdef L_frexpf
float frexpf (float x, int *exp)
{
	return (float) frexp( (double)x, exp );
}
#endif


#ifdef L_hypotf
float hypotf (float x, float y)
{
	return (float) hypot( (double)x, (double)y );
}
#endif


#ifdef L_ilogbf
int ilogbf (float x)
{
	return (float) ilogb( (double)x );
}
#endif


#ifdef L_ldexpf
float ldexpf (float x, int exp)
{
	return (float) ldexp( (double)x, exp );
}
#endif


#ifdef L_lgammaf
float lgammaf (float x)
{
	return (float) lgamma( (double)x );
}
#endif


#if 0
#ifdef L_llrintf
long long llrintf (float x)
{
	return (float) llrint( (double)x );
}
#endif


#ifdef L_llroundf
long long llroundf (float x)
{
	return (float) llround( (double)x );
}
#endif
#endif

#ifdef L_log10f
float log10f (float x)
{
	return (float) log10( (double)x );
}
#endif


#ifdef L_log1pf
float log1pf (float x)
{
	return (float) log1p( (double)x );
}
#endif


#if 0
#ifdef L_log2f
float log2f (float x)
{
	return (float) log2( (double)x );
}
#endif
#endif


#ifdef L_logbf
float logbf (float x)
{
	return (float) logb( (double)x );
}
#endif


#ifdef L_logf
float logf (float x)
{
	return (float) log( (double)x );
}
#endif


#if 0
#ifdef L_lrintf
long lrintf (float x)
{
	return (float) lrint( (double)x );
}
#endif


#ifdef L_lroundf
long lroundf (float x)
{
	return (float) lround( (double)x );
}
#endif
#endif


#ifdef L_modff
float modff (float x, float *iptr)
{
	double y, result;
	result = modf ( x, &y );
	*iptr = (float)y;
	return (float) result;

}
#endif


#if 0
#ifdef L_nearbyintf
float nearbyintf (float x)
{
	return (float) nearbyint( (double)x );
}
#endif
#endif


#ifdef L_nextafterf
float nextafterf (float x, float y)
{
	return (float) nextafter( (double)x, (double)y );
}
#endif


#if 0
#ifdef L_nexttowardf
float nexttowardf (float x, long double y)
{
	return (float) nexttoward( (double)x, (double)y );
}
#endif
#endif

#ifdef L_powf
float powf (float x, float y)
{
	return (float) pow( (double)x, (double)y );
}
#endif


#ifdef L_remainderf
float remainderf (float x, float y)
{
	return (float) remainder( (double)x, (double)y );
}
#endif


#if 0
#ifdef L_remquof
float remquof (float x, float y, int *quo)
{
	return (float) remquo( (double)x, (double)y, quo );
}
#endif
#endif


#ifdef L_rintf
float rintf (float x)
{
	return (float) rint( (double)x );
}
#endif


#if 0
#ifdef L_roundf
float roundf (float x)
{
	return (float) round( (double)x );
}
#endif


#ifdef L_scalblnf
float scalblnf (float x, long exp)
{
	return (float) scalbln( (double)x, exp );
}
#endif
#endif


#ifdef L_scalbnf
float scalbnf (float x, int exp)
{
	return (float) scalbn( (double)x, exp );
}
#endif


#ifdef L_sinf
float sinf (float x)
{
	return (float) sin( (double)x );
}
#endif


#ifdef L_sinhf
float sinhf (float x)
{
	return (float) sinh( (double)x );
}
#endif


#ifdef L_sqrtf
float sqrtf (float x)
{
	return (float) sqrt( (double)x );
}
#endif


#ifdef L_tanf
float tanf (float x)
{
	return (float) tan( (double)x );
}
#endif


#ifdef L_tanhf
float tanhf (float x)
{
	return (float) tanh( (double)x );
}
#endif


#if 0
#ifdef L_tgammaf
float tgammaf (float x)
{
	return (float) tgamma( (double)x );
}
#endif


#ifdef L_truncf
float truncf (float x)
{
	return (float) trunc( (double)x );
}
#endif
#endif