Blame math/w_lgammal_r_compat.c

Packit 6c4009
/* w_lgammal_r.c -- long double version of w_lgamma_r.c.
Packit 6c4009
 * Conversion to long double by Ulrich Drepper,
Packit 6c4009
 * Cygnus Support, drepper@cygnus.com.
Packit 6c4009
 */
Packit 6c4009
Packit 6c4009
/*
Packit 6c4009
 * ====================================================
Packit 6c4009
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
Packit 6c4009
 *
Packit 6c4009
 * Developed at SunPro, a Sun Microsystems, Inc. business.
Packit 6c4009
 * Permission to use, copy, modify, and distribute this
Packit 6c4009
 * software is freely granted, provided that this notice
Packit 6c4009
 * is preserved.
Packit 6c4009
 * ====================================================
Packit 6c4009
 */
Packit 6c4009
Packit 6c4009
/*
Packit 6c4009
 * wrapper long double lgammal_r(long double x, int *signgamp)
Packit 6c4009
 */
Packit 6c4009
Packit 6c4009
#include <math.h>
Packit 6c4009
#include <math_private.h>
Packit 6c4009
#include <math-svid-compat.h>
Packit 6c4009
#include <libm-alias-ldouble.h>
Packit 6c4009
Packit 6c4009
Packit 6c4009
#if LIBM_SVID_COMPAT
Packit 6c4009
long double
Packit 6c4009
__lgammal_r(long double x, int *signgamp)
Packit 6c4009
{
Packit 6c4009
	long double y = __ieee754_lgammal_r(x,signgamp);
Packit 6c4009
	if(__builtin_expect(!isfinite(y), 0)
Packit 6c4009
	   && isfinite(x) && _LIB_VERSION != _IEEE_)
Packit 6c4009
		return __kernel_standard(x, x,
Packit 6c4009
					 __floorl(x)==x&&x<=0.0
Packit 6c4009
					 ? 215 /* lgamma pole */
Packit 6c4009
					 : 214); /* lgamma overflow */
Packit 6c4009
Packit 6c4009
	return y;
Packit 6c4009
}
Packit 6c4009
libm_alias_ldouble_r (__lgamma, lgamma, _r)
Packit 6c4009
#endif