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: j1(x),y1(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 j1l = __j1l
  40 #pragma weak y1l = __y1l
  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 invsqrtpi= 5.641895835477562869480794515607725858441e-0001L,
  53 tpi     = 0.636619772367581343075535053490057448L;
  54 
  55 static GENERIC pone(GENERIC);
  56 static GENERIC qone(GENERIC);
  57 
  58 static const GENERIC r0[7] = {
  59   -6.249999999999999999999999999999999627320e-0002L,
  60    1.940606727194041716205384618494641565464e-0003L,

  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: j1(x),y1(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 __j1l = j1l
  40 #pragma weak __y1l = y1l
  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 invsqrtpi= 5.641895835477562869480794515607725858441e-0001L,
  53 tpi     = 0.636619772367581343075535053490057448L;
  54 
  55 static GENERIC pone(GENERIC);
  56 static GENERIC qone(GENERIC);
  57 
  58 static const GENERIC r0[7] = {
  59   -6.249999999999999999999999999999999627320e-0002L,
  60    1.940606727194041716205384618494641565464e-0003L,