Update.
[glibc.git] / math / w_j1f.c
blob2a7c8db819cf11aad8b1765296e63de7b22bc09d
1 /* w_j1f.c -- float version of w_j1.c.
2 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3 */
5 /*
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
9 * Developed at SunPro, a Sun Microsystems, Inc. business.
10 * Permission to use, copy, modify, and distribute this
11 * software is freely granted, provided that this notice
12 * is preserved.
13 * ====================================================
16 #if defined(LIBM_SCCS) && !defined(lint)
17 static char rcsid[] = "$NetBSD: w_j1f.c,v 1.3 1995/05/10 20:49:17 jtc Exp $";
18 #endif
21 * wrapper of j1f,y1f
24 #include "math.h"
25 #include "math_private.h"
27 #ifdef __STDC__
28 float j1f(float x) /* wrapper j1f */
29 #else
30 float j1f(x) /* wrapper j1f */
31 float x;
32 #endif
34 #ifdef _IEEE_LIBM
35 return __ieee754_j1f(x);
36 #else
37 float z;
38 z = __ieee754_j1f(x);
39 if(_LIB_VERSION == _IEEE_ || __isnanf(x) ) return z;
40 if(fabsf(x)>(float)X_TLOSS) {
41 /* j1(|x|>X_TLOSS) */
42 return (float)__kernel_standard((double)x,(double)x,136);
43 } else
44 return z;
45 #endif
48 #ifdef __STDC__
49 float y1f(float x) /* wrapper y1f */
50 #else
51 float y1f(x) /* wrapper y1f */
52 float x;
53 #endif
55 #ifdef _IEEE_LIBM
56 return __ieee754_y1f(x);
57 #else
58 float z;
59 z = __ieee754_y1f(x);
60 if(_LIB_VERSION == _IEEE_ || __isnanf(x) ) return z;
61 if(x <= (float)0.0){
62 if(x==(float)0.0)
63 /* d= -one/(x-x); */
64 return (float)__kernel_standard((double)x,(double)x,110);
65 else
66 /* d = zero/(x-x); */
67 return (float)__kernel_standard((double)x,(double)x,111);
69 if(x>(float)X_TLOSS) {
70 /* y1(x>X_TLOSS) */
71 return (float)__kernel_standard((double)x,(double)x,137);
72 } else
73 return z;
74 #endif