src/_fastcorr.c
int compute_bessel_zeros(double nu, int N, double*xi){
  int i;
  for(i = 0; i < N; i++){
    if (nu == 0){
      xi[i] = gsl_sf_bessel_zero_J0(i+1);
    }
    else if (nu == 1){
      xi[i] = gsl_sf_bessel_zero_J1(i+1);
    }
    else{//nu is not 0 or 1
      xi[i] = gsl_sf_bessel_zero_Jnu(nu, i+1);
    }
  }
  return 0;
}
/*Compute pi times the zeros of J_nu where nu is an integer.
 */
Compute pi times the zeros of J_nu where nu is an integer.
Compute pi times the zeros of J_nu where nu is an integer .
int compute_bessel_zeros(double nu, int N, double*xi){
  int i;
  for(i = 0; i < N; i++){
    if (nu == 0){
      xi[i] = gsl_sf_bessel_zero_J0(i+1);
    }
    else if (nu == 1){
      xi[i] = gsl_sf_bessel_zero_J1(i+1);
    }
    else{
      xi[i] = gsl_sf_bessel_zero_Jnu(nu, i+1);
    }
  }
  return 0;
}
int compute_bessel_zeros ( double nu , int N , double * xi ) { int i ; for ( i = 0 ; i < N ; i ++ ) { if ( nu == 0 ) { xi [ i ] = gsl_sf_bessel_zero_J0 ( i + 1 ) ; } else if ( nu == 1 ) { xi [ i ] = gsl_sf_bessel_zero_J1 ( i + 1 ) ; } else { xi [ i ] = gsl_sf_bessel_zero_Jnu ( nu , i + 1 ) ; } } return 0 ; }
Compute pi times the zeros of J_nu where nu is an integer.
Compute pi times the zeros of J_nu where nu is an integer .
//nu is not 0 or 1