Subversion Repositories Kolibri OS

Rev

Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
4973 right-hear 1
/* Copyright (C) 1994 DJ Delorie, see COPYING.DJ for details */
2
/* @(#)w_j1.c 5.1 93/09/24 */
3
/*
4
 * ====================================================
5
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6
 *
7
 * Developed at SunPro, a Sun Microsystems, Inc. business.
8
 * Permission to use, copy, modify, and distribute this
9
 * software is freely granted, provided that this notice
10
 * is preserved.
11
 * ====================================================
12
 */
13
 
14
#if defined(LIBM_SCCS) && !defined(lint)
15
static char rcsid[] = "$Id: w_j1.c,v 1.4 1994/08/10 20:34:35 jtc Exp $";
16
#endif
17
 
18
/*
19
 * wrapper of j1,y1
20
 */
21
 
22
#include "math.h"
23
#include "math_private.h"
24
 
25
#ifdef __STDC__
26
	double j1(double x)		/* wrapper j1 */
27
#else
28
	double j1(x)			/* wrapper j1 */
29
	double x;
30
#endif
31
{
32
#ifdef _IEEE_LIBM
33
	return __ieee754_j1(x);
34
#else
35
	double z;
36
	z = __ieee754_j1(x);
37
	if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z;
38
	if(fabs(x)>X_TLOSS) {
39
	        return __kernel_standard(x,x,36); /* j1(|x|>X_TLOSS) */
40
	} else
41
	    return z;
42
#endif
43
}
44
 
45
#ifdef __STDC__
46
	double y1(double x)		/* wrapper y1 */
47
#else
48
	double y1(x)			/* wrapper y1 */
49
	double x;
50
#endif
51
{
52
#ifdef _IEEE_LIBM
53
	return __ieee754_y1(x);
54
#else
55
	double z;
56
	z = __ieee754_y1(x);
57
	if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z;
58
        if(x <= 0.0){
59
                if(x==0.0)
60
                    /* d= -one/(x-x); */
61
                    return __kernel_standard(x,x,10);
62
                else
63
                    /* d = zero/(x-x); */
64
                    return __kernel_standard(x,x,11);
65
        }
66
	if(x>X_TLOSS) {
67
	        return __kernel_standard(x,x,37); /* y1(x>X_TLOSS) */
68
	} else
69
	    return z;
70
#endif
71
}