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 "libm_synonyms.h"



  44 #include <math.h>
  45 #if defined(__SUNPRO_C)
  46 #include <sunmath.h>
  47 #endif
  48 
  49 #define GENERIC long double
  50 static GENERIC 
  51 zero    = 0.0L,
  52 small   = 1.0e-9L,
  53 tiny    = 1.0e-38L,
  54 one     = 1.0L,
  55 five    = 5.0L,
  56 eight   = 8.0L,
  57 invsqrtpi= 5.641895835477562869480794515607725858441e-0001L,
  58 tpi     = 0.636619772367581343075535053490057448L;
  59 
  60 static GENERIC pzero(), qzero();
  61 static GENERIC r0[7] = {
  62   -2.499999999999999999999999999999998934492e-0001L,
  63    1.272657927360049786327618451133763714880e-0002L,

  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 "libm_synonyms.h"
  44 
  45 #include "longdouble.h"
  46 
  47 #include <math.h>
  48 #if defined(__SUNPRO_C)
  49 #include <sunmath.h>
  50 #endif
  51 
  52 #define GENERIC long double
  53 static GENERIC 
  54 zero    = 0.0L,
  55 small   = 1.0e-9L,
  56 tiny    = 1.0e-38L,
  57 one     = 1.0L,
  58 five    = 5.0L,
  59 eight   = 8.0L,
  60 invsqrtpi= 5.641895835477562869480794515607725858441e-0001L,
  61 tpi     = 0.636619772367581343075535053490057448L;
  62 
  63 static GENERIC pzero(), qzero();
  64 static GENERIC r0[7] = {
  65   -2.499999999999999999999999999999998934492e-0001L,
  66    1.272657927360049786327618451133763714880e-0002L,