]> git.ipfire.org Git - thirdparty/glibc.git/blame - math/w_hypotl_compat.c
Fix the inaccuracy of j0f/j1f/y0f/y1f [BZ #14469, #14470, #14471, #14472]
[thirdparty/glibc.git] / math / w_hypotl_compat.c
CommitLineData
ee188d55
RM
1/* w_hypotl.c -- long double version of w_hypot.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 hypotl(x,y)
19 */
20
9d13fb24 21#include <math.h>
9277c064 22#include <math_private.h>
813378e9 23#include <math-svid-compat.h>
92892fdb 24#include <libm-alias-ldouble.h>
ee188d55
RM
25
26
4f3647e4 27#if LIBM_SVID_COMPAT
0ac5ae23
UD
28long double
29__hypotl(long double x, long double y)
ee188d55 30{
ee188d55
RM
31 long double z;
32 z = __ieee754_hypotl(x,y);
d81f90cc
WD
33 if(__builtin_expect(!isfinite(z), 0)
34 && isfinite(x) && isfinite(y) && _LIB_VERSION != _IEEE_)
41bf21a1 35 return __kernel_standard_l(x, y, 204); /* hypot overflow */
0ac5ae23
UD
36
37 return z;
ee188d55 38}
92892fdb 39libm_alias_ldouble (__hypot, hypot)
4f3647e4 40#endif