* sysdeps/unix/sysv/linux/powerpc/powerpc64/sysdep.h
[glibc.git] / math / w_j1.c
blobdb8c4ce333143377985b3a53a522d31c1f07a411
1 /* @(#)w_j1.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_j1.c,v 1.6 1995/05/10 20:49:15 jtc Exp $";
15 #endif
18 * wrapper of j1,y1
21 #include <math.h>
22 #include "math_private.h"
24 #ifdef __STDC__
25 double j1(double x) /* wrapper j1 */
26 #else
27 double j1(x) /* wrapper j1 */
28 double x;
29 #endif
31 #ifdef _IEEE_LIBM
32 return __ieee754_j1(x);
33 #else
34 double z;
35 z = __ieee754_j1(x);
36 if(_LIB_VERSION == _IEEE_ || __isnan(x) ) return z;
37 if(fabs(x)>X_TLOSS) {
38 return __kernel_standard(x,x,36); /* j1(|x|>X_TLOSS) */
39 } else
40 return z;
41 #endif
43 #ifdef NO_LONG_DOUBLE
44 strong_alias (j1, j1l)
45 #endif
48 #ifdef __STDC__
49 double y1(double x) /* wrapper y1 */
50 #else
51 double y1(x) /* wrapper y1 */
52 double x;
53 #endif
55 #ifdef _IEEE_LIBM
56 return __ieee754_y1(x);
57 #else
58 double z;
59 z = __ieee754_y1(x);
60 if(_LIB_VERSION == _IEEE_ || __isnan(x) ) return z;
61 if(x <= 0.0){
62 if(x==0.0)
63 /* d= -one/(x-x); */
64 return __kernel_standard(x,x,10);
65 else
66 /* d = zero/(x-x); */
67 return __kernel_standard(x,x,11);
69 if(x>X_TLOSS) {
70 return __kernel_standard(x,x,37); /* y1(x>X_TLOSS) */
71 } else
72 return z;
73 #endif
75 #ifdef NO_LONG_DOUBLE
76 strong_alias (y1, y1l)
77 #endif