Update.
[glibc.git] / math / w_j0.c
blob5a018760bb98d919c48c0c3e0c5ee8a7af3f9bfd
1 /* @(#)w_j0.c 5.1 93/09/24 */
2 /*
3 * ====================================================
4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 * Developed at SunPro, a Sun Microsystems, Inc. business.
7 * Permission to use, copy, modify, and distribute this
8 * software is freely granted, provided that this notice
9 * is preserved.
10 * ====================================================
13 #if defined(LIBM_SCCS) && !defined(lint)
14 static char rcsid[] = "$NetBSD: w_j0.c,v 1.6 1995/05/10 20:49:11 jtc Exp $";
15 #endif
18 * wrapper j0(double x), y0(double x)
21 #include "math.h"
22 #include "math_private.h"
24 #ifdef __STDC__
25 double j0(double x) /* wrapper j0 */
26 #else
27 double j0(x) /* wrapper j0 */
28 double x;
29 #endif
31 #ifdef _IEEE_LIBM
32 return __ieee754_j0(x);
33 #else
34 double z = __ieee754_j0(x);
35 if(_LIB_VERSION == _IEEE_ || __isnan(x)) return z;
36 if(fabs(x)>X_TLOSS) {
37 return __kernel_standard(x,x,34); /* j0(|x|>X_TLOSS) */
38 } else
39 return z;
40 #endif
42 #ifdef NO_LONG_DOUBLE
43 strong_alias (j0, j0l)
44 #endif
47 #ifdef __STDC__
48 double y0(double x) /* wrapper y0 */
49 #else
50 double y0(x) /* wrapper y0 */
51 double x;
52 #endif
54 #ifdef _IEEE_LIBM
55 return __ieee754_y0(x);
56 #else
57 double z;
58 z = __ieee754_y0(x);
59 if(_LIB_VERSION == _IEEE_ || __isnan(x) ) return z;
60 if(x <= 0.0){
61 if(x==0.0)
62 /* d= -one/(x-x); */
63 return __kernel_standard(x,x,8);
64 else
65 /* d = zero/(x-x); */
66 return __kernel_standard(x,x,9);
68 if(x>X_TLOSS) {
69 return __kernel_standard(x,x,35); /* y0(x>X_TLOSS) */
70 } else
71 return z;
72 #endif
74 #ifdef NO_LONG_DOUBLE
75 strong_alias (y0, y0l)
76 #endif