1/* w_j0f.c -- float version of w_j0.c.
2 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3 */
4
5/*
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 *
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 * ====================================================
14 */
15
16#include <sys/cdefs.h>
17#if defined(LIBM_SCCS) && !defined(lint)
18__RCSID("$NetBSD: w_j0f.c,v 1.5 1999/07/02 15:37:44 simonb Exp $");
19#endif
20
21/*
22 * wrapper j0f(float x), y0f(float x)
23 */
24
25#include "math.h"
26#include "math_private.h"
27
28float
29j0f(float x)		/* wrapper j0f */
30{
31#ifdef _IEEE_LIBM
32	return __ieee754_j0f(x);
33#else
34	float z = __ieee754_j0f(x);
35	if(_LIB_VERSION == _IEEE_ || isnanf(x)) return z;
36	if(fabsf(x)>(float)X_TLOSS) {
37		/* j0f(|x|>X_TLOSS) */
38	        return (float)__kernel_standard((double)x,(double)x,134);
39	} else
40	    return z;
41#endif
42}
43
44float
45y0f(float x)		/* wrapper y0f */
46{
47#ifdef _IEEE_LIBM
48	return __ieee754_y0f(x);
49#else
50	float z;
51	z = __ieee754_y0f(x);
52	if(_LIB_VERSION == _IEEE_ || isnanf(x) ) return z;
53        if(x <= (float)0.0){
54                if(x==(float)0.0)
55                    /* d= -one/(x-x); */
56                    return (float)__kernel_standard((double)x,(double)x,108);
57                else
58                    /* d = zero/(x-x); */
59                    return (float)__kernel_standard((double)x,(double)x,109);
60        }
61	if(x>(float)X_TLOSS) {
62		/* y0(x>X_TLOSS) */
63	        return (float)__kernel_standard((double)x,(double)x,135);
64	} else
65	    return z;
66#endif
67}
68