gsl 1.14.7 → 1.15.3

data/ext/multimin.c

@@ -403,14 +403,12 @@ void rb_gsl_multimin_function_fdf_df(const gsl_vector *x, void *p,
 void rb_gsl_multimin_function_fdf_fdf(const gsl_vector *x, void *p, 
 				      double *f, gsl_vector *g)
 {
-  VALUE vx, vf, vg, proc_fdf, proc_f, proc_df, vp, ary, result;
+  VALUE vx, vg, proc_f, proc_df, vp, ary, result;
   vx = Data_Wrap_Struct(cgsl_vector, 0, NULL, (gsl_vector *) x);
   vg = Data_Wrap_Struct(cgsl_vector, 0, NULL, g);
-  vf = rb_float_new(*f);
   ary = (VALUE) p;
   proc_f = rb_ary_entry(ary, 0);
   proc_df = rb_ary_entry(ary, 1);
-  proc_fdf = rb_ary_entry(ary, 2);
   vp = rb_ary_entry(ary, RARRAY_LEN(ary)-1);
   if (NIL_P(vp)) {
     result = rb_funcall(proc_f, RBGSL_ID_call, 1, vx);

data/ext/odeiv.c

@@ -81,7 +81,8 @@ static int calc_jac(double t, const double y[], double *dfdy, double dfdt[], voi
 static int calc_func(double t, const double y[], double dydt[], void *data)
 {
   VALUE ary, params, proc;
-  VALUE result;
+  // local variable "result" declared and set, but never used
+  //VALUE result;
   VALUE vy, vdydt;
   gsl_vector_view ytmp, dydttmp;
   size_t dim;
@@ -100,9 +101,9 @@ static int calc_func(double t, const double y[], double dydt[], void *data)
   vy = Data_Wrap_Struct(cgsl_vector_view_ro, 0, NULL, &ytmp);
   vdydt = Data_Wrap_Struct(cgsl_vector_view, 0, NULL, &dydttmp);
 
-  if (NIL_P(params)) result = rb_funcall((VALUE) proc, RBGSL_ID_call, 3, rb_float_new(t), 
+  if (NIL_P(params)) /*result =*/ rb_funcall((VALUE) proc, RBGSL_ID_call, 3, rb_float_new(t), 
 					 vy, vdydt);
-  else result = rb_funcall((VALUE) proc, RBGSL_ID_call, 4, rb_float_new(t), vy, vdydt, params);
+  else /*result =*/ rb_funcall((VALUE) proc, RBGSL_ID_call, 4, rb_float_new(t), vy, vdydt, params);
  
   return GSL_SUCCESS;
 }
@@ -110,7 +111,9 @@ static int calc_func(double t, const double y[], double dydt[], void *data)
 static int calc_jac(double t, const double y[], double *dfdy, double dfdt[], void *data)
 {
   VALUE params, proc, ary;
-  VALUE result, vdfdt;
+  VALUE vdfdt;
+  // local variable "result" declared and set, but never used
+  //VALUE result;
   VALUE vy, vmjac;
   gsl_vector_view ytmp, dfdttmp;
   gsl_matrix_view mv;
@@ -133,9 +136,9 @@ static int calc_jac(double t, const double y[], double *dfdy, double dfdt[], voi
   vy = Data_Wrap_Struct(cgsl_vector_view_ro, 0, NULL, &ytmp);
   vmjac = Data_Wrap_Struct(cgsl_matrix_view, 0, NULL, &mv);
   vdfdt = Data_Wrap_Struct(cgsl_vector_view, 0, NULL, &dfdttmp);
-  if (NIL_P(params)) result = rb_funcall((VALUE) proc, RBGSL_ID_call, 4, rb_float_new(t),
+  if (NIL_P(params)) /*result =*/ rb_funcall((VALUE) proc, RBGSL_ID_call, 4, rb_float_new(t),
 					 vy, vmjac, vdfdt);
-  else result = rb_funcall((VALUE) proc, RBGSL_ID_call, 5, rb_float_new(t), 
+  else /*result =*/ rb_funcall((VALUE) proc, RBGSL_ID_call, 5, rb_float_new(t), 
 			   vy, vmjac, vdfdt, params);
   return GSL_SUCCESS;
 }

data/ext/poly_source.c

@@ -466,9 +466,16 @@ static VALUE FUNCTION(rb_gsl_poly,solve_quadratic)(int argc, VALUE *argv, VALUE
 	     "wrong number of arguments (3 numbers or 1 array or 1 vector)");
     break;
   }
-  r = gsl_vector_alloc(2);
-  gsl_vector_set(r, 0, x0);
-  gsl_vector_set(r, 1, x1);
+  // If n == 0, we want to return an empty gsl_vector, but gsl_vectors can'y be
+  // empty, so just return an empty Array.
+  if(n == 0) {
+    return rb_ary_new();
+  }
+  r = gsl_vector_alloc(n);
+  switch(n) {
+    case 2: gsl_vector_set(r, 1, x1); /* fall through */
+    case 1: gsl_vector_set(r, 0, x0);
+  }
   return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, r);
 }
 
@@ -507,9 +514,16 @@ static VALUE FUNCTION(rb_gsl_poly,complex_solve_quadratic)(int argc, VALUE *argv
 	     "wrong number of arguments (3 numbers or 1 array or 1 vector)");
     break;
   }
-  r = gsl_vector_complex_alloc(2);
-  gsl_vector_complex_set(r, 0, z0);
-  gsl_vector_complex_set(r, 1, z1);
+  // If n == 0, we want to return an empty gsl_vector, but gsl_vectors can'y be
+  // empty, so just return an empty Array.
+  if(n == 0) {
+    return rb_ary_new();
+  }
+  r = gsl_vector_complex_alloc(n);
+  switch(n) {
+    case 2: gsl_vector_complex_set(r, 1, z1); /* fall through */
+    case 1: gsl_vector_complex_set(r, 0, z0);
+  }
   return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, r);
 }
 
@@ -547,10 +561,12 @@ static VALUE FUNCTION(rb_gsl_poly,solve_cubic)(int argc, VALUE *argv, VALUE obj)
 	     "wrong number of arguments (3 numbers or 1 array or 1 vector)");
     break;
   }
-  r = gsl_vector_alloc(3);
-  gsl_vector_set(r, 0, x0);
-  gsl_vector_set(r, 1, x1);
-  gsl_vector_set(r, 2, x2);
+  r = gsl_vector_alloc(n);
+  switch(n) {
+    case 3: gsl_vector_set(r, 2, x2); /* fall through */
+    case 2: gsl_vector_set(r, 1, x1); /* fall through */
+    case 1: gsl_vector_set(r, 0, x0);
+  }
   return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, r);
 }
 
@@ -588,10 +604,12 @@ static VALUE FUNCTION(rb_gsl_poly,complex_solve_cubic)(int argc, VALUE *argv, VA
 	     "wrong number of arguments (3 numbers or 1 array or 1 vector)");
     break;
   }
-  r = gsl_vector_complex_alloc(3);
-  gsl_vector_complex_set(r, 0, z0);
-  gsl_vector_complex_set(r, 1, z1);
-  gsl_vector_complex_set(r, 2, z2);
+  r = gsl_vector_complex_alloc(n);
+  switch(n) {
+    case 3: gsl_vector_complex_set(r, 2, z2); /* fall through */
+    case 2: gsl_vector_complex_set(r, 1, z1); /* fall through */
+    case 1: gsl_vector_complex_set(r, 0, z0);
+  }
   return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, r);
 }
 
@@ -633,11 +651,13 @@ static VALUE FUNCTION(rb_gsl_poly,solve_quartic)(int argc, VALUE *argv, VALUE ob
 	     "wrong number of arguments (3 numbers or 1 array or 1 vector)");
     break;
   }
-  r = gsl_vector_alloc(4);
-  gsl_vector_set(r, 0, x0);
-  gsl_vector_set(r, 1, x1);
-  gsl_vector_set(r, 2, x2);
-  gsl_vector_set(r, 3, x3);
+  r = gsl_vector_alloc(n);
+  switch(n) {
+    case 4: gsl_vector_set(r, 3, x3); /* fall through */
+    case 3: gsl_vector_set(r, 2, x2); /* fall through */
+    case 2: gsl_vector_set(r, 1, x1); /* fall through */
+    case 1: gsl_vector_set(r, 0, x0);
+  }
   return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, r);
 }
 
@@ -678,11 +698,13 @@ static VALUE FUNCTION(rb_gsl_poly,complex_solve_quartic)(int argc, VALUE *argv,
 	     "wrong number of arguments (3 numbers or 1 array or 1 vector)");
     break;
   }
-  r = gsl_vector_complex_alloc(4);
-  gsl_vector_complex_set(r, 0, z0);
-  gsl_vector_complex_set(r, 1, z1);
-  gsl_vector_complex_set(r, 2, z2);
-  gsl_vector_complex_set(r, 3, z3);
+  r = gsl_vector_complex_alloc(n);
+  switch(n) {
+    case 4: gsl_vector_complex_set(r, 0, z0); /* fall through */
+    case 3: gsl_vector_complex_set(r, 1, z1); /* fall through */
+    case 2: gsl_vector_complex_set(r, 2, z2); /* fall through */
+    case 1: gsl_vector_complex_set(r, 3, z3);
+  }
   return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, r);
 }
 
@@ -696,7 +718,9 @@ VALUE FUNCTION(rb_gsl_poly,complex_solve)(int argc, VALUE *argv, VALUE obj)
   gsl_vector *a = NULL, *z = NULL;
   gsl_complex c;
   gsl_vector_complex *r = NULL;
-  int status, i, flag = 0;
+  int i, flag = 0;
+  // local variable "status" declared and set, but never used
+  //int status;
 
   switch (argc) {
   case 1:
@@ -747,7 +771,7 @@ VALUE FUNCTION(rb_gsl_poly,complex_solve)(int argc, VALUE *argv, VALUE obj)
     w = gsl_poly_complex_workspace_alloc(size);
     flag = 1;
   }
-  status = gsl_poly_complex_solve(a->data, size, w, z->data);
+  /*status =*/ gsl_poly_complex_solve(a->data, size, w, z->data);
 
   if (flag == 1) gsl_poly_complex_workspace_free(w);
   gsl_vector_free(a);
@@ -783,15 +807,19 @@ static VALUE FUNCTION(rb_gsl_poly,solve_quadratic2)(VALUE obj)
   d = a1*a1 - 4.0*a2*a0;    /* determinant */
   if (d >= 0.0) {
     n = gsl_poly_solve_quadratic(a2, a1, a0, &x0, &x1);
-    r = gsl_vector_alloc(2);
-    gsl_vector_set(r, 0, x0);
-    gsl_vector_set(r, 1, x1);
+    r = gsl_vector_alloc(n);
+    switch(n) {
+      case 2: gsl_vector_set(r, 1, x1); /* fall through */
+      case 1: gsl_vector_set(r, 0, x0);
+    }
     return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, r);
   } else {
     n = gsl_poly_complex_solve_quadratic(a2, a1, a0, &z0, &z1);
-    r2 = gsl_vector_complex_alloc(2);
-    gsl_vector_complex_set(r2, 0, z0);
-    gsl_vector_complex_set(r2, 1, z1);
+    r2 = gsl_vector_complex_alloc(n);
+    switch(n) {
+      case 2: gsl_vector_complex_set(r2, 1, z1); /* fall through */
+      case 1: gsl_vector_complex_set(r2, 0, z0);
+    }
     return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, r2);
   }
 
@@ -803,7 +831,6 @@ static VALUE FUNCTION(rb_gsl_poly,complex_solve_quadratic2)(VALUE obj)
   GSL_TYPE(gsl_vector) *v = NULL;
   gsl_vector_complex *r = NULL;
   double a2, a1, a0;
-  double d;                /* determinant */
   int n;
   gsl_complex z0, z1;
   Data_Get_Struct(obj, GSL_TYPE(gsl_vector), v);
@@ -813,11 +840,12 @@ static VALUE FUNCTION(rb_gsl_poly,complex_solve_quadratic2)(VALUE obj)
   a2 = FUNCTION(gsl_vector,get)(v, 2);   /* coefficients */
   a1 = FUNCTION(gsl_vector,get)(v, 1);
   a0 = FUNCTION(gsl_vector,get)(v, 0);
-  d = a1*a1 - 4.0*a2*a0;    /* determinant */
   n = gsl_poly_complex_solve_quadratic(a2, a1, a0, &z0, &z1);
-  r = gsl_vector_complex_alloc(2);
-  gsl_vector_complex_set(r, 0, z0);
-  gsl_vector_complex_set(r, 1, z1);
+  r = gsl_vector_complex_alloc(n);
+  switch(n) {
+    case 2: gsl_vector_complex_set(r, 1, z1); /* fall through */
+    case 1: gsl_vector_complex_set(r, 0, z0);
+  }
   return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, r);
 }
 
@@ -839,10 +867,12 @@ static VALUE FUNCTION(rb_gsl_poly,solve_cubic2)(VALUE obj)
   a1 = FUNCTION(gsl_vector,get)(v, 1)/a3;
   a0 = FUNCTION(gsl_vector,get)(v, 0)/a3;
   n = gsl_poly_solve_cubic(a2, a1, a0, &x0, &x1, &x2);
-  r = gsl_vector_alloc(3);
-  gsl_vector_set(r, 0, x0);
-  gsl_vector_set(r, 1, x1);
-  gsl_vector_set(r, 2, x2);
+  r = gsl_vector_alloc(n);
+  switch(n) {
+    case 3: gsl_vector_set(r, 2, x2); /* fall through */
+    case 2: gsl_vector_set(r, 1, x1); /* fall through */
+    case 1: gsl_vector_set(r, 0, x0);
+  }
   return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, r);
 }
 
@@ -862,10 +892,12 @@ static VALUE FUNCTION(rb_gsl_poly,complex_solve_cubic2)(VALUE obj)
   a1 = FUNCTION(gsl_vector,get)(v, 1)/a3;
   a0 = FUNCTION(gsl_vector,get)(v, 0)/a3;
   n = gsl_poly_complex_solve_cubic(a2, a1, a0, &z0, &z1, &z2);
-  r = gsl_vector_complex_alloc(3);
-  gsl_vector_complex_set(r, 0, z0);
-  gsl_vector_complex_set(r, 1, z1);
-  gsl_vector_complex_set(r, 2, z2);
+  r = gsl_vector_complex_alloc(n);
+  switch(n) {
+    case 3: gsl_vector_complex_set(r, 2, z2); /* fall through */
+    case 2: gsl_vector_complex_set(r, 1, z1); /* fall through */
+    case 1: gsl_vector_complex_set(r, 0, z0);
+  }
   return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, r);
 }
 
@@ -934,7 +966,9 @@ VALUE FUNCTION(rb_gsl_poly,complex_solve2)(int argc, VALUE *argv, VALUE obj)
   gsl_poly_complex_workspace *w = NULL;
   gsl_complex c;
   gsl_vector_complex *r = NULL;
-  int n, flag = 0;
+  int flag = 0;
+  // local variable "status" declared and set, but never used
+  //int status;
 
   Data_Get_Struct(obj, GSL_TYPE(gsl_vector), v);
   size = v->size;
@@ -953,7 +987,7 @@ VALUE FUNCTION(rb_gsl_poly,complex_solve2)(int argc, VALUE *argv, VALUE obj)
     flag = 1;
   }
   
-  n = gsl_poly_complex_solve(a->data, size, w, z->data);
+  /*status =*/ gsl_poly_complex_solve(a->data, size, w, z->data);
 
   r = gsl_vector_complex_alloc(size - 1);
   for (i = 0; i < size - 1; i++) {
@@ -1584,16 +1618,12 @@ static VALUE rb_gsl_poly_wfit(int argc, VALUE *argv, VALUE obj)
   gsl_multifit_linear_workspace *space = NULL;
   gsl_matrix *X = NULL, *cov = NULL;
   gsl_vector *x, *y = NULL, *w, *c = NULL;
-  gsl_vector_view xx, yy, ww;
   size_t order, i, j;
   double chisq, val;
   int status, flag = 0;
   VALUE vc, vcov;
   if (argc != 4 && argc != 5) 
     rb_raise(rb_eArgError, "wrong number of arguments (%d for 4 or 5)", argc);
-  x = &xx.vector;
-  w = &ww.vector;
-  y = &yy.vector;
   Data_Get_Vector(argv[0], x);
   Data_Get_Vector(argv[1], w);
   Data_Get_Vector(argv[2], y);
@@ -1615,7 +1645,7 @@ static VALUE rb_gsl_poly_wfit(int argc, VALUE *argv, VALUE obj)
       gsl_matrix_set(X, i, j, val);
     }
   }
-  status = gsl_multifit_linear(X, y, c, cov, &chisq, space);
+  status = gsl_multifit_wlinear(X, w, y, c, cov, &chisq, space);
   if (flag == 1) gsl_multifit_linear_free(space);
   vc = Data_Wrap_Struct(cgsl_poly, 0, gsl_vector_free, c);
   vcov = Data_Wrap_Struct(cgsl_matrix, 0, gsl_matrix_free, cov);

data/ext/sf_bessel.c

@@ -288,14 +288,16 @@ static VALUE rb_gsl_sf_bessel_jl_e(VALUE obj, VALUE n, VALUE x)
 static VALUE rb_gsl_sf_bessel_xl_array(VALUE obj, VALUE n1, VALUE x,
 				       int (*f)(int, double, double[]))
 {
-  int nmax, n, status;
+  int nmax, n;
+  // local variable "status" declared and set, but never used
+  //int status;
   gsl_vector *v = NULL;
   CHECK_FIXNUM(n1);
   Need_Float(x);
   nmax = FIX2INT(n1);
   n = nmax  + 1;
   v = gsl_vector_alloc(n);
-  status = (*f)(nmax, NUM2DBL(x), v->data);
+  /*status =*/ (*f)(nmax, NUM2DBL(x), v->data);
   return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
 }
 
@@ -465,7 +467,9 @@ static VALUE rb_gsl_sf_bessel_Jnu_e(VALUE obj, VALUE nu, VALUE x)
 static VALUE rb_gsl_sf_bessel_sequence_Jnu_e(int argc, VALUE *argv, VALUE obj)
 {
   size_t i, size;
-  int status, flag = 0;
+  // local variable "status" declared and set, but never used
+  //int status;
+  int flag = 0;
   gsl_vector *v = NULL;
   gsl_mode_t mode;
   char c;
@@ -515,7 +519,7 @@ static VALUE rb_gsl_sf_bessel_sequence_Jnu_e(int argc, VALUE *argv, VALUE obj)
     flag = 0;
     break;
   }
-  status = gsl_sf_bessel_sequence_Jnu_e(NUM2DBL(nu), mode, size, v->data);
+  /*status =*/ gsl_sf_bessel_sequence_Jnu_e(NUM2DBL(nu), mode, size, v->data);
   if (flag == 1) return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
   else return ary;
 }

data/ext/sf_coulomb.c

@@ -35,11 +35,12 @@ static VALUE rb_gsl_sf_hydrogenicR_e(VALUE obj, VALUE n, VALUE l,
 {
   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(l);
   Need_Float(Z); Need_Float(r);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_hydrogenicR_e(FIX2INT(n), FIX2INT(l),
+  /*status =*/ gsl_sf_hydrogenicR_e(FIX2INT(n), FIX2INT(l),
 			    NUM2DBL(Z), NUM2DBL(r), rslt);
   return v;
 }
@@ -87,7 +88,7 @@ static VALUE rb_gsl_sf_coulomb_wave_FG_array(VALUE obj, VALUE Lmin, VALUE kmax,
 {
   double F_exponent, G_exponent;
   int status;
-  size_t size, stride;
+  size_t size;
   gsl_vector *vf = NULL, *vg = NULL;
   VALUE fary, gary;
   CHECK_FIXNUM(kmax);
@@ -95,7 +96,6 @@ static VALUE rb_gsl_sf_coulomb_wave_FG_array(VALUE obj, VALUE Lmin, VALUE kmax,
   size = FIX2INT(kmax);
   vf = gsl_vector_alloc(size);
   vg = gsl_vector_alloc(size);
-  stride = vf->stride;
 
   status = gsl_sf_coulomb_wave_FG_array(NUM2DBL(Lmin), size, NUM2DBL(eta), 
 				   NUM2DBL(x), vf->data, vg->data,
@@ -112,7 +112,7 @@ static VALUE rb_gsl_sf_coulomb_wave_FGp_array(VALUE obj, VALUE Lmin, VALUE kmax,
 {
   double F_exponent, G_exponent;
   int status;
-  size_t size, stride;
+  size_t size;
   gsl_vector *vf = NULL, *vg = NULL, *vfp = NULL, *vgp = NULL;
   VALUE fary, gary, fpary, gpary;
   CHECK_FIXNUM(kmax);
@@ -122,7 +122,6 @@ static VALUE rb_gsl_sf_coulomb_wave_FGp_array(VALUE obj, VALUE Lmin, VALUE kmax,
   vfp = gsl_vector_alloc(size);
   vg = gsl_vector_alloc(size);
   vgp = gsl_vector_alloc(size);
-  stride = vf->stride;
 
   status = gsl_sf_coulomb_wave_FGp_array(NUM2DBL(Lmin), size, NUM2DBL(eta), 
 				     NUM2DBL(x), vf->data, vfp->data, 
@@ -165,12 +164,13 @@ static VALUE rb_gsl_sf_coulomb_CL_array(VALUE obj, VALUE Lmin, VALUE kmax,
 {
   gsl_vector *v = NULL;
   size_t size;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(kmax);
   Need_Float(Lmin); Need_Float(eta); 
   size = FIX2INT(kmax);
   v = gsl_vector_alloc(size);
-  status = gsl_sf_coulomb_CL_array(NUM2DBL(Lmin), size, NUM2DBL(eta), v->data);
+  /*status =*/ gsl_sf_coulomb_CL_array(NUM2DBL(Lmin), size, NUM2DBL(eta), v->data);
   return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
 }
 

data/ext/sf_coupling.c

@@ -25,11 +25,12 @@ static VALUE rb_gsl_sf_coupling_3j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALU
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
   CHECK_FIXNUM(two_ma); CHECK_FIXNUM(two_mb); CHECK_FIXNUM(two_mc);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_coupling_3j_e(FIX2INT(two_ja), FIX2INT(two_jb),
+  /*status =*/ gsl_sf_coupling_3j_e(FIX2INT(two_ja), FIX2INT(two_jb),
 			  FIX2INT(two_jc), FIX2INT(two_ma),
 			  FIX2INT(two_mb), FIX2INT(two_mc),
 			  rslt);
@@ -50,11 +51,12 @@ static VALUE rb_gsl_sf_coupling_6j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALU
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
   CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_coupling_6j_e(FIX2INT(two_ja), FIX2INT(two_jb),
+  /*status =*/ gsl_sf_coupling_6j_e(FIX2INT(two_ja), FIX2INT(two_jb),
 			  FIX2INT(two_jc), FIX2INT(two_jd),
 			  FIX2INT(two_je), FIX2INT(two_jf),
 			  rslt);
@@ -83,12 +85,13 @@ static VALUE rb_gsl_sf_coupling_9j_e(VALUE obj, VALUE two_ja, VALUE two_jb,
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
   CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
   CHECK_FIXNUM(two_jg); CHECK_FIXNUM(two_jh); CHECK_FIXNUM(two_ji);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = gsl_sf_coupling_9j_e(FIX2INT(two_ja), FIX2INT(two_jb),
+  /*status =*/ gsl_sf_coupling_9j_e(FIX2INT(two_ja), FIX2INT(two_jb),
 			   FIX2INT(two_jc), FIX2INT(two_jd),
 			   FIX2INT(two_je), FIX2INT(two_jf),
 			   FIX2INT(two_jg), FIX2INT(two_jh),