gsl 1.14.7 → 1.15.3

data/ext/ieee.c

@@ -21,6 +21,7 @@ static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
 {
 #ifdef RUBY_1_9_LATER
   rb_io_t *fptr = NULL;
+  double ftmp;
 #else
   OpenFile *fptr = NULL;
 #endif
@@ -61,7 +62,8 @@ static VALUE rb_gsl_ieee_fprintf_double(int argc, VALUE *argv, VALUE obj)
     rb_raise(rb_eTypeError, "wrong argument type %s (Float expected)",
 	     rb_class2name(CLASS_OF(vtmp)));
 #ifdef RUBY_1_9_LATER
-  gsl_ieee_fprintf_double(fp, &(RFLOAT_VALUE(vtmp)));
+  ftmp = RFLOAT_VALUE(vtmp);
+  gsl_ieee_fprintf_double(fp, &ftmp);
 #else
   gsl_ieee_fprintf_double(fp, &(RFLOAT(vtmp)->value));
 #endif

data/ext/integration.c

@@ -267,7 +267,9 @@ static VALUE rb_gsl_integration_qng(int argc, VALUE *argv, VALUE obj)
   double result, abserr;
   size_t neval;
   gsl_function *F = NULL;
-  int status, itmp;
+  int status;
+  // local variable 'itmp' declared and set, but never used
+  //int itmp;
 
   if (argc < 1) rb_raise(rb_eArgError, 
 			 "wrong number of arguments (%d for >= 1)", argc);
@@ -276,10 +278,10 @@ static VALUE rb_gsl_integration_qng(int argc, VALUE *argv, VALUE obj)
   case T_MODULE:  case T_CLASS:  case T_OBJECT:
     CHECK_FUNCTION(argv[0]);
     Data_Get_Struct(argv[0], gsl_function, F);
-    itmp = get_a_b_epsabs_epsrel(argc, argv, 1, &a, &b, &epsabs, &epsrel);
+    /*itmp =*/ get_a_b_epsabs_epsrel(argc, argv, 1, &a, &b, &epsabs, &epsrel);
     break;
   default:
-    itmp = get_a_b_epsabs_epsrel(argc, argv, 0, &a, &b, &epsabs, &epsrel);
+    /*itmp =*/ get_a_b_epsabs_epsrel(argc, argv, 0, &a, &b, &epsabs, &epsrel);
     Data_Get_Struct(obj, gsl_function, F);
     break;
   }

data/ext/linalg.c

@@ -666,7 +666,6 @@ static VALUE rb_gsl_linalg_QR_LQ_decomposition(int argc, VALUE *argv, VALUE obj,
   gsl_vector *tau = NULL;
   int (*fdecomp)(gsl_matrix *, gsl_vector *);
   int itmp, status;
-  size_t size;
   VALUE vtau, mdecomp, omatrix;
 
   switch (TYPE(obj)) {
@@ -712,7 +711,6 @@ static VALUE rb_gsl_linalg_QR_LQ_decomposition(int argc, VALUE *argv, VALUE obj,
     rb_raise(rb_eRuntimeError, "unknown operation");
     break;
   }
-  size = m->size1;
   switch (argc - itmp) {
   case 0:
     tau = gsl_vector_alloc(GSL_MIN(mtmp->size1, mtmp->size2));
@@ -931,7 +929,7 @@ static VALUE rb_gsl_linalg_QR_LQ_lssolve(int argc, VALUE *argv, VALUE obj, int f
   gsl_matrix *m = NULL;
   gsl_vector *b = NULL, *x = NULL, *tau = NULL, *r = NULL;
   VALUE omatrix;
-  int flagm = 0, flagt = 0, flagb = 0, flagx = 0, flagr = 0, itmp, status;
+  int flagm = 0, flagt = 0, flagb = 0, itmp, status;
   size_t size;
   int (*fdecomp)(gsl_matrix*, gsl_vector*);
   int (*flssolve)(const gsl_matrix*, const gsl_vector*, const gsl_vector*, gsl_vector*,
@@ -991,22 +989,17 @@ static VALUE rb_gsl_linalg_QR_LQ_lssolve(int argc, VALUE *argv, VALUE obj, int f
   case 2:
     CHECK_VECTOR(argv[argc-2]);
     Data_Get_Struct(argv[argc-2], gsl_vector, x);
-    flagx = 0;
     CHECK_VECTOR(argv[argc-1]);
     Data_Get_Struct(argv[argc-1], gsl_vector, r);
-    flagr = 0;
     break;
   case 1:
     CHECK_VECTOR(argv[argc-1]);
     Data_Get_Struct(argv[argc-1], gsl_vector, x);
-    flagx = 0;
     r = gsl_vector_alloc(x->size);
-    flagr = 1;
     break;
   case 0:
     x = gsl_vector_alloc(m->size1);
     r = gsl_vector_alloc(m->size1);
-    flagx = 1; flagr = 1;
     break;
   default:
     rb_raise(rb_eArgError, "wrong number of arguments");
@@ -2133,7 +2126,7 @@ static VALUE rb_gsl_linalg_QRLQPT_RLsvx(int argc, VALUE *argv, VALUE obj, int fl
   gsl_matrix *QR = NULL;
   gsl_vector *b = NULL;
   gsl_permutation *p = NULL;
-  int itmp, flagb = 0;
+  int itmp;
   VALUE vtmp, klass;
   int (*fsvx)(const gsl_matrix*, const gsl_permutation*, gsl_vector*);
   switch (flag) {
@@ -2174,7 +2167,6 @@ static VALUE rb_gsl_linalg_QRLQPT_RLsvx(int argc, VALUE *argv, VALUE obj, int fl
   itmp++;
   if (TYPE(argv[itmp]) == T_ARRAY) {
     b = make_cvector_from_rarray(argv[itmp]);
-    flagb = 1;
   } else {
     CHECK_VECTOR(argv[itmp]);
     Data_Get_Struct(argv[itmp], gsl_vector, b);
@@ -2853,7 +2845,8 @@ static VALUE rb_gsl_linalg_bidiag_decomp(int argc, VALUE *argv, VALUE obj)
   gsl_matrix *A = NULL, *Atmp = NULL;
   gsl_vector *tau_U = NULL, *tau_V = NULL;
   size_t size0;
-  int status;
+  // local variable "status" was defined and set, but never used
+  //int status;
   VALUE vu, vv, vA;
 
   switch (TYPE(obj)) {
@@ -2870,7 +2863,7 @@ static VALUE rb_gsl_linalg_bidiag_decomp(int argc, VALUE *argv, VALUE obj)
   size0 = GSL_MIN(A->size1, A->size2);
   tau_U = gsl_vector_alloc(size0);
   tau_V = gsl_vector_alloc(size0-1);
-  status = gsl_linalg_bidiag_decomp(A, tau_U, tau_V);
+  /*status =*/ gsl_linalg_bidiag_decomp(A, tau_U, tau_V);
   vA = Data_Wrap_Struct(cgsl_matrix, 0, gsl_matrix_free, A);
   vu = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, tau_U);
   vv = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, tau_V);
@@ -3320,9 +3313,10 @@ static VALUE rb_gsl_linalg_balance_columns_bang(int argc, VALUE *argv, VALUE obj
   gsl_matrix *A = NULL;
   gsl_vector *D = NULL;
   VALUE mat, vec;
-  int status;
+  // local variable "status" was defined and set, but never used
+  //int status;
   rb_gsl_linalg_balance_columns_init(argc, argv, obj, &mat, &vec, &A, &D);
-  status = gsl_linalg_balance_columns(A, D);
+  /*status =*/ gsl_linalg_balance_columns(A, D);
   return rb_ary_new3(2, mat, vec);
 }
 
@@ -3331,11 +3325,12 @@ static VALUE rb_gsl_linalg_balance_columns(int argc, VALUE *argv, VALUE obj)
   gsl_matrix *A = NULL, *Anew;
   gsl_vector *D = NULL;
   VALUE mat, vec;
-  int status;
+  // local variable "status" was defined and set, but never used
+  //int status;
   rb_gsl_linalg_balance_columns_init(argc, argv, obj, &mat, &vec, &A, &D);
   Anew = make_matrix_clone(A);
   mat = Data_Wrap_Struct(cgsl_matrix, 0, gsl_matrix_free, Anew);
-  status = gsl_linalg_balance_columns(Anew, D);
+  /*status =*/ gsl_linalg_balance_columns(Anew, D);
   return rb_ary_new3(2, mat, vec);
 }
 

data/ext/linalg_complex.c

@@ -362,7 +362,7 @@ static VALUE rb_gsl_linalg_complex_LU_lndet(int argc, VALUE *argv, VALUE obj)
   gsl_matrix_complex *m = NULL, *mtmp = NULL;
   gsl_permutation *p = NULL;
   double lndet;
-  int flagm = 0, signum, itmp;
+  int flagm = 0, signum;
   switch (TYPE(obj)) {
   case T_MODULE:
   case T_CLASS:
@@ -376,7 +376,6 @@ static VALUE rb_gsl_linalg_complex_LU_lndet(int argc, VALUE *argv, VALUE obj)
     } else {
       mtmp = m;
     }
-    itmp = 1;
     break;
   default:
     Data_Get_Struct(obj, gsl_matrix_complex, m);
@@ -387,7 +386,6 @@ static VALUE rb_gsl_linalg_complex_LU_lndet(int argc, VALUE *argv, VALUE obj)
     } else {
       mtmp = m;
     }
-    itmp = 0;
   }
   if (flagm == 1) {
     p = gsl_permutation_alloc(m->size1);

data/ext/matrix_complex.c

@@ -199,7 +199,7 @@ static VALUE rb_gsl_matrix_complex_eye(int argc, VALUE *argv, VALUE klass)
 {
   size_t n, i;
   gsl_matrix_complex *m = NULL;
-  gsl_complex z, *p = &z;
+  gsl_complex z, *pz = &z;
   switch (argc) {
   case 1:
     CHECK_FIXNUM(argv[0]);
@@ -223,7 +223,8 @@ static VALUE rb_gsl_matrix_complex_eye(int argc, VALUE *argv, VALUE klass)
       break;
     default:
       if (rb_obj_is_kind_of(argv[1], cgsl_complex)) {
-	Data_Get_Struct(argv[1], gsl_complex, p);
+	Data_Get_Struct(argv[1], gsl_complex, pz);
+        z = *pz;
       } else {
 	rb_raise(rb_eTypeError, 
 		 "wrong argument type %s", rb_class2name(CLASS_OF(argv[1])));
@@ -422,6 +423,7 @@ static VALUE rb_gsl_matrix_complex_set_row(int argc, VALUE *argv, VALUE obj)
     default:
       CHECK_COMPLEX(argv[k]);
       Data_Get_Struct(argv[k], gsl_complex, pz);
+      z = *pz;
       break;
     }
     gsl_matrix_complex_set(A, i, k-1, z);
@@ -448,6 +450,7 @@ static VALUE rb_gsl_matrix_complex_set_col(int argc, VALUE *argv, VALUE obj)
     default:
       CHECK_COMPLEX(argv[k]);
       Data_Get_Struct(argv[k], gsl_complex, pz);
+      z = *pz;
       break;
     }
     gsl_matrix_complex_set(A, k-1, j, z);
@@ -1519,8 +1522,10 @@ static VALUE rb_gsl_matrix_complex_indgen_singleton(int argc, VALUE *argv, VALUE
   return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, mnew);
 }
 
-
-static int gsl_matrix_complex_equal(const gsl_matrix_complex *m1,
+// Starting with version 1.15, GSL provides a gsl_matrix_complex_equal
+// function, but it only determines absolute equality (i.e. is has no epsilon
+// argument).
+static int gsl_matrix_complex_equal_eps(const gsl_matrix_complex *m1,
   const gsl_matrix_complex *m2, double eps)
 {
   gsl_complex z1, z2;
@@ -1555,7 +1560,7 @@ static VALUE rb_gsl_matrix_complex_equal(int argc, VALUE *argv, VALUE obj)
   Data_Get_Struct(obj, gsl_matrix_complex, m1);
   CHECK_MATRIX_COMPLEX(argv[0]);
   Data_Get_Struct(argv[0], gsl_matrix_complex, m2);
-  ret = gsl_matrix_complex_equal(m1, m2, eps);
+  ret = gsl_matrix_complex_equal_eps(m1, m2, eps);
   if (ret == 1) return Qtrue;
   else return Qfalse;
 }

data/ext/matrix_int.c

@@ -71,13 +71,13 @@ enum {
   GSL_MATRIX_INT_DIV,
 };
 
-static VALUE rb_gsl_matrix_int_operation1(VALUE obj, VALUE other, int flag);
 static VALUE rb_gsl_matrix_int_operation1(VALUE obj, VALUE other, int flag)
 {
   gsl_matrix_int *a, *anew, *b;
   gsl_vector_int *vi, *vinew;
   double bval;
-  int result;
+  // local variable "result" declared and set, but never used
+  //int result;
   Data_Get_Struct(obj, gsl_matrix_int, a);
   switch (TYPE(other)) {
   case T_FIXNUM:
@@ -86,16 +86,16 @@ static VALUE rb_gsl_matrix_int_operation1(VALUE obj, VALUE other, int flag)
     anew = make_matrix_int_clone(a);
     switch (flag) {
     case GSL_MATRIX_INT_ADD: 
-      result = gsl_matrix_int_add_constant(anew, bval);
+      /*result =*/ gsl_matrix_int_add_constant(anew, bval);
       break;
     case GSL_MATRIX_INT_SUB: 
-      result = gsl_matrix_int_add_constant(anew, -bval);
+      /*result =*/ gsl_matrix_int_add_constant(anew, -bval);
       break;
     case GSL_MATRIX_INT_MUL: 
-      result = gsl_matrix_int_scale(anew, bval);
+      /*result =*/ gsl_matrix_int_scale(anew, bval);
       break;
     case GSL_MATRIX_INT_DIV: 
-      result = gsl_matrix_int_scale(anew, 1.0/bval);
+      /*result =*/ gsl_matrix_int_scale(anew, 1.0/bval);
       break;
     default:
       break;
@@ -109,16 +109,16 @@ static VALUE rb_gsl_matrix_int_operation1(VALUE obj, VALUE other, int flag)
       Data_Get_Struct(other, gsl_matrix_int, b);
       switch (flag) {
       case GSL_MATRIX_INT_ADD:
-	result = gsl_matrix_int_add(anew, b);
+	/*result =*/ gsl_matrix_int_add(anew, b);
 	break;
       case GSL_MATRIX_INT_SUB:
-	result = gsl_matrix_int_sub(anew, b);
+	/*result =*/ gsl_matrix_int_sub(anew, b);
 	break;
       case GSL_MATRIX_INT_MUL:
-	result = gsl_matrix_int_mul_elements(anew, b);
+	/*result =*/ gsl_matrix_int_mul_elements(anew, b);
 	break;
       case GSL_MATRIX_INT_DIV:
-	result = gsl_matrix_int_div_elements(anew, b);
+	/*result =*/ gsl_matrix_int_div_elements(anew, b);
 	break;
       default:
 	break;

data/ext/matrix_source.c

@@ -875,14 +875,15 @@ static VALUE FUNCTION(rb_gsl_matrix,to_s)(VALUE obj)
   char buf[32], format[32], format2[32];
   size_t i, j;
   VALUE str;
-  BASE x, min;
+  BASE x;
   int dig = 8;
 #ifdef BASE_INT
+  BASE min;
   BASE max;
 #endif
   Data_Get_Struct(obj, GSL_TYPE(gsl_matrix), m);
-  min = FUNCTION(gsl_matrix,min)(m);
 #ifdef BASE_INT
+  min = FUNCTION(gsl_matrix,min)(m);
   max = gsl_matrix_int_max(m);
   dig = (int) GSL_MAX(fabs(max),fabs(min));
   if (dig > 0) dig = ceil(log10(dig+1e-10));

data/ext/multifit.c

@@ -354,14 +354,15 @@ static VALUE rb_gsl_multifit_fdfsolver_gradient(int argc, VALUE *argv, VALUE obj
 {
   gsl_multifit_fdfsolver *solver = NULL;
   gsl_vector *g = NULL;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   Data_Get_Struct(obj, gsl_multifit_fdfsolver, solver);
   if (argc == 1) {
     Data_Get_Vector(argv[0], g);
     return INT2FIX(gsl_multifit_gradient(solver->J, solver->f, g));
   } else {
     g = gsl_vector_alloc(solver->x->size);
-    status = gsl_multifit_gradient(solver->J, solver->f, g);
+    /*status =*/ gsl_multifit_gradient(solver->J, solver->f, g);
     return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, g);
   }
 }
@@ -371,7 +372,8 @@ static VALUE rb_gsl_multifit_fdfsolver_covar(int argc, VALUE *argv, VALUE obj)
   gsl_multifit_fdfsolver *solver = NULL;
   gsl_matrix *covar = NULL;
   double epsrel;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   if (argc < 1) rb_raise(rb_eArgError, "too few arguments");
   Need_Float(argv[0]);
   Data_Get_Struct(obj, gsl_multifit_fdfsolver, solver);
@@ -379,7 +381,7 @@ static VALUE rb_gsl_multifit_fdfsolver_covar(int argc, VALUE *argv, VALUE obj)
   switch (argc) {
   case 1:
     covar = gsl_matrix_alloc(solver->x->size, solver->x->size);
-    status = gsl_multifit_covar(solver->J, epsrel, covar);
+    /*status =*/ gsl_multifit_covar(solver->J, epsrel, covar);
     return Data_Wrap_Struct(cgsl_matrix, 0, gsl_matrix_free, covar);
     break;
   case 2:
@@ -442,13 +444,14 @@ static VALUE rb_gsl_multifit_gradient(int argc, VALUE *argv, VALUE obj)
 {
   gsl_vector *g = NULL, *f = NULL;
   gsl_matrix *J = NULL;
-  int status;
+  // local variable "status" declared and set, but never used
+  //int status;
   switch (argc) {
   case 2:
     Data_Get_Matrix(argv[0], J);
     Data_Get_Vector(argv[1], f);  
     g = gsl_vector_alloc(f->size);
-    status = gsl_multifit_gradient(J, f, g);
+    /*status =*/ gsl_multifit_gradient(J, f, g);
     return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, g);
     break;  
   case 3:
@@ -474,7 +477,7 @@ static VALUE rb_gsl_multifit_covar(int argc, VALUE *argv, VALUE obj)
     Data_Get_Matrix(argv[0], J);
     epsrel = NUM2DBL(argv[1]);
     covar = gsl_matrix_alloc(J->size2, J->size2);
-    status = gsl_multifit_covar(J, epsrel, covar);
+    /*status =*/ gsl_multifit_covar(J, epsrel, covar);
     return Data_Wrap_Struct(cgsl_matrix, 0, gsl_matrix_free, covar);
     break;
   case 3:
@@ -783,10 +786,9 @@ static int Gaussian_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A, x0, var, xi, yy, wi;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   var = gsl_vector_get(v, 3);
   x0 = gsl_vector_get(v, 2);
@@ -854,10 +856,9 @@ static int Gaussian_2peaks_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A1, x01, var1, A2, x02, var2, xi, yy, wi;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   A1 = gsl_vector_get(v, 1);
   x01 = gsl_vector_get(v, 2);
@@ -925,10 +926,9 @@ static int Exponential_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A, b, xi, yy, wi;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   A = gsl_vector_get(v, 1);
   b = gsl_vector_get(v, 2);
@@ -991,10 +991,9 @@ static int DblExponential_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A1, b1, A2, b2, xi, yy1, yy2, wi;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   A1 = gsl_vector_get(v, 1);
   b1 = gsl_vector_get(v, 2);
@@ -1059,10 +1058,9 @@ static int Lorentzian_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A, B, x0, xi, yy, wi;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   A = gsl_vector_get(v, 1);
   x0 = gsl_vector_get(v, 2);
@@ -1124,10 +1122,9 @@ static int Sin_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A, fc, phi, xi, ys, yc, wi;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   A = gsl_vector_get(v, 1);
   fc = gsl_vector_get(v, 2);
@@ -1190,10 +1187,9 @@ static int Hill_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double y0, m, xhalf, r, yy, xi, wi, a;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   y0 = gsl_vector_get(v, 0);
   m = gsl_vector_get(v, 1);
@@ -1257,10 +1253,9 @@ static int Sigmoid_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double m, x0, r, xi, wi, a, yy;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   m = gsl_vector_get(v, 1);
   x0 = gsl_vector_get(v, 2);
@@ -1322,10 +1317,9 @@ static int Power_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A, r, xi, wi, a;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   A = gsl_vector_get(v, 1);
   r = gsl_vector_get(v, 2);
@@ -1387,10 +1381,9 @@ static int Lognormal_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A, x0, width, xi, wi, a, b;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   A = gsl_vector_get(v, 1);
   x0 = gsl_vector_get(v, 2);
@@ -1449,10 +1442,9 @@ static int Rayleigh_df(const gsl_vector *v, void *data, gsl_matrix *J)
   double A, var, xi, yy, wi;
   size_t i;
   struct fitting_xydata *xydata;
-  gsl_vector *x, *y, *w;
+  gsl_vector *x, *w;
   xydata = (struct fitting_xydata*) data;
   x = xydata->x;
-  y = xydata->y;
   w = xydata->w;
   var = gsl_vector_get(v, 1);
   A = gsl_vector_get(v, 0);