gsl 1.14.7 → 1.15.3

data/ext/sf_dilog.c

@@ -26,11 +26,12 @@ static VALUE rb_gsl_sf_complex_dilog_e(VALUE obj, VALUE r, VALUE theta)
 {
   gsl_sf_result *re, *im;
   VALUE vre, vim;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(r); Need_Float(theta);
   vre = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, re);
   vim = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, im);
-  status = gsl_sf_complex_dilog_e(NUM2DBL(r), NUM2DBL(theta), re, im);
+  /*status =*/ gsl_sf_complex_dilog_e(NUM2DBL(r), NUM2DBL(theta), re, im);
   return rb_ary_new3(2, vre, vim);
 }
 

data/ext/sf_elementary.c

@@ -16,10 +16,11 @@ static VALUE rb_gsl_sf_multiply_e(VALUE obj, VALUE x, VALUE y)
 {
   gsl_sf_result *r;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x); Need_Float(y);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r);
-  status = gsl_sf_multiply_e(NUM2DBL(x), NUM2DBL(y), r);
+  /*status =*/ gsl_sf_multiply_e(NUM2DBL(x), NUM2DBL(y), r);
   return v;
 }
 
@@ -28,11 +29,12 @@ static VALUE rb_gsl_sf_multiply_err_e(VALUE obj, VALUE x, VALUE dx,
 {
   gsl_sf_result *r;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x); Need_Float(y);
   Need_Float(dx); Need_Float(dy);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r);
-  status = gsl_sf_multiply_err_e(NUM2DBL(x), NUM2DBL(dx), 
+  /*status =*/ gsl_sf_multiply_err_e(NUM2DBL(x), NUM2DBL(dx),
 				 NUM2DBL(y), NUM2DBL(dy), r);
   return v;
 }

data/ext/sf_elljac.c

@@ -15,9 +15,10 @@
 static VALUE rb_gsl_sf_elljac_e(VALUE obj, VALUE n, VALUE m)
 {
   double sn, cn, dn;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(n); Need_Float(m);
-  status = gsl_sf_elljac_e(NUM2DBL(n), NUM2DBL(m), &sn, &cn, &dn);
+  /*status =*/ gsl_sf_elljac_e(NUM2DBL(n), NUM2DBL(m), &sn, &cn, &dn);
   return rb_ary_new3(3, rb_float_new(sn), 
 		     rb_float_new(cn), rb_float_new(dn));
 }

data/ext/sf_exp.c

@@ -30,10 +30,11 @@ static VALUE rb_gsl_sf_exp_e10_e(VALUE obj, VALUE x)
 {
   gsl_sf_result_e10 *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
-  status = gsl_sf_exp_e10_e(NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_exp_e10_e(NUM2DBL(x), rslt);
   return v;
 }
 
@@ -52,10 +53,11 @@ static VALUE rb_gsl_sf_exp_mult_e10_e(VALUE obj, VALUE x, VALUE y)
 {
   gsl_sf_result_e10 *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x);   Need_Float(y);
   v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
-  status = gsl_sf_exp_mult_e10_e(NUM2DBL(x), NUM2DBL(y), rslt);
+  /*status =*/ gsl_sf_exp_mult_e10_e(NUM2DBL(x), NUM2DBL(y), rslt);
   return v;
 }
 
@@ -108,10 +110,11 @@ static VALUE rb_gsl_sf_exp_err_e10_e(VALUE obj, VALUE x, VALUE dx)
 {
   gsl_sf_result_e10 *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x);   Need_Float(dx);
   v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
-  status = gsl_sf_exp_err_e10_e(NUM2DBL(x), NUM2DBL(dx), rslt);
+  /*status =*/ gsl_sf_exp_err_e10_e(NUM2DBL(x), NUM2DBL(dx), rslt);
   return v;
 }
 
@@ -120,11 +123,12 @@ static VALUE rb_gsl_sf_exp_mult_err_e(VALUE obj, VALUE x, VALUE dx,
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x);   Need_Float(y);
   Need_Float(dx);   Need_Float(dy);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_exp_mult_err_e(NUM2DBL(x), NUM2DBL(dx), NUM2DBL(y), NUM2DBL(dy), rslt);
+  /*status =*/ gsl_sf_exp_mult_err_e(NUM2DBL(x), NUM2DBL(dx), NUM2DBL(y), NUM2DBL(dy), rslt);
   return v;
 }
 
@@ -133,11 +137,12 @@ static VALUE rb_gsl_sf_exp_mult_err_e10_e(VALUE obj, VALUE x, VALUE dx,
 {
   gsl_sf_result_e10 *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x);   Need_Float(y);
   Need_Float(dx);   Need_Float(dy);
   v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
-  status = gsl_sf_exp_mult_err_e10_e(NUM2DBL(x), NUM2DBL(dx), NUM2DBL(y), NUM2DBL(dy), rslt);
+  /*status =*/ gsl_sf_exp_mult_err_e10_e(NUM2DBL(x), NUM2DBL(dx), NUM2DBL(y), NUM2DBL(dy), rslt);
   return v;
 }
 

data/ext/sf_gamma.c

@@ -35,11 +35,12 @@ static VALUE rb_gsl_sf_lngamma_sgn_e(VALUE obj, VALUE x)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   double sgn;
   Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_lngamma_sgn_e(NUM2DBL(x), rslt, &sgn);
+  /*status =*/ gsl_sf_lngamma_sgn_e(NUM2DBL(x), rslt, &sgn);
   return rb_ary_new3(2, v, rb_float_new(sgn));
 }
 
@@ -165,10 +166,11 @@ static VALUE rb_gsl_sf_lnchoose_e(VALUE obj, VALUE n, VALUE m)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(n); CHECK_FIXNUM(m);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_lnchoose_e(FIX2INT(n), FIX2INT(m), rslt);
+  /*status =*/ gsl_sf_lnchoose_e(FIX2INT(n), FIX2INT(m), rslt);
   return v;
 }
 
@@ -197,10 +199,11 @@ static VALUE rb_gsl_sf_lnpoch_sgn_e(VALUE obj, VALUE a, VALUE x)
   gsl_sf_result *rslt = NULL;
   VALUE v;
   double sgn;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(a); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_lnpoch_sgn_e(NUM2DBL(a), NUM2DBL(x), rslt, &sgn);
+  /*status =*/ gsl_sf_lnpoch_sgn_e(NUM2DBL(a), NUM2DBL(x), rslt, &sgn);
   return rb_ary_new3(2, v, rb_float_new(sgn));
 }
 

data/ext/sf_gegenbauer.c

@@ -50,11 +50,12 @@ static VALUE rb_gsl_sf_gegenpoly_n_e(VALUE obj, VALUE n, VALUE lambda, VALUE x)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(n);
   Need_Float(lambda); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_gegenpoly_n_e(FIX2INT(n), NUM2DBL(lambda), NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_gegenpoly_n_e(FIX2INT(n), NUM2DBL(lambda), NUM2DBL(x), rslt);
 
   return v;
 }

data/ext/sf_hyperg.c

@@ -31,10 +31,11 @@ static VALUE rb_gsl_sf_hyperg_1F1_int_e(VALUE obj, VALUE m, VALUE n, VALUE x)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(m); CHECK_FIXNUM(n); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_hyperg_1F1_int_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_hyperg_1F1_int_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
   return v;
 }
 
@@ -57,10 +58,11 @@ static VALUE rb_gsl_sf_hyperg_U_int_e(VALUE obj, VALUE m, VALUE n, VALUE x)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(m); CHECK_FIXNUM(n); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_hyperg_U_int_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_hyperg_U_int_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
   return v;
 }
 
@@ -68,10 +70,11 @@ static VALUE rb_gsl_sf_hyperg_U_int_e10_e(VALUE obj, VALUE m, VALUE n, VALUE x)
 {
   gsl_sf_result_e10 *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(m); CHECK_FIXNUM(n); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
-  status = gsl_sf_hyperg_U_int_e10_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_hyperg_U_int_e10_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
   return v;
 }
 
@@ -121,10 +124,11 @@ static VALUE rb_gsl_sf_hyperg_2F1_conj_e(VALUE obj, VALUE aR, VALUE aI, VALUE c,
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(aR); Need_Float(aI);  Need_Float(c); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_hyperg_2F1_conj_e(NUM2DBL(aR), NUM2DBL(aI), NUM2DBL(c), NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_hyperg_2F1_conj_e(NUM2DBL(aR), NUM2DBL(aI), NUM2DBL(c), NUM2DBL(x), rslt);
   return v;
 }
 
@@ -137,10 +141,11 @@ static VALUE rb_gsl_sf_hyperg_2F1_renorm_e(VALUE obj, VALUE a, VALUE b, VALUE c,
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(a); Need_Float(b);  Need_Float(c); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_hyperg_2F1_renorm_e(NUM2DBL(a), NUM2DBL(b), NUM2DBL(c), NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_hyperg_2F1_renorm_e(NUM2DBL(a), NUM2DBL(b), NUM2DBL(c), NUM2DBL(x), rslt);
   return v;
 }
 
@@ -153,10 +158,11 @@ static VALUE rb_gsl_sf_hyperg_2F1_conj_renorm_e(VALUE obj, VALUE aR, VALUE aI, V
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(aR); Need_Float(aI);  Need_Float(c); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_hyperg_2F1_conj_renorm_e(NUM2DBL(aR), NUM2DBL(aI), NUM2DBL(c), NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_hyperg_2F1_conj_renorm_e(NUM2DBL(aR), NUM2DBL(aI), NUM2DBL(c), NUM2DBL(x), rslt);
   return v;
 }
 

data/ext/sf_laguerre.c

@@ -78,11 +78,12 @@ static VALUE rb_gsl_sf_laguerre_n_e(VALUE obj, VALUE n, VALUE a, VALUE x)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(n);
   Need_Float(a); Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_laguerre_n_e(FIX2INT(n), NUM2DBL(a), NUM2DBL(x), rslt);
+  /*status =*/ gsl_sf_laguerre_n_e(FIX2INT(n), NUM2DBL(a), NUM2DBL(x), rslt);
   return v;
 }
 

data/ext/sf_legendre.c

@@ -255,12 +255,13 @@ static VALUE rb_gsl_sf_legendre_H3d_e(VALUE obj,VALUE l,  VALUE lambda, VALUE et
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(l);
   Need_Float(lambda);
   Need_Float(eta);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_legendre_H3d_e(FIX2INT(l), NUM2DBL(lambda), NUM2DBL(eta), rslt);
+  /*status =*/ gsl_sf_legendre_H3d_e(FIX2INT(l), NUM2DBL(lambda), NUM2DBL(eta), rslt);
   return v;
 }
 

data/ext/sf_log.c

@@ -47,7 +47,8 @@ static VALUE rb_gsl_sf_complex_log_e(int argc, VALUE *argv, VALUE obj)
   gsl_complex *z = NULL;
   VALUE vlnr, vtheta;
   double re, im;
-  int status;
+  // local variable "status" was defined and set, but never used
+  //int status;
   switch (argc) {
   case 1:
     CHECK_COMPLEX(argv[0]);
@@ -65,7 +66,7 @@ static VALUE rb_gsl_sf_complex_log_e(int argc, VALUE *argv, VALUE obj)
   }
   vlnr = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt1);
   vtheta = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt2);
-  status = gsl_sf_complex_log_e(re, im, rslt1, rslt2);
+  /*status =*/ gsl_sf_complex_log_e(re, im, rslt1, rslt2);
   return rb_ary_new3(2, vlnr, vtheta);
 }
 

data/ext/sf_power.c

@@ -27,11 +27,12 @@ static VALUE rb_gsl_sf_pow_int_e(VALUE obj, VALUE x, VALUE n)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x);
   CHECK_FIXNUM(n);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_pow_int_e(NUM2DBL(x), FIX2INT(n), rslt);
+  /*status =*/ gsl_sf_pow_int_e(NUM2DBL(x), FIX2INT(n), rslt);
   return v;
 }
 

data/ext/sf_trigonometric.c

@@ -63,7 +63,8 @@ static VALUE rb_gsl_sf_complex_XXX_e(int argc, VALUE *argv, VALUE obj,
   gsl_complex *z;
   double re, im;
   VALUE v1, v2;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   switch (argc) {
   case 1:
     CHECK_COMPLEX(argv[0]);
@@ -82,7 +83,7 @@ static VALUE rb_gsl_sf_complex_XXX_e(int argc, VALUE *argv, VALUE obj,
   }
   v1 = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r1);
   v2 = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r2);
-  status = (*f)(re, im, r1, r2);
+  /*status =*/ (*f)(re, im, r1, r2);
   return rb_ary_new3(2, v1, v2);
 }
 
@@ -157,11 +158,12 @@ static VALUE rb_gsl_sf_sin_err_e(VALUE obj, VALUE x, VALUE dx)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x);
   Need_Float(dx);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_sin_err_e(NUM2DBL(x), NUM2DBL(dx), rslt);
+  /*status =*/ gsl_sf_sin_err_e(NUM2DBL(x), NUM2DBL(dx), rslt);
   return v;
 }
 
@@ -169,11 +171,12 @@ static VALUE rb_gsl_sf_cos_err_e(VALUE obj, VALUE x, VALUE dx)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Need_Float(x);
   Need_Float(dx);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_cos_err_e(NUM2DBL(x), NUM2DBL(dx), rslt);
+  /*status =*/ gsl_sf_cos_err_e(NUM2DBL(x), NUM2DBL(dx), rslt);
   return v;
 }
 

data/ext/signal.c

@@ -105,7 +105,7 @@ static VALUE rb_gsl_fft_conv_corr(int argc, VALUE *argv, VALUE obj,
         enum FFT_CONV_CORR flag2)
 {
   double *data1, *data2, *data3;
-  size_t stride1, stride2, stride3 = 1, size1, size2;
+  size_t stride1, stride2, size1, size2;
 #ifdef HAVE_NARRAY_H
   int naflag1, naflag2, shape;
 #else
@@ -170,14 +170,12 @@ static VALUE rb_gsl_fft_conv_corr(int argc, VALUE *argv, VALUE obj,
       break;
     }
     data3 = v->data;
-    stride3 = 1;
     break;
   case 1:
 #ifdef HAVE_NARRAY_H
     shape = (int) size1;
     ary = na_make_object(NA_DFLOAT, 1, &shape, cNArray);
     data3 = NA_PTR_TYPE(ary, double*);
-    stride3 = 1;
 #endif
     break;
   default:

data/ext/siman.c

@@ -569,7 +569,7 @@ static VALUE rb_gsl_siman_solver_solve(VALUE obj, VALUE rng,
   siman_metric *metric = NULL;
   siman_print *print = NULL;
   gsl_vector *vtmp = NULL;
-  gsl_siman_params_t *params = NULL, ppp;
+  gsl_siman_params_t *params = NULL;
   int flag = 0;
   /*  Data_Get_Struct(obj, siman_solver, ss);*/
   CHECK_VECTOR(vx0p);
@@ -622,7 +622,6 @@ static VALUE rb_gsl_siman_solver_solve(VALUE obj, VALUE rng,
   ss->proc_metric  = metric->proc;
 
   gsl_vector_memcpy(ss->vx, vtmp);
-  ppp = *params;
 
   if (NIL_P(vprint)) {
     gsl_siman_solve(r, ss, rb_gsl_siman_Efunc_t, 

data/ext/stats.c

@@ -45,6 +45,7 @@ static VALUE rb_gsl_stats_XXX1(int argc, VALUE *argv, VALUE obj,
 {
   size_t stride, size;
   double *data = NULL;
+  if (argc < 1) rb_raise(rb_eArgError, "too few arguments");
   data = get_vector_stats2(argc, argv, obj, &stride, &size);
   return rb_float_new((*f)(data, stride, size, NUM2DBL(argv[argc-1])));
 }

data/ext/vector_complex.c

@@ -1986,7 +1986,10 @@ static VALUE rb_gsl_vector_complex_zip(int argc, VALUE *argv, VALUE obj)
   return ary;
 }
 
-static int gsl_vector_complex_equal(const gsl_vector_complex *v1,
+// Starting with version 1.15, GSL provides a gsl_vector_complex_equal
+// function, but it only determines absolute equality (i.e. is has no epsilon
+// argument).
+static int gsl_vector_complex_equal_eps(const gsl_vector_complex *v1,
   const gsl_vector_complex *v2, double eps)
 {
   gsl_complex z1, z2;
@@ -2019,7 +2022,7 @@ static VALUE rb_gsl_vector_complex_equal(int argc, VALUE *argv, VALUE obj)
   Data_Get_Struct(obj, gsl_vector_complex, v1);
   CHECK_VECTOR_COMPLEX(argv[0]);
   Data_Get_Struct(argv[0], gsl_vector_complex, v2);
-  ret = gsl_vector_complex_equal(v1, v2, eps);
+  ret = gsl_vector_complex_equal_eps(v1, v2, eps);
   if (ret == 1) return Qtrue;
   else return Qfalse;
 }