rb-gsl 1.15.3.1 → 1.15.3.2

data/ext/nmf.c

@@ -32,8 +32,8 @@ void pp(const gsl_matrix *m)
 {
   int r, c;
 
-  for(r=0; r<m->size1; r++) {
-    for(c=0; c<m->size2; c++) {
+  for(r=0; r < (int) m->size1; r++) {
+    for(c=0; c < (int) m->size2; c++) {
       printf(" %.2f", gsl_matrix_get(m, r, c));
     }
     printf("\n");
@@ -47,9 +47,9 @@ double difcost(const gsl_matrix *a, const gsl_matrix *b)
   int i, j;
   double dif=0, d;
 
-  for (i=0; i < a->size1; i++)
+  for (i=0; i < (int) a->size1; i++)
   {
-    for (j=0; j < a->size2; j++)
+    for (j=0; j < (int) a->size2; j++)
     {
       d = gsl_matrix_get(a, i, j) - gsl_matrix_get(b, i, j);
       dif += d*d;
@@ -66,9 +66,9 @@ static void initmatrix(gsl_matrix *m, double min, double max)
 
   srand(time(NULL));
 
-  for(i=0; i < m->size1; i++)
+  for(i=0; i < (int) m->size1; i++)
   {
-    for(j=0; j < m->size2; j++)
+    for(j=0; j < (int) m->size2; j++)
     {
       val = min + (int) (max * (rand() / (RAND_MAX + min)));
       gsl_matrix_set(m, i, j, val);

data/ext/ntuple.c

@@ -191,7 +191,7 @@ static VALUE rb_gsl_ntuple_select_fn_set_f(int argc, VALUE *argv, VALUE obj)
     CHECK_PROC(argv[0]);
     rb_ary_store(ary, 0, argv[0]);
     ary2 = rb_ary_new2(argc-1);
-    for (i = 1; i < argc; i++) rb_ary_store(ary2, i-1, argv[i]);
+    for (i = 1; (int) i < argc; i++) rb_ary_store(ary2, i-1, argv[i]);
     rb_ary_store(ary, 1, ary2);
     break;
   }
@@ -236,7 +236,7 @@ static VALUE rb_gsl_ntuple_select_fn_set_params(int argc, VALUE *argv, VALUE obj
     rb_ary_store(ary, 1, argv[0]);
   } else {
     ary2 = rb_ary_new2(argc);
-    for (i = 0; i < argc; i++) rb_ary_store(ary2, i, argv[i]);
+    for (i = 0; (int) i < argc; i++) rb_ary_store(ary2, i, argv[i]);
     rb_ary_store(ary, 1, ary2);
   }
   return obj;
@@ -310,7 +310,7 @@ static VALUE rb_gsl_ntuple_value_fn_set_f(int argc, VALUE *argv, VALUE obj)
     CHECK_PROC(argv[0]);
     rb_ary_store(ary, 0, argv[0]);
     ary2 = rb_ary_new2(argc-1);
-    for (i = 1; i < argc; i++) rb_ary_store(ary2, i-1, argv[i]);
+    for (i = 1; (int) i < argc; i++) rb_ary_store(ary2, i-1, argv[i]);
     rb_ary_store(ary, 1, ary2);
     break;
   }
@@ -353,7 +353,7 @@ static VALUE rb_gsl_ntuple_value_fn_set_params(int argc, VALUE *argv, VALUE obj)
     rb_ary_store(ary, 1, argv[0]);
   } else {
     ary2 = rb_ary_new2(argc);
-    for (i = 0; i < argc; i++) rb_ary_store(ary2, i, argv[i]);
+    for (i = 0; (int) i < argc; i++) rb_ary_store(ary2, i, argv[i]);
     rb_ary_store(ary, 1, ary2);
   }
   return obj;

data/ext/odeiv.c

@@ -223,7 +223,7 @@ static void set_sys(int argc, VALUE *argv, gsl_odeiv_system *sys)
     break;
   default:
     vparams = rb_ary_new2(argc-itmp);
-    for (i = itmp, j = 0; i < argc; i++, j++) rb_ary_store(vparams, j, argv[i]);
+    for (i = itmp, j = 0; (int) i < argc; i++, j++) rb_ary_store(vparams, j, argv[i]);
   }
   dimension = FIX2INT(dim);
   sys->dimension = dimension;
@@ -251,7 +251,7 @@ static VALUE rb_gsl_odeiv_system_set_params(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     vparams = rb_ary_new2(argc);
-    for (i = 0; i < argc; i++) rb_ary_store(vparams, i, argv[i]);
+    for (i = 0; (int) i < argc; i++) rb_ary_store(vparams, i, argv[i]);
   }
   //  rb_ary_delete_at(ary, 3);
   rb_ary_store(ary, 3, vparams); 

data/ext/permutation.c

@@ -101,7 +101,7 @@ static VALUE rb_gsl_permutation_get(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     bnew = gsl_permutation_alloc(argc);
-    for (j = 0; j < argc; j++) {
+    for (j = 0; j < (size_t) argc; j++) {
       i = FIX2INT(argv[j]);
       if (i < 0) k = b->size + i; else k = i;
       bnew->data[j] = b->data[k];

data/ext/poly2.c

@@ -37,7 +37,7 @@ static gsl_poly_int* mygsl_poly_hermite(int n1)
     p2 = gsl_vector_int_calloc(n1 + 1);
     memcpy(p1->data, coef2, 3*sizeof(int));
     memcpy(p2->data, coef1, 2*sizeof(int));
-    for (n = 2; n < n1; n++) {
+    for (n = 2; (int) n < n1; n++) {
       gsl_vector_int_memcpy(p0, p1);
       mygsl_vector_int_shift_scale2(p0, n);
       gsl_vector_int_scale(p2, 2*n);
@@ -77,7 +77,7 @@ static gsl_poly_int* mygsl_poly_cheb(int n1)
     p2 = gsl_vector_int_calloc(n1 + 1);
     memcpy(p1->data, coef2, 3*sizeof(int));
     memcpy(p2->data, coef1, 2*sizeof(int));
-    for (n = 2; n < n1; n++) {
+    for (n = 2; (int) n < n1; n++) {
       gsl_vector_int_memcpy(p0, p1);
       mygsl_vector_int_shift_scale2(p0, n);
 
@@ -116,7 +116,7 @@ static gsl_poly_int* mygsl_poly_chebII(int n1)
     p2 = gsl_vector_int_calloc(n1 + 1);
     memcpy(p1->data, coef2, 3*sizeof(int));
     memcpy(p2->data, coef1, 2*sizeof(int));
-    for (n = 2; n < n1; n++) {
+    for (n = 2; (int) n < n1; n++) {
       gsl_vector_int_memcpy(p0, p1);
       mygsl_vector_int_shift_scale2(p0, n);
       gsl_vector_int_sub(p0, p2);
@@ -148,7 +148,7 @@ static gsl_poly_int* mygsl_poly_laguerre(int n)
     break;
   default:
     k = gsl_sf_fact(n);
-    for (m = 0; m <= n; m++) {
+    for (m = 0; (int) m <= n; m++) {
       val = k*k/gsl_sf_fact(n-m)/gsl_pow_2(gsl_sf_fact(m));
       if (m%2 == 1) val *= -1;
       gsl_vector_int_set(p0, m, val);
@@ -164,7 +164,7 @@ static gsl_poly_int* mygsl_poly_bessel(int n)
   gsl_vector_int *p0;
   if (n < 0) rb_raise(rb_eArgError, "order must be >= 0");
   p0 = gsl_vector_int_calloc(n + 1);
-  for (k = 0; k <= n; k++) {
+  for (k = 0; (int) k <= n; k++) {
     gsl_vector_int_set(p0, k, gsl_sf_fact(n+k)/gsl_sf_fact(n-k)/gsl_sf_fact(k)/((int) pow(2, k)));
   }
   return p0;
@@ -191,7 +191,7 @@ static gsl_poly_int* mygsl_poly_bell(int n1)
   default:
     p1 = gsl_vector_int_calloc(n1 + 1);
     memcpy(p1->data, coef2, 3*sizeof(int));
-    for (n = 2; n < n1; n++) {
+    for (n = 2; (int) n < n1; n++) {
       gsl_vector_int_memcpy(p0, p1);
       mygsl_vector_int_shift(p0, n);
       for (j = 0; j < n; j++) {

data/ext/poly_source.c

@@ -186,7 +186,7 @@ static VALUE rb_gsl_poly_eval_singleton(VALUE klass, VALUE a, VALUE x)
 #ifdef GSL_1_11_LATER
 static VALUE rb_gsl_complex_poly_complex_eval(VALUE a, VALUE b)
 {
-  gsl_vector_complex *coef, *zb, *vnew;
+  gsl_vector_complex *coef = NULL, *zb, *vnew;
   gsl_complex *zc;
   gsl_complex z, *zx, *res;
   VALUE ret;
@@ -220,7 +220,7 @@ static VALUE rb_gsl_complex_poly_complex_eval(VALUE a, VALUE b)
     break;
   case T_ARRAY:
     ret = rb_ary_new2(RARRAY_LEN(b));
-    for (i = 0; i < RARRAY_LEN(b); i++) {
+    for (i = 0; (int) i < RARRAY_LEN(b); i++) {
       Data_Get_Struct(rb_ary_entry(b, i), gsl_complex, zx);
       res = (gsl_complex*) malloc(sizeof(gsl_complex));
       *res = gsl_complex_poly_complex_eval(zc, N, *zx);
@@ -282,7 +282,7 @@ static VALUE FUNCTION(rb_gsl_poly,eval)(VALUE obj, VALUE xx)
     break;
   case T_ARRAY:
     ary = rb_ary_new2(RARRAY_LEN(xx));
-    for (i = 0; i < RARRAY_LEN(xx); i++) {
+    for (i = 0; (int) i < RARRAY_LEN(xx); i++) {
       x = rb_ary_entry(xx, i);
       Need_Float(x);
       rb_ary_store(ary, i, rb_float_new(FUNCTION(gsl_poly,eval)(p->data, p->size, NUM2DBL(x))));
@@ -385,7 +385,7 @@ static VALUE FUNCTION(rb_gsl_poly,eval2)(int argc, VALUE *argv, VALUE obj)
     break;
   case T_ARRAY:
     ary = rb_ary_new2(RARRAY_LEN(xx));
-    for (i = 0; i < RARRAY_LEN(xx); i++) {
+    for (i = 0; (int) i < RARRAY_LEN(xx); i++) {
       x = rb_ary_entry(xx, i);
       Need_Float(x);
       rb_ary_store(ary, i, rb_float_new(FUNCTION(gsl_poly,eval)(p->data, size, NUM2DBL(x))));
@@ -1265,12 +1265,15 @@ GSL_TYPE(gsl_vector)* FUNCTION(gsl_poly,deconv_vector)(const GSL_TYPE(gsl_vector
   FUNCTION(gsl_vector,set)(vnew, n-1, FUNCTION(gsl_vector,get)(c2, c2->size-1)/aa);
   for (i = n - 2, k = 1; k < n; i--, k++) {
     x = FUNCTION(gsl_vector,get)(c2, c2->size-1-k);
-    for (j = n-1; j >= 0; j--) {
+    for (j = n-1;; j--) {
       z = FUNCTION(gsl_vector,get)(vnew, j);
       jj = c2->size-1-k-j;
-      if (jj > k || jj < 0) continue;
+      //if (jj > k || jj < 0) continue;
+      if (jj <= k) {
       y = FUNCTION(gsl_vector,get)(a2, jj);
       x -= y*z;
+      }
+      if (j == 0) break;
     }
     FUNCTION(gsl_vector,set)(vnew, i, x/aa);
   }

data/ext/sf.c

@@ -775,22 +775,22 @@ VALUE rb_gsl_sf_eval_double_m(double (*func)(double, gsl_mode_t), VALUE argv, VA
   VALUE ary, xx;
   size_t i, k, n;
   double val;
-  gsl_mode_t mode;
-  char c;
+  /*gsl_mode_t mode;
+  char c;*/
 #ifdef HAVE_NARRAY_H
   double *ptr1, *ptr2;
   struct NARRAY *na;
 #endif
   switch (TYPE(m)) {
   case T_STRING:
-    c = tolower(NUM2CHR(m));
+    /*c = tolower(NUM2CHR(m));
     if (c == 'd') mode = GSL_PREC_DOUBLE;
     else if (c == 's') mode = GSL_PREC_SINGLE;
     else if (c == 'a') mode = GSL_PREC_APPROX;
-    else mode = GSL_PREC_DOUBLE;
+    else mode = GSL_PREC_DOUBLE;*/
     break;
   case T_FIXNUM:
-    mode = FIX2INT(m);
+    /*mode = FIX2INT(m);*/
     break;
   default:
     rb_raise(rb_eArgError, "wrong type argument %s (String or Fixnum expected)",
@@ -859,19 +859,19 @@ VALUE rb_gsl_sf_eval_double2_m(double (*func)(double, double, gsl_mode_t),
   VALUE ary, xx;
   size_t i, k, n;
   double val, xx2;
-  gsl_mode_t mode;
-  char c;
+  /*gsl_mode_t mode;
+  char c;*/
 #ifdef HAVE_NARRAY_H
   double *ptr1, *ptr2;
   struct NARRAY *na;
 #endif
   Need_Float(x2);
   xx2 = NUM2DBL(x2);
-  c = tolower(NUM2CHR(m));
+  /*c = tolower(NUM2CHR(m));
   if (c == 'd') mode = GSL_PREC_DOUBLE;
   else if (c == 's') mode = GSL_PREC_SINGLE;
   else if (c == 'a') mode = GSL_PREC_APPROX;
-  else mode = GSL_PREC_DOUBLE;
+  else mode = GSL_PREC_DOUBLE;*/
   if (CLASS_OF(argv) == rb_cRange) argv = rb_gsl_range2ary(argv);
   switch (TYPE(argv)) {
   case T_FLOAT:
@@ -935,8 +935,8 @@ VALUE rb_gsl_sf_eval_double3_m(double (*func)(double, double, double, gsl_mode_t
   VALUE ary, xx;
   size_t i, k, n;
   double val, xx2, xx3;
-  gsl_mode_t mode;
-  char c;
+  /*gsl_mode_t mode;
+  char c;*/
 #ifdef HAVE_NARRAY_H
   double *ptr1, *ptr2;
   struct NARRAY *na;
@@ -944,11 +944,11 @@ VALUE rb_gsl_sf_eval_double3_m(double (*func)(double, double, double, gsl_mode_t
   Need_Float(x2); Need_Float(x3);
   xx2 = NUM2DBL(x2);
   xx3 = NUM2DBL(x3);
-  c = tolower(NUM2CHR(m));
+  /*c = tolower(NUM2CHR(m));
   if (c == 'd') mode = GSL_PREC_DOUBLE;
   else if (c == 's') mode = GSL_PREC_SINGLE;
   else if (c == 'a') mode = GSL_PREC_APPROX;
-  else mode = GSL_PREC_DOUBLE;
+  else mode = GSL_PREC_DOUBLE;*/
   if (CLASS_OF(argv) == rb_cRange) argv = rb_gsl_range2ary(argv);
   switch (TYPE(argv)) {
   case T_FLOAT:
@@ -1013,19 +1013,19 @@ VALUE rb_gsl_sf_eval_double4_m(double (*func)(double, double, double, double,
   VALUE ary, xx;
   size_t i, k, n;
   double val, xx2, xx3, xx4;
-  gsl_mode_t mode;
-  char c;
+  /*gsl_mode_t mode;
+  char c;*/
 #ifdef HAVE_NARRAY_H
   double *ptr1, *ptr2;
   struct NARRAY *na;
 #endif
   Need_Float(x2); Need_Float(x3); Need_Float(x4);
   xx2 = NUM2DBL(x2);  xx3 = NUM2DBL(x3); xx4 = NUM2DBL(x4);
-  c = tolower(NUM2CHR(m));
+  /*c = tolower(NUM2CHR(m));
   if (c == 'd') mode = GSL_PREC_DOUBLE;
   else if (c == 's') mode = GSL_PREC_SINGLE;
   else if (c == 'a') mode = GSL_PREC_APPROX;
-  else mode = GSL_PREC_DOUBLE;
+  else mode = GSL_PREC_DOUBLE;*/
   if (CLASS_OF(argv) == rb_cRange) argv = rb_gsl_range2ary(argv);
   switch (TYPE(argv)) {
   case T_FLOAT:
@@ -1086,10 +1086,9 @@ VALUE rb_gsl_sf_eval_e(int (*func)(double, gsl_sf_result*), VALUE x)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2DBL(x), rslt);
+  (*func)(NUM2DBL(x), rslt);
   return v;
 }
 
@@ -1097,10 +1096,9 @@ VALUE rb_gsl_sf_eval_e_int(int (*func)(int, gsl_sf_result*), VALUE x)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   CHECK_FIXNUM(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2INT(x), rslt);
+  (*func)(NUM2INT(x), rslt);
   return v;
 }
 
@@ -1108,9 +1106,8 @@ VALUE rb_gsl_sf_eval_e_uint(int (*func)(unsigned int, gsl_sf_result*), VALUE x)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2UINT(x), rslt);
+  (*func)(NUM2UINT(x), rslt);
   return v;
 }
 
@@ -1119,10 +1116,9 @@ VALUE rb_gsl_sf_eval_e_int_uint(int (*func)(int, unsigned int, gsl_sf_result*),
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   CHECK_FIXNUM(n);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(FIX2INT(n), NUM2UINT(x), rslt);
+  (*func)(FIX2INT(n), NUM2UINT(x), rslt);
   return v;
 }
 
@@ -1131,10 +1127,9 @@ VALUE rb_gsl_sf_eval_e_double_uint(int (*func)(double, unsigned int, gsl_sf_resu
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   Need_Float(y);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2DBL(y), NUM2UINT(x), rslt);
+  (*func)(NUM2DBL(y), NUM2UINT(x), rslt);
   return v;
 }
 
@@ -1143,11 +1138,10 @@ VALUE rb_gsl_sf_eval_e_int_double(int (*func)(int, double, gsl_sf_result*),
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   CHECK_FIXNUM(n);
   Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(FIX2INT(n), NUM2DBL(x), rslt);
+  (*func)(FIX2INT(n), NUM2DBL(x), rslt);
   return v;
 }
 
@@ -1156,11 +1150,10 @@ VALUE rb_gsl_sf_eval_e_int_double2(int (*func)(int, double, double, gsl_sf_resul
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   CHECK_FIXNUM(n);
   Need_Float(x1); Need_Float(x2);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(FIX2INT(n), NUM2DBL(x1), NUM2DBL(x2), rslt);
+  (*func)(FIX2INT(n), NUM2DBL(x1), NUM2DBL(x2), rslt);
   return v;
 }
 
@@ -1170,11 +1163,10 @@ VALUE rb_gsl_sf_eval_e_int_int_double(int (*func)(int, int, double, gsl_sf_resul
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   CHECK_FIXNUM(n1); CHECK_FIXNUM(n2);
   Need_Float(x);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(FIX2INT(n1), FIX2INT(n2), NUM2DBL(x), rslt);
+  (*func)(FIX2INT(n1), FIX2INT(n2), NUM2DBL(x), rslt);
   return v;
 }
 
@@ -1183,10 +1175,9 @@ VALUE rb_gsl_sf_eval_e_double2(int (*func)(double, double, gsl_sf_result*),
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   Need_Float(x1); Need_Float(x2); 
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2DBL(x1), NUM2DBL(x2), rslt);
+  (*func)(NUM2DBL(x1), NUM2DBL(x2), rslt);
   return v;
 }
 
@@ -1196,10 +1187,9 @@ VALUE rb_gsl_sf_eval_e_double3(int (*func)(double, double, double, gsl_sf_result
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   Need_Float(x1); Need_Float(x2); Need_Float(x3);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2DBL(x1), NUM2DBL(x2),NUM2DBL(x3), rslt);
+  (*func)(NUM2DBL(x1), NUM2DBL(x2),NUM2DBL(x3), rslt);
   return v;
 }
 
@@ -1210,7 +1200,6 @@ VALUE rb_gsl_sf_eval_e_m(int (*func)(double, gsl_mode_t, gsl_sf_result*),
   char c;
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   Need_Float(x);
   switch (TYPE(m)) {
   case T_STRING:
@@ -1229,7 +1218,7 @@ VALUE rb_gsl_sf_eval_e_m(int (*func)(double, gsl_mode_t, gsl_sf_result*),
     break;
   }
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2DBL(x), mode, rslt);
+  (*func)(NUM2DBL(x), mode, rslt);
   return v;
 }
 
@@ -1241,7 +1230,6 @@ VALUE rb_gsl_sf_eval_e_double2_m(int (*func)(double, double, gsl_mode_t, gsl_sf_
   char c;
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   Need_Float(x1);  Need_Float(x2);
   switch (TYPE(m)) {
   case T_STRING:
@@ -1260,7 +1248,7 @@ VALUE rb_gsl_sf_eval_e_double2_m(int (*func)(double, double, gsl_mode_t, gsl_sf_
     break;
   }
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2DBL(x1), NUM2DBL(x2), mode, rslt);
+  (*func)(NUM2DBL(x1), NUM2DBL(x2), mode, rslt);
   return v;
 }
 
@@ -1271,7 +1259,6 @@ VALUE rb_gsl_sf_eval_e_double3_m(int (*func)(double, double, double, gsl_mode_t,
   char c;
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   Need_Float(x1); Need_Float(x2); Need_Float(x3);
   switch (TYPE(m)) {
   case T_STRING:
@@ -1290,7 +1277,7 @@ VALUE rb_gsl_sf_eval_e_double3_m(int (*func)(double, double, double, gsl_mode_t,
     break;
   }
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2DBL(x1), NUM2DBL(x2),NUM2DBL(x3), mode, rslt);
+  (*func)(NUM2DBL(x1), NUM2DBL(x2),NUM2DBL(x3), mode, rslt);
   return v;
 }
 
@@ -1302,7 +1289,6 @@ VALUE rb_gsl_sf_eval_e_double4_m(int (*func)(double, double, double, double, gsl
   char c;
   gsl_sf_result *rslt = NULL;
   VALUE v;
-  int status;
   Need_Float(x1); Need_Float(x2); Need_Float(x3); Need_Float(x4);
   switch (TYPE(m)) {
   case T_STRING:
@@ -1321,7 +1307,7 @@ VALUE rb_gsl_sf_eval_e_double4_m(int (*func)(double, double, double, double, gsl
     break;
   }    
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  status = (*func)(NUM2DBL(x1), NUM2DBL(x2),NUM2DBL(x3), NUM2DBL(x4), mode, rslt);
+  (*func)(NUM2DBL(x1), NUM2DBL(x2),NUM2DBL(x3), NUM2DBL(x4), mode, rslt);
   return v;
 }