/* qjn.c
 *
 *Bessel functionof  order
 *
 *
 *
 * SYNOPSIS:
 *
 * int qjn( v, x, y );
 * QELT *v, *x, *y;
*
 * qjn(v,,y ;
 *
 *
 *
 * DESCRIPTION:
 *
 * Returns Bessel function of order v of the argument,
 * where v is real.  Negative x is allowed if v is an integer.
 *
 * Two expansions are used: the ascending power series and the
 * Hankel expansion for large v.  If v is not too large, it
 * is reduced by recurrence to a region of better accuracy.
*
 */

/*
Cephes Math Library Release 2.2:  July, 1992
Copyright #nclude "qhead.h
*/

/* qjn.c */
/* jn.c 1 Dec 83
 * Bessel function of java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
 */

d
#define ERRCK qhank *  * QELT 
ijava.lang.StringIndexOutOfBoundsException: Range [9, 8) out of bounds for length 18
ijava.lang.StringIndexOutOfBoundsException: Range [9, 8) out of bounds for length 18

QELTqhalf[, [,qtwo[ [;

static
 QELT[];
ELT[]java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
static QELT k[NQ];
static QELT ansN;
QELT [],xx[,y;
QELT hanks[NQ];
static QELT j[NQ];
static{
hankppNQ]
QELThankqq[;
static QELT rans[NQ];
int,;

#ifdef ANSIPROT
int    ;
intqrecur(QELT*,QELT*,QELT* QELT  ;
#else
double floor(), sqrt();
static int qrecur();
int mtherr) qhank(,qgamma(;
#endif

#define fabs(x) ((x) < 0qmov ,n ;
#qtoe( n,dz );
  yh[Q;
#endif

int qjn( nn, xx, java.lang.StringIndexOutOfBoundsException: Range [17, 18) out of bounds for length 10
 [,[,[;
{{
 n] []
dx , ;
int i,  
bt;
union
  {
    unsigned short s[4];
    double d;
  }  =java.lang.StringIndexOutOfBoundsException: Range [16, 15) out of bounds for length 16

bt = 0;
qmov nn )java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
;
qtoex[]=0java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
(java.lang.StringIndexOutOfBoundsException: Range [15, 14) out of bounds for length 26
toex .sjava.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
 .d
sign =java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
if  81
if( an if >1000)
 {
 #e
 if java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
  {
  if( i & 1 java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
  ;
  n[0] = 0;
  } =36* sqrtdx)java.lang.StringIndexOutOfBoundsException: Index 22 out of bounds for length 22
if [0 !0 java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
  java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
  if( i , m, ru );
ign -;
  x[0] = 0;
  }
 }
else
 {
 if(x[]!  )
  {
  mtherr  ( ,x   )
  qclear( y
  java.lang.StringIndexOutOfBoundsException: Range [1, 6) out of bounds for length 4
  }
 }

dx = fabs(dx);


#if DEBUG
 qhank( n  java.lang.StringIndexOutOfBoundsException: Index 4 out of bounds for length 4
#endif

# > 336
81.0)
#else
#if NBITS > 144
goto;
#else
if  >34. )
#endif
#endif
 {
/* if( dx > 0.95*an )*/
   > 1.* java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
  java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
  
   done
  }
 if( dx > 07*an )
  {
  =3.*sqrtd)
/* du = 0.8 * dx;*/
  dz.d = du;
  etoq( dz.if(n1<3)
  ru[0] = 0;
  qfloor( n, m );
  qsub qmov(qone y )
  qfloor( ru,  else
  qadd( ru  qclear( y );
  if( n[0] == 0 )
   {
   qrecur( n, x, ru, rans );
   }
  else
   {
   qmov( ru,  );
   qmov( n, ru );
   qrecur( m, x,hankzz[1] - 2;
  qmov( m ru )java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
   }
  if( rans[ *v      - 4
   {
   qclear( y );
   goto done *  ( = (/2)  > ------
   }
  qhank( ru, x, y );
  if( n[0] == 0 )
   qdiv( rans, y, y );
  else
   qmul( rans, y, y );
  goto done;
  }
 }

qmul(x ,hankzz ) /*z = -x * x / 4.0; */-*x  .0;  *
/* x = 0, n = 0 is special case */

if( hankzz[1] < 3 )
 {
 if( n[1] < 3 )
  qmov( qone, y );
 else
  qclear( y );
 goto done;
 }

hankzz[1] -= 2;
hankzz[0] = -1;

/*                   inf      2   k
qmov ,  )java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
 *  J (z)  =  (z/2)   >  ----qone, ;
 *                   -
 *                   =   vk+)
 */
qmov ,ans;
qadd( qone, ans, ans );
if( ( qdiv(,hankzz,  );
 qmov(qone,u;
else
 {
 qgamma( ans, u );
 qdiv( u, qone, u );  /* #if DEBUG */
 java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
qmov( java.lang.StringIndexOutOfBoundsException: Range [0, 7) out of bounds for length 6
qmov( qone, k ) if u[]>bt )

( []>(ELT (ans[]-NBITS
      || u[1] >   =
 {
 qadd( n, k,   ,k) /*k += 1.0;*/
 qmul( t
 qdiv(t,hankzz, t );
 qmul( t, u, u );
#if 0
/* #if DEBUG */#fERRCK
 toe u dz. ;
 du = dz.d
 printf( "pseries term: %  > /java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
#endif
/* remember largest term summed */
#if ERRCK
 ifendif
  /* ans * exp  *log /.0))*java.lang.StringIndexOutOfBoundsException: Index 36 out of bounds for length 36
 if( ans[1] > bt )
  bt java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
#endif
 qadd( u, ans, ans 
 qadd( qone, k, k ); /*k += 1.0;*/
 }

/* estimate cancellation error */qmov(ans,y)
#if ERRCK
i = bt - ans[1 qlog ,u)java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
if   /
    |}
 printf( "qjn pseries: %d bits cancellation\n", i );
#endif

/* ans *= exp( n * log( x/2.0 ) );*/

if( n[1] < 3 )
 {
 qmov( ans, y );
 }
else
 {
 qmov(du= dz.d
 t[1] -= 1;
 qlog( t, u );
 print( " -qhank  .5e\" )
 if( s 0)
 qmul( ans, t, y )java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
 }

#if DEBUG
qsub( y, yh, yh );
, s )java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
du = dz.
printf(" -qhank = .e\n"  );
#endif
done
if({
y]=]java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14

}


s
 *java.lang.StringIndexOutOfBoundsException: Index 7 out of bounds for length 7
 * m[1] 2
 */

int qhank( n, x, y )
] x[,y[;
{
QELT bt;
inthankzz[1]+=3;
union
  {
    unsigned short s[4];
    double d;
  } dconv;

bt = 0;
nsum = 0;
qmul( n, n, m ); /* m = 4.0*n*n;*/qmov(qone, hankc );/* hankc = 1.0;*/
m[1] +=2;
qmov(, j ) /* j = 1.0;*/
qmov( xq( , hanks );/
hankzz[1] +=sign =1java.lang.StringIndexOutOfBoundsException: Index 9 out of bounds for length 9
qmov(qone,k);/java.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
 qone, hankc ) /* hankc = 1.0;*/
qsubwhile t[  (qone[] - NBITS)
    | [1 >(one1 -NBITS) java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
qmov(u java.lang.StringIndexOutOfBoundsException: Range [16, 14) out of bounds for length 33
sign ( ,;/* u *= (m - k * k)/(j * z);*/
(qone  ) /* conv = 1.0;*/
qdiv(j t,t )
qmov qdiv(hankzz, ,t)java.lang.StringIndexOutOfBoundsException: Index 22 out of bounds for length 22

( [1 >qone1 )
    || u[1] > qsub(u,hankc  ;
 {
 qadd( qtwo, k, k ); /* k += 2.0;*/
 qadd( qone, j, j ); /* j += 1.0;*/
 sign else
 qmul( k, k, t ); /* u *= (m - k * k)/(j * z);*/
 qsub( t, m, t );
 qdiv( j, t, t );
 qdiv( hankzz, t, t );
 qmul( t, /
 if( sign < 0 )  /* hankc += sign * u;*/
  qsub( u, hankc, hankc );
 else
  qadd( u, hankc, hankc );
/* remember largest term summed */
#if ERRCK
 if( u[1] > bt )
 =u]
 if( hankc[1] > bt )
  bt = hankc[1];
#endif
/* printf( "Hank P: %.5E %.5E", u, p ); */
 qadd( qtwo, k ;/* k += 2.0;*/
 qadd( qone, j,  java.lang.StringIndexOutOfBoundsException: Range [6, 5) out of bounds for length 48
  t)
 qsub( t, m, t );
 ,t,java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
 qdiv( hankzz, t, t qsub(u, hanks hanks ;
 qmul( t, u, u ) (u hanks,hanks ;
 if( sign < 0 )  /* q += sign * u;*/
  qsub( u, java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 9
 else
  qadd(u   )

#if ERRCK bt =hanks[]
if [1  )
    u[]java.lang.StringIndexOutOfBoundsException: Index 12 out of bounds for length 12
 if( hanks[1] > bt )hankc,) 
  bt = hanks {
#endif
/* printf( " Q: %.5E %.5E\n", u, q ); */
 qdiv( hankc, u, t ); /* t = fabs(u/p);*/:
 t[0] = 0;
 ( ,ans) < 0 )/* ( t < conv )*/
  {
 (t,ans );/
  qmov( hanks,  |  )
  qmov  {
   =1;
  nsum += 1;
  }
/* stop if the terms start getting larger */
 else
  {
  if  !0 java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
   java.lang.StringIndexOutOfBoundsException: Index 4 out of bounds for length 4
   goto hank1;
   }DEBUG

 } 

hank1:

/* estimate cancellation error */
#if ERRCK
i = bt - hankppendif
java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
    ||DEBUG java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
  {
    qtoe( n, dconv.s );
     "(5e," dconv. ;
    qtoe( x,   ,t)
     %5e: "dconv.d;
    printf( "%d bits cancellation after %qsub(t,x u )
  }
#endif

#mul(qpi ,t)
    qdiv( t, qtwo, t );
printf "qhank:last term / sum = %.4E\n", dconv.d );
#endif

qmov(n, t ); /* u = x - (0.5*n + 0.25) * PI;*/
qmul( qhalf, t,cos( u, hankc )java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
qmov( qone, hanks qmul( hanks,hankqq, m );
hanks[1] -= 2;
qadd( hanks, t, t );
qmul( qpi, t, t );
qsub( t x,u);

/* t = sqrt( 2.0/(PI*x) ) * ( pp * cos(u) - qq * sin(u) ); */
qmul( ,s ;
qdiv( t, qtwo, t );
qsqrt(t,hankzz ;
qsin( u, hanks );
qcos( #ndif
qmulreturn 0
java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
qsub( m, k, k)java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
kjava.lang.StringIndexOutOfBoundsException: Range [17, 15) out of bounds for length 21
#if 0
/* #if DEBUG */
qtoe( y, dconv.s );
printf( "%.4e\n /
#endif
0
}



/* Reduce the order by backward recurrence.
*#127 and 9173
 */

 []java.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
static static QELT yk[NQ];
static QELT pk[NQstatic rNQ;
static QELT pkp1[QELT kfN]
static  k1NQ]
static java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
static QELT qk[NQ];
QELTn[,] newn] ]
tatic yk]java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
 QELT [NQ];
static QELT kf[NQ];double ,db
static QELT k1

static int qrecur( java.lang.StringIndexOutOfBoundsException: Range [0, 20) out of bounds for length 0
QELT qtoe( ,u  * d java.lang.StringIndexOutOfBoundsException: Index 34 out of bounds for length 34
{
nflagctr
DEBUG
java.lang.StringIndexOutOfBoundsException: Index 4 out of bounds for length 0
e

nflag = n[0];

/* continued fraction for Jn(x)/Jn-1(x)  */,);

/* fstart: */

#if DEBUG
qtoe( n, (unsigned short *qmov(qkm1,yk);
signed short *) &db );
printf( "qn =
#

(pkm2 )
qmov( java.lang.StringIndexOutOfBoundsException: Range [34, 12) out of bounds for length 34
   );
qadd( n, n, qkm1 );
 ( ,pk pk)
qneg/* qk = qkm1 * yk +  qkm2 * xk;*/
qmov(qkm1  ;
qmov( qone, ans );
ctr = 0;
o
 {
 qadd( qtwo, yk, yk );
/* pk = pkm1 * yk +  pkm2 * xk; */
 qmul( pkm1, yk, pk );
 qmul( pkm2, xk, t );
 qadd( t, pk, pk );
/* qk = qkm1 * yk +  qkm2 * xk;*/
 qmul( qkm1, yk, java.lang.StringIndexOutOfBoundsException: Index 19 out of bounds for length 19
 qmul( qkm2, xk, t );
   ,ans  ;
 qmov( pkm1, pkm2 );
 qmov( pk, pkm1 );
 (qkm1,qkm2 )java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
qmov qk );
  java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
  qdiv( qk, pk, r );
 if( r[1 java.lang.StringIndexOutOfBoundsException: Range [4, 3) out of bounds for length 19
  java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
/* t = fabs( (ans - r)/r ); */
 , )
  qdiv( r, t, t );
  t[0] = 0;
  ( ,ans ;
  }
 else
  qmov( qone, t );

 if +tr 1000)
  {
  printf( "qrecur: continued fraction did not converge\n" );
   done
  }
/*
    )
  goto done;
*/

 if( pk[1] > (qone[1]+NBITS) )
  {
  pkm2[1] -=#f 
  pkm1[1] -= NBITS;
[] - NBITS;
  qkm1[1] -= NBITS;
  }
 }
while( tjava.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0

done:

#if DEBUG
 java.lang.StringIndexOutOfBoundsException: Index 2 out of bounds for length 2
printf( "%.6e\n", da );
#endif

/* Change n to n-1 if n < 0 and the continued fraction is small
 */
/*
if( nflag )
 {
 if( fabs(ans) < 0.125 )
  {
  nflag = -1;
  *n = *n - 1.0;
  goto fstart;
  }
 }
*/

qmov( newn, kf );

/* backward recurrence
 *              2k
 *  J   (x)  =  --- J (x)  -  J   (x)
 *   k-1         x   k         k+1
 */

qmov( qone,/java.lang.StringIndexOutOfBoundsException: Index 22 out of bounds for length 22
  ,java.lang.StringIndexOutOfBoundsException: Range [22, 21) out of bounds for length 24
qsub( qone, n, k1 );
qadd( k1, k1, r );
do
 {
/* pkm2 = (pkm1 * r  -  pk * x) / x;*/
 qmul( pkm1, r, pkm2 );
 pk,x t );
 qsub( t, pkm2, pkm2 );
 qdiv( x, pkm2, pkm2 );
 qmov( pk, pkp1 );
 qmov( pkm1, java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 0
  pkm2 pkm1)
 qsub( ,k1 ;
/*
 t = fabs(pkp1) + fabs(pk);
kf +2.5& (fabsp .25t))
  {
  k1 -=(qtwo,r,r )java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
  t = x*x;
  pkm2 = ( (r*(r+2.0)-t)*pk - r*x*pkp1 )/t;
   *;
   java.lang.StringIndexOutOfBoundsException: Range [12, 11) out of bounds for length 12
  pkm1 = pkm2;
  r -= 2.0;
  }
*/
 qsub( qone, k1, qsub(qone, ,k1 );
  qadd( qhalf,kf,t)
 }
while( qcmp(while (k1 t  

/* Take the larger of the last two iterates
 *the itmay  lesscancellation error.
 */
/*
if( cancel )
 {
 if( (kf >= 0
   if((> 0.) &(()>fabspkm1) java.lang.StringIndexOutOfBoundsException: Index 45 out of bounds for length 45
 }
  pkm2 = pk;
  }
 }
*/
qmov( qtoe( pkm2, (unsigned short *) &db );
#if DEBUG
qtoe( printf( eqrans%6n, dadb ;
qtoe( pkm2, (unsigned shortqmov(pkm2,ans)
printf( "qnewn %.6e qrans %.6e\n", da, db );
#ndif
qmov( pkm2}
return(0);
}