19 * CDDL HEADER END
20 */
21
22 /*
23 * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
24 */
25 /*
26 * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
27 * Use is subject to license terms.
28 */
29
30 /*
31 * Floating point Bessel's function of the first and second kinds
32 * of order zero: j0(x),y0(x);
33 *
34 * Special cases:
35 * y0(0)=y1(0)=yn(n,0) = -inf with division by zero signal;
36 * y0(-ve)=y1(-ve)=yn(n,-ve) are NaN with invalid signal.
37 */
38
39 #pragma weak j0l = __j0l
40 #pragma weak y0l = __y0l
41
42 #include "libm.h"
43 #include "longdouble.h"
44
45 #define GENERIC long double
46 static const GENERIC
47 zero = 0.0L,
48 small = 1.0e-9L,
49 tiny = 1.0e-38L,
50 one = 1.0L,
51 five = 5.0L,
52 eight = 8.0L,
53 invsqrtpi= 5.641895835477562869480794515607725858441e-0001L,
54 tpi = 0.636619772367581343075535053490057448L;
55
56 static GENERIC pzero(GENERIC);
57 static GENERIC qzero(GENERIC);
58
59 static GENERIC r0[7] = {
60 -2.499999999999999999999999999999998934492e-0001L,
|
19 * CDDL HEADER END
20 */
21
22 /*
23 * Copyright 2011 Nexenta Systems, Inc. All rights reserved.
24 */
25 /*
26 * Copyright 2006 Sun Microsystems, Inc. All rights reserved.
27 * Use is subject to license terms.
28 */
29
30 /*
31 * Floating point Bessel's function of the first and second kinds
32 * of order zero: j0(x),y0(x);
33 *
34 * Special cases:
35 * y0(0)=y1(0)=yn(n,0) = -inf with division by zero signal;
36 * y0(-ve)=y1(-ve)=yn(n,-ve) are NaN with invalid signal.
37 */
38
39 #pragma weak __j0l = j0l
40 #pragma weak __y0l = y0l
41
42 #include "libm.h"
43 #include "longdouble.h"
44
45 #define GENERIC long double
46 static const GENERIC
47 zero = 0.0L,
48 small = 1.0e-9L,
49 tiny = 1.0e-38L,
50 one = 1.0L,
51 five = 5.0L,
52 eight = 8.0L,
53 invsqrtpi= 5.641895835477562869480794515607725858441e-0001L,
54 tpi = 0.636619772367581343075535053490057448L;
55
56 static GENERIC pzero(GENERIC);
57 static GENERIC qzero(GENERIC);
58
59 static GENERIC r0[7] = {
60 -2.499999999999999999999999999999998934492e-0001L,
|