.
[glibc/history.git] / math / w_j1l.c
blob7bd31ad9951d4d77c92eaba2ce4bdc19af370f7e
1 /* w_j1l.c -- long double version of w_j1.c.
2 * Conversion to long double by Ulrich Drepper,
3 * Cygnus Support, drepper@cygnus.com.
4 */
6 /*
7 * ====================================================
8 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
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 * ====================================================
17 #if defined(LIBM_SCCS) && !defined(lint)
18 static char rcsid[] = "$NetBSD: $";
19 #endif
22 * wrapper of j1l,y1l
25 #include <math.h>
26 #include "math_private.h"
28 #ifdef __STDC__
29 long double __j1l(long double x) /* wrapper j1l */
30 #else
31 long double __j1l(x) /* wrapper j1l */
32 long double x;
33 #endif
35 #ifdef _IEEE_LIBM
36 return __ieee754_j1l(x);
37 #else
38 long double z;
39 z = __ieee754_j1l(x);
40 if(_LIB_VERSION == _IEEE_ || __isnanl(x) ) return z;
41 if(fabsl(x)>X_TLOSS) {
42 return __kernel_standard(x,x,236); /* j1(|x|>X_TLOSS) */
43 } else
44 return z;
45 #endif
47 weak_alias (__j1l, j1l)
49 #ifdef __STDC__
50 long double __y1l(long double x) /* wrapper y1l */
51 #else
52 long double __y1l(x) /* wrapper y1l */
53 long double x;
54 #endif
56 #ifdef _IEEE_LIBM
57 return __ieee754_y1l(x);
58 #else
59 long double z;
60 z = __ieee754_y1l(x);
61 if(_LIB_VERSION == _IEEE_ || __isnanl(x) ) return z;
62 if(x <= 0.0){
63 if(x==0.0)
64 /* d= -one/(x-x); */
65 return __kernel_standard(x,x,210);
66 else
67 /* d = zero/(x-x); */
68 return __kernel_standard(x,x,211);
70 if(x>X_TLOSS) {
71 return __kernel_standard(x,x,237); /* y1(x>X_TLOSS) */
72 } else
73 return z;
74 #endif
76 weak_alias (__y1l, y1l)