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 j1 = __j1
  40 #pragma weak y1 = __y1
  41 
  42 #include "libm.h"
  43 #include "libm_synonyms.h"
  44 #include "libm_protos.h"
  45 #include <math.h>
  46 #include <values.h>
  47 
  48 #define GENERIC double
  49 static const GENERIC
  50 zero    = 0.0,
  51 small   = 1.0e-5,
  52 tiny    = 1.0e-20,
  53 one     = 1.0,
  54 invsqrtpi = 5.641895835477562869480794515607725858441e-0001,
  55 tpi     = 0.636619772367581343075535053490057448;
  56 
  57 static GENERIC pone(GENERIC), qone(GENERIC);
  58 static const GENERIC r0[4] = {
  59         -6.250000000000002203053200981413218949548e-0002,
  60         1.600998455640072901321605101981501263762e-0003,
  61         -1.963888815948313758552511884390162864930e-0005,
  62         8.263917341093549759781339713418201620998e-0008,
  63 };

  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 __j1 = j1
  40 #pragma weak __y1 = y1
  41 
  42 #include "libm.h"

  43 #include "libm_protos.h"
  44 #include <math.h>
  45 #include <values.h>
  46 
  47 #define GENERIC double
  48 static const GENERIC
  49 zero    = 0.0,
  50 small   = 1.0e-5,
  51 tiny    = 1.0e-20,
  52 one     = 1.0,
  53 invsqrtpi = 5.641895835477562869480794515607725858441e-0001,
  54 tpi     = 0.636619772367581343075535053490057448;
  55 
  56 static GENERIC pone(GENERIC), qone(GENERIC);
  57 static const GENERIC r0[4] = {
  58         -6.250000000000002203053200981413218949548e-0002,
  59         1.600998455640072901321605101981501263762e-0003,
  60         -1.963888815948313758552511884390162864930e-0005,
  61         8.263917341093549759781339713418201620998e-0008,
  62 };