src/_fastcorr.c
int compute_weights(double nu, int N, double*xi, double*w){
  int i;
  double two_over_pi2 = 2./(PI*PI);
  for(i = 0; i < N; i++){
    if (nu == 0){
      w[i] = two_over_pi2 / (xi[i] * gsl_sf_bessel_J1(PI*xi[i]));
    }else{ //nu > 0
      w[i] = two_over_pi2 / (xi[i] * gsl_sf_bessel_Jnu(nu+1, PI*xi[i]));
    }
  }
  return 0;
}
/*Compute weights of the quadrature routine.
 */
Compute weights of the quadrature routine.
Compute weights of the quadrature routine .
int compute_weights(double nu, int N, double*xi, double*w){
  int i;
  double two_over_pi2 = 2./(PI*PI);
  for(i = 0; i < N; i++){
    if (nu == 0){
      w[i] = two_over_pi2 / (xi[i] * gsl_sf_bessel_J1(PI*xi[i]));
    }else{ 
      w[i] = two_over_pi2 / (xi[i] * gsl_sf_bessel_Jnu(nu+1, PI*xi[i]));
    }
  }
  return 0;
}
int compute_weights ( double nu , int N , double * xi , double * w ) { int i ; double two_over_pi2 = 2. / ( PI * PI ) ; for ( i = 0 ; i < N ; i ++ ) { if ( nu == 0 ) { w [ i ] = two_over_pi2 / ( xi [ i ] * gsl_sf_bessel_J1 ( PI * xi [ i ] ) ) ; } else { w [ i ] = two_over_pi2 / ( xi [ i ] * gsl_sf_bessel_Jnu ( nu + 1 , PI * xi [ i ] ) ) ; } } return 0 ; }
Compute weights of the quadrature routine.
Compute weights of the quadrature routine .
//nu > 0