]> git.ipfire.org Git - thirdparty/glibc.git/blame - math/w_lgammal_r_compat.c
Fix the inaccuracy of j0f/j1f/y0f/y1f [BZ #14469, #14470, #14471, #14472]
[thirdparty/glibc.git] / math / w_lgammal_r_compat.c
CommitLineData
ee188d55
RM
1/* w_lgammal_r.c -- long double version of w_lgamma_r.c.
2 * Conversion to long double by Ulrich Drepper,
3 * Cygnus Support, drepper@cygnus.com.
4 */
5
6/*
7 * ====================================================
8 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
9 *
10 * Developed at SunPro, a Sun Microsystems, Inc. business.
11 * Permission to use, copy, modify, and distribute this
12 * software is freely granted, provided that this notice
13 * is preserved.
14 * ====================================================
15 */
16
ee188d55
RM
17/*
18 * wrapper long double lgammal_r(long double x, int *signgamp)
19 */
20
9d13fb24 21#include <math.h>
9277c064 22#include <math_private.h>
813378e9 23#include <math-svid-compat.h>
a8dce619 24#include <libm-alias-ldouble.h>
ee188d55
RM
25
26
4f3647e4 27#if LIBM_SVID_COMPAT
0ac5ae23
UD
28long double
29__lgammal_r(long double x, int *signgamp)
ee188d55 30{
0ac5ae23 31 long double y = __ieee754_lgammal_r(x,signgamp);
d81f90cc
WD
32 if(__builtin_expect(!isfinite(y), 0)
33 && isfinite(x) && _LIB_VERSION != _IEEE_)
0ac5ae23 34 return __kernel_standard(x, x,
e44acb20 35 floorl(x)==x&&x<=0.0
0ac5ae23
UD
36 ? 215 /* lgamma pole */
37 : 214); /* lgamma overflow */
38
39 return y;
ee188d55 40}
a8dce619 41libm_alias_ldouble_r (__lgamma, lgamma, _r)
4f3647e4 42#endif