10  * See the License for the specific language governing permissions
  11  * and limitations under the License.
  12  *
  13  * When distributing Covered Code, include this CDDL HEADER in each
  14  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
  15  * If applicable, add the following below this CDDL HEADER, with the
  16  * fields enclosed by brackets "[]" replaced with your own identifying
  17  * information: Portions Copyright [yyyy] [name of copyright owner]
  18  *
  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 #if defined(ELFOBJ)
  31 #pragma weak jnl = __jnl
  32 #pragma weak ynl = __ynl
  33 #endif
  34 
  35 /*
  36  * floating point Bessel's function of the 1st and 2nd kind
  37  * of order n: jn(n,x),yn(n,x);
  38  *
  39  * Special cases:
  40  *      y0(0)=y1(0)=yn(n,0) = -inf with division by zero signal;
  41  *      y0(-ve)=y1(-ve)=yn(n,-ve) are NaN with invalid signal.
  42  * Note 2. About jn(n,x), yn(n,x)
  43  *      For n=0, j0(x) is called,
  44  *      for n=1, j1(x) is called,
  45  *      for n<x, forward recursion us used starting
  46  *      from values of j0(x) and j1(x).
  47  *      for n>x, a continued fraction approximation to
  48  *      j(n,x)/j(n-1,x) is evaluated and then backward
  49  *      recursion is used starting from a supposed value
  50  *      for j(n,x). The resulting value of j(0,x) is
  51  *      compared with the actual value to correct the
  52  *      supposed value of j(n,x).
  53  *

  10  * See the License for the specific language governing permissions
  11  * and limitations under the License.
  12  *
  13  * When distributing Covered Code, include this CDDL HEADER in each
  14  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
  15  * If applicable, add the following below this CDDL HEADER, with the
  16  * fields enclosed by brackets "[]" replaced with your own identifying
  17  * information: Portions Copyright [yyyy] [name of copyright owner]
  18  *
  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 #pragma weak __jnl = jnl
  31 #pragma weak __ynl = ynl


  32 
  33 /*
  34  * floating point Bessel's function of the 1st and 2nd kind
  35  * of order n: jn(n,x),yn(n,x);
  36  *
  37  * Special cases:
  38  *      y0(0)=y1(0)=yn(n,0) = -inf with division by zero signal;
  39  *      y0(-ve)=y1(-ve)=yn(n,-ve) are NaN with invalid signal.
  40  * Note 2. About jn(n,x), yn(n,x)
  41  *      For n=0, j0(x) is called,
  42  *      for n=1, j1(x) is called,
  43  *      for n<x, forward recursion us used starting
  44  *      from values of j0(x) and j1(x).
  45  *      for n>x, a continued fraction approximation to
  46  *      j(n,x)/j(n-1,x) is evaluated and then backward
  47  *      recursion is used starting from a supposed value
  48  *      for j(n,x). The resulting value of j(0,x) is
  49  *      compared with the actual value to correct the
  50  *      supposed value of j(n,x).
  51  *