src/_fastcorr.c
int compute_spherical_bessel_zeros(int nu, int N, double*xi){
  int i;
  for(i = 0; i < N; i++){
    //Note: GSL returns PI the root we want
    if (nu == 0){
      xi[i] = i+1.;
    }else{//nu > 0
      xi[i] = inv_PI * gsl_sf_bessel_zero_Jnu((nu+1)/2, 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_spherical_bessel_zeros(int nu, int N, double*xi){
  int i;
  for(i = 0; i < N; i++){
    if (nu == 0){
      xi[i] = i+1.;
    }else{
      xi[i] = inv_PI * gsl_sf_bessel_zero_Jnu((nu+1)/2, i+1);
    }
  }
  return 0;
}
int compute_spherical_bessel_zeros ( int nu , int N , double * xi ) { int i ; for ( i = 0 ; i < N ; i ++ ) { if ( nu == 0 ) { xi [ i ] = i + 1. ; } else { xi [ i ] = inv_PI * gsl_sf_bessel_zero_Jnu ( ( nu + 1 ) / 2 , 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 .
//Note: GSL returns PI the root we want //nu > 0