rb-gsl 1.15.3.1 → 1.15.3.2

data/ext/histogram.c

@@ -118,7 +118,8 @@ static VALUE rb_gsl_histogram_alloc_from_file(VALUE klass, VALUE name)
   sprintf(buf, "wc %s", filename);
   fp = popen(buf, "r");
   if (fp == NULL) rb_raise(rb_eIOError, "popen failed.");
-  fgets(buf, 1024, fp);
+  if (fgets(buf, 1024, fp) == NULL)
+    rb_sys_fail(0);
   pclose(fp);
   sscanf(buf, "%d", &nn);
   n = (size_t) nn;      /* vector length */
@@ -355,7 +356,7 @@ static VALUE rb_gsl_histogram_accumulate(int argc, VALUE *argv, VALUE obj)
   Data_Get_Struct(obj, gsl_histogram, h);
   if (TYPE(argv[0]) == T_ARRAY) {
     //    for (i = 0; i < RARRAY(argv[0])->len; i++)
-    for (i = 0; i < RARRAY_LEN(argv[0]); i++)
+    for (i = 0; (int) i < RARRAY_LEN(argv[0]); i++)
       gsl_histogram_accumulate(h, NUM2DBL(rb_ary_entry(argv[0], i)), weight);
   } else if (VECTOR_P(argv[0])) {
     Data_Get_Struct(argv[0], gsl_vector, v);
@@ -1556,13 +1557,13 @@ void mygsl_histogram_integrate(const gsl_histogram *h, gsl_histogram *hi,
 {
   size_t i;
   if (iend >= istart) {
-    if (istart < 0) istart = 0;
+    //if (istart < 0) istart = 0;
     if (iend >= h->n) iend = h->n-1;
     hi->bin[istart] = h->bin[istart];
     for (i = istart+1; i <= iend; i++) hi->bin[i] = hi->bin[i-1] + h->bin[i];
   } else {
     if (istart >= h->n) istart = h->n-1;
-    if (iend < 0) iend = 0;
+    //if (iend < 0) iend = 0;
     hi->bin[istart] = h->bin[istart];
     for (i = istart-1; i >= iend; i--) {
       hi->bin[i] = hi->bin[i+1] + h->bin[i];
@@ -1769,7 +1770,7 @@ static VALUE rb_gsl_histogram_reverse(VALUE obj)
 static double histogram_percentile(const gsl_histogram *h, double f)
 {
   double sum = gsl_histogram_sum(h), sf;
-  double val, s = 0, x;
+  double val = 0, s = 0, x;
   double ri, ri1;
   size_t i;
   sf = sum * f;
@@ -1806,7 +1807,7 @@ static VALUE rb_gsl_histogram_median(VALUE obj)
 static double histogram_percentile_inv(const gsl_histogram *h, double x)
 {
   double sum = gsl_histogram_sum(h);
-  double val, s = 0;
+  double val = 0, s = 0;
   double ri, ri1, q;
   size_t i;
 

data/ext/histogram2d.c

@@ -857,17 +857,17 @@ static gsl_histogram2d* mygsl_histogram2d_calloc_integrate(const gsl_histogram2d
   if (flag == -1) {
     hi->bin[n-1] = h->bin[n-1];
     i = nx - 1;
-    for (j = ny-2, k = 0; j >= 0; j--, k++) {
+    for (j = ny-2, k = 0;; j--, k++) {
       hi->bin[n-1-k] = gsl_histogram2d_get(hi, i, j+1) + gsl_histogram2d_get(h, i, j);
       if (j == 0) break;
     }
     j = ny - 1;
-    for (i = nx-2; i >= 0; i--) {
+    for (i = nx-2;; i--) {
       hi->bin[i*ny + j] = gsl_histogram2d_get(hi, i+1, j) + gsl_histogram2d_get(h, i, j);
       if (i == 0) break;
     }
-    for (i = nx-2; i >= 0; i--) {
-      for (j = ny-2; j >= 0; j--) {
+    for (i = nx-2;; i--) {
+      for (j = ny-2;; j--) {
 	hi->bin[i*ny+j] = ((gsl_histogram2d_get(hi, i+1, j) 
 			    + gsl_histogram2d_get(hi, i, j+1))
 			   - gsl_histogram2d_get(hi, i+1, j+1))

data/ext/interp.c

@@ -220,7 +220,7 @@ static VALUE rb_gsl_interp_evaluate(VALUE obj, VALUE xxa, VALUE yya, VALUE xx,
       ptrz = (double*) na->ptr;
       ary = na_make_object(NA_DFLOAT, na->rank, na->shape, CLASS_OF(xx));
       ptr = NA_PTR_TYPE(ary, double*);
-      for (i = 0; i < na->total; i++)
+      for (i = 0; (int) i < na->total; i++)
 	ptr[i] = (*eval)(rgi->p, ptrx, ptry, ptrz[i], rgi->a);
       return ary;
     }

data/ext/linalg.c

@@ -71,7 +71,7 @@ static VALUE rb_gsl_linalg_LU_decomposition(int argc, VALUE *argv, VALUE obj, in
   Data_Get_Struct(omatrix, gsl_matrix, mtmp);
   if (flag == LINALG_DECOMP_BANG) {
     m = mtmp;
-    RBASIC(omatrix)->klass = cgsl_matrix_LU;
+    RBGSL_SET_CLASS(omatrix, cgsl_matrix_LU);
     objm = omatrix;
   } else {
     m = make_matrix_clone(mtmp);
@@ -692,7 +692,7 @@ static VALUE rb_gsl_linalg_QR_LQ_decomposition(int argc, VALUE *argv, VALUE obj,
     fdecomp = &gsl_linalg_QR_decomp;
     m = mtmp;
     mdecomp = omatrix;
-    RBASIC(mdecomp)->klass = cgsl_matrix_QR;
+    RBGSL_SET_CLASS(mdecomp, cgsl_matrix_QR);
     break;
 #ifdef GSL_1_6_LATER
   case LINALG_LQ_DECOMP:
@@ -704,7 +704,7 @@ static VALUE rb_gsl_linalg_QR_LQ_decomposition(int argc, VALUE *argv, VALUE obj,
     fdecomp = &gsl_linalg_LQ_decomp;
     m = mtmp;
     mdecomp = omatrix;
-    RBASIC(mdecomp)->klass = cgsl_matrix_LQ;
+    RBGSL_SET_CLASS(mdecomp, cgsl_matrix_LQ);
     break;
 #endif
   default:
@@ -731,7 +731,7 @@ static VALUE rb_gsl_linalg_QR_LQ_decomposition(int argc, VALUE *argv, VALUE obj,
       vtau = Data_Wrap_Struct(cgsl_vector_tau, 0, gsl_vector_free, tau);
       return rb_ary_new3(2, mdecomp, vtau);
     } else {
-      RBASIC(argv[itmp])->klass = cgsl_vector_tau;
+      RBGSL_SET_CLASS(argv[itmp], cgsl_vector_tau);
       return mdecomp;
     }
     break;
@@ -740,7 +740,7 @@ static VALUE rb_gsl_linalg_QR_LQ_decomposition(int argc, VALUE *argv, VALUE obj,
    if (argc == itmp) {
       return Data_Wrap_Struct(cgsl_vector_tau, 0, gsl_vector_free, tau);
     } else {
-      RBASIC(argv[itmp])->klass = cgsl_vector_tau;
+      RBGSL_SET_CLASS(argv[itmp], cgsl_vector_tau);
       return INT2FIX(status);
     }
     break;
@@ -1628,14 +1628,14 @@ static VALUE rb_gsl_linalg_QRLQPT_decomp_bang(int argc, VALUE *argv, VALUE obj,
   norm = gsl_vector_alloc(size0);
   switch (flag) {
   case LINALG_QRPT:
-    RBASIC(vA)->klass = cgsl_matrix_QRPT;
+    RBGSL_SET_CLASS(vA, cgsl_matrix_QRPT);
     vtau = Data_Wrap_Struct(cgsl_vector_tau, 0, gsl_vector_free, tau);
     vp = Data_Wrap_Struct(cgsl_permutation, 0, gsl_permutation_free, p);
     gsl_linalg_QRPT_decomp(A, tau, p, &signum, norm);
     break;
 #ifdef GSL_1_6_LATER
   case LINALG_PTLQ:
-    RBASIC(vA)->klass = cgsl_matrix_PTLQ;
+    RBGSL_SET_CLASS(vA, cgsl_matrix_PTLQ);
     vtau = Data_Wrap_Struct(cgsl_vector_tau, 0, gsl_vector_free, tau);
     vp = Data_Wrap_Struct(cgsl_permutation, 0, gsl_permutation_free, p);
     gsl_linalg_PTLQ_decomp(A, tau, p, &signum, norm);
@@ -3454,7 +3454,7 @@ static VALUE rb_gsl_linalg_hessenberg_unpack_accum(int argc, VALUE *argv, VALUE
   gsl_matrix *H = NULL, *V = NULL;
   gsl_vector *tau = NULL;
   size_t i;
-  VALUE val;
+  VALUE val = Qnil;
   switch (argc) {
   case 2:
     /* nothing to do */
@@ -3493,7 +3493,7 @@ static VALUE rb_gsl_linalg_hesstri_decomp(int argc, VALUE *argv, VALUE module)
   gsl_matrix *A = NULL, *B = NULL, *Anew, *Bnew;
   gsl_matrix *U = NULL, *V = NULL;
   gsl_vector *work = NULL;
-  VALUE vH, vR, vU, vV, ary;
+  VALUE vH, vR, vU = Qnil, vV = Qnil, ary;
   int flag = 0;
   switch (argc) {
   case 2:
@@ -3546,7 +3546,7 @@ static VALUE rb_gsl_linalg_hesstri_decomp_bang(int argc, VALUE *argv, VALUE modu
   gsl_matrix *A = NULL, *B = NULL;
   gsl_matrix *U = NULL, *V = NULL;
   gsl_vector *work = NULL;
-  VALUE vH, vR, vU, vV, ary;
+  VALUE vH, vR, vU = Qnil, vV = Qnil, ary;
   int flag = 0;
   switch (argc) {
   case 2:
@@ -3665,6 +3665,9 @@ void Init_gsl_linalg(VALUE module)
   VALUE mgsl_linalg_tridiag;
   VALUE mgsl_linalg_HH;
   VALUE mgsl_linalg_Householder;
+#ifdef GSL_1_9_LATER
+  VALUE mhessen;
+#endif
 
   mgsl_linalg = rb_define_module_under(module, "Linalg");
   mgsl_linalg_LU = rb_define_module_under(mgsl_linalg, "LU");
@@ -4019,7 +4022,6 @@ void Init_gsl_linalg(VALUE module)
 #endif
 
 #ifdef GSL_1_9_LATER
-  VALUE mhessen;
   mhessen = rb_define_module_under(mgsl_linalg, "Hessenberg");
   rb_define_module_function(mhessen, "decomp", rb_gsl_linalg_hessenberg_decomp, 1);
   rb_define_module_function(mgsl_linalg, "heesenberg_decomp", rb_gsl_linalg_hessenberg_decomp, 1);  

data/ext/linalg_complex.c

@@ -51,8 +51,8 @@ VALUE rb_gsl_linalg_complex_LU_decomp(int argc, VALUE *argv, VALUE obj)
   case 0:
     p = gsl_permutation_alloc(size);
     gsl_linalg_complex_LU_decomp(m, p, &signum);
-    if (itmp == 1) RBASIC(argv[0])->klass = cgsl_matrix_complex_LU;
-    else RBASIC(obj)->klass = cgsl_matrix_complex_LU;
+    if (itmp == 1) RBGSL_SET_CLASS(argv[0], cgsl_matrix_complex_LU);
+    else RBGSL_SET_CLASS(obj, cgsl_matrix_complex_LU);
     obj2 = Data_Wrap_Struct(cgsl_permutation, 0, gsl_permutation_free, p);
     return rb_ary_new3(2, obj2, INT2FIX(signum));
     break;
@@ -60,8 +60,8 @@ VALUE rb_gsl_linalg_complex_LU_decomp(int argc, VALUE *argv, VALUE obj)
     CHECK_PERMUTATION(argv[itmp]);
     Data_Get_Struct(argv[itmp], gsl_permutation, p);
     gsl_linalg_complex_LU_decomp(m, p, &signum);
-    if (itmp == 1) RBASIC(argv[0])->klass = cgsl_matrix_complex_LU;
-    else RBASIC(obj)->klass = cgsl_matrix_complex_LU;
+    if (itmp == 1) RBGSL_SET_CLASS(argv[0], cgsl_matrix_complex_LU);
+    else RBGSL_SET_CLASS(obj, cgsl_matrix_complex_LU);
     return INT2FIX(signum);
     break;
   default:
@@ -617,6 +617,9 @@ void Init_gsl_linalg_complex(VALUE module)
 {
   VALUE mgsl_linalg_complex;
   VALUE mgsl_linalg_complex_LU;
+#ifdef GSL_1_10_LATER
+  VALUE mgsl_linalg_complex_chol, mgsl_linalg_complex_Householder;
+#endif
 
   mgsl_linalg_complex = rb_define_module_under(module, "Complex");
   mgsl_linalg_complex_LU = rb_define_module_under(mgsl_linalg_complex, "LU");
@@ -701,7 +704,6 @@ void Init_gsl_linalg_complex(VALUE module)
 #endif
 
 #ifdef GSL_1_10_LATER
-  VALUE mgsl_linalg_complex_chol;
   mgsl_linalg_complex_chol = rb_define_module_under(mgsl_linalg_complex, "Cholesky");
   cgsl_matrix_complex_C = rb_define_class_under(mgsl_linalg_complex_chol, "CholeskyMatrix", cgsl_matrix_complex);    
   rb_define_singleton_method(mgsl_linalg_complex_chol, "decomp", rb_gsl_linalg_cholesky_decomp, -1);
@@ -713,7 +715,7 @@ void Init_gsl_linalg_complex(VALUE module)
   rb_define_method(cgsl_matrix_complex, "cholesky_svx", rb_gsl_linalg_cholesky_svx, -1);
   rb_define_method(cgsl_matrix_complex_C, "svx", rb_gsl_linalg_cholesky_svx, -1);
 
-  VALUE mgsl_linalg_complex_Householder = rb_define_module_under(mgsl_linalg_complex, "Householder");
+  mgsl_linalg_complex_Householder = rb_define_module_under(mgsl_linalg_complex, "Householder");
  rb_define_singleton_method(mgsl_linalg_complex, "householder_transform", 
 			     rb_gsl_linalg_complex_householder_transform, -1);
   rb_define_singleton_method(mgsl_linalg_complex_Householder, "transform", 

data/ext/math.c

@@ -213,7 +213,7 @@ static VALUE rb_gsl_math_eval2(double (*func)(const double, const double), VALUE
     Check_Type(yy, T_ARRAY);
     size = RARRAY_LEN(xx);
     //    if (size != RARRAY(yy)->len) rb_raise(rb_eRuntimeError, "array sizes are different.");
-    if (size != RARRAY_LEN(yy)) rb_raise(rb_eRuntimeError, "array sizes are different.");
+    if ((int) size != RARRAY_LEN(yy)) rb_raise(rb_eRuntimeError, "array sizes are different.");
     ary = rb_ary_new2(size);
     for (i = 0; i < size; i++) {
       x = rb_ary_entry(xx, i);

data/ext/matrix_complex.c

@@ -343,7 +343,7 @@ static VALUE rb_gsl_matrix_complex_set(int argc, VALUE *argv, VALUE obj)
     row_set_argv[0] = INT2FIX(0);
     row_set_argv[1] = INT2FIX(m->size2);
 
-    for(k = 0; k < argc && k < m->size1; k++) {
+    for(k = 0; (int) k < argc && k < m->size1; k++) {
       vv = gsl_matrix_complex_row(m, k);
       rb_gsl_vector_complex_set_subvector(2, row_set_argv, &vv.vector, argv[k]);
     }
@@ -374,7 +374,7 @@ static VALUE rb_gsl_matrix_complex_set(int argc, VALUE *argv, VALUE obj)
         rb_gsl_vector_complex_set_subvector(2, row_set_argv, &vv.vector, other);
       } else {
         // m[...] = [[row0], [row1], ...] # multiple rows
-        if(n1 != RARRAY_LEN(other)) {
+        if((int) n1 != RARRAY_LEN(other)) {
           rb_raise(rb_eRangeError, "row counts do not match (%d != %d)",
 		   (int) n1, (int) RARRAY_LEN(other));
         }
@@ -414,7 +414,7 @@ static VALUE rb_gsl_matrix_complex_set_row(int argc, VALUE *argv, VALUE obj)
   CHECK_FIXNUM(argv[0]);
   Data_Get_Struct(obj, gsl_matrix_complex, A);
   i = FIX2INT(argv[0]);
-  for (k = 1; k < argc; k++) {
+  for (k = 1; (int) k < argc; k++) {
     if (k-1 >= A->size1) break;
     switch (TYPE(argv[k])) {
     case T_ARRAY:
@@ -442,7 +442,7 @@ static VALUE rb_gsl_matrix_complex_set_col(int argc, VALUE *argv, VALUE obj)
   Data_Get_Struct(obj, gsl_matrix_complex, A);
   j = FIX2INT(argv[0]);
   for (k = 1; k < argc; k++) {
-    if (k-1 >= A->size2) break;
+    if (k-1 >= (int) A->size2) break;
     switch (TYPE(argv[k])) {
     case T_ARRAY:
       z = ary2complex(argv[k]);
@@ -632,13 +632,13 @@ static VALUE rb_gsl_matrix_complex_to_s(int argc, VALUE *argv, VALUE obj)
           "%s[ %4.3e %4.3e ]", (j==0) ? "" : " ", GSL_REAL(z), GSL_IMAG(z));
       rb_str_cat(str, buf, strlen(buf));
       // if too many cols
-      if (j >= max_cols-1 && j != m->size2-1) {
+      if ((int) j >= max_cols-1 && j != m->size2-1) {
         rb_str_cat(str, " ...", 4);
         break;
       }
     }
     // if too many rows
-    if (i >= max_rows-1 && i != m->size1-1) {
+    if ((int) i >= max_rows-1 && i != m->size1-1) {
       rb_str_cat(str, "\n  ...", 6);
       break;
     }

data/ext/matrix_source.c

@@ -349,7 +349,7 @@ static GSL_TYPE(gsl_matrix)* FUNCTION(cr_matrix,from_ranges)(int argc, VALUE *ar
   FUNCTION(get_range,beg_en_n)(argv[0], &beg, &en, &n, &step);
   m = FUNCTION(gsl_matrix,calloc)(argc, n);
   FUNCTION(set_ptr_data,by_range)(m->data, n, argv[0]);
-  for (i = 1; i < argc; i++) {
+  for (i = 1; (int) i < argc; i++) {
     if (CLASS_OF(argv[i]) != rb_cRange) 
       rb_raise(rb_eTypeError, "wrong argument type %s (Range expected)",
 	       rb_class2name(CLASS_OF(argv[i])));
@@ -392,11 +392,11 @@ GSL_TYPE(gsl_matrix)* FUNCTION(gsl_matrix,alloc_from_arrays)(int argc, VALUE *ar
   n = RARRAY_LEN(argv[0]);
   m = FUNCTION(gsl_matrix,alloc)(argc, n);
   if (m == NULL) rb_raise(rb_eNoMemError, "gsl_matrix_alloc failed");
-  for (i = 0; i < argc; i++) {
+  for (i = 0; (int) i < argc; i++) {
     if (CLASS_OF(argv[i]) == rb_cRange) argv[i] = rb_gsl_range2ary(argv[i]);
     else Check_Type(argv[i], T_ARRAY);
     for (j = 0; j < n; j++) {
-      if (j >= RARRAY_LEN(argv[i])) FUNCTION(gsl_matrix,set)(m, i, j, 0);
+      if ((int) j >= RARRAY_LEN(argv[i])) FUNCTION(gsl_matrix,set)(m, i, j, 0);
       else FUNCTION(gsl_matrix,set)(m, i, j, NUMCONV2(rb_ary_entry(argv[i], j)));
     }
   }
@@ -413,7 +413,7 @@ GSL_TYPE(gsl_matrix)* FUNCTION(gsl_matrix,alloc_from_vectors)(int argc, VALUE *a
   Data_Get_Struct(argv[0], GSL_TYPE(gsl_vector), v);
   m = FUNCTION(gsl_matrix,alloc)(argc, v->size);
   if (m == NULL) rb_raise(rb_eNoMemError, "gsl_matrix_alloc failed");
-  for (i = 0; i < argc; i++) {
+  for (i = 0; (int) i < argc; i++) {
     CHECK_VEC(argv[i]);
     Data_Get_Struct(argv[i], GSL_TYPE(gsl_vector), v);
     FUNCTION(gsl_matrix,set_row)(m, i, v);
@@ -432,7 +432,7 @@ GSL_TYPE(gsl_matrix)* FUNCTION(gsl_matrix,alloc_from_colvectors)(int argc, VALUE
 //  m = FUNCTION(gsl_matrix,alloc)(argc, v->size);
   m = FUNCTION(gsl_matrix,alloc)(v->size, argc);  
   if (m == NULL) rb_raise(rb_eNoMemError, "gsl_matrix_alloc failed");
-  for (i = 0; i < argc; i++) {
+  for (i = 0; (int) i < argc; i++) {
     CHECK_VEC(argv[i]);
     Data_Get_Struct(argv[i], GSL_TYPE(gsl_vector), v);
     FUNCTION(gsl_matrix,set_col)(m, i, v);
@@ -547,7 +547,7 @@ static VALUE FUNCTION(rb_gsl_matrix,diagonal_singleton)(int argc, VALUE *argv, V
     break;
   default:
     m = FUNCTION(gsl_matrix,calloc)(argc, argc);
-    for (i = 0; i < argc; i++) {
+    for (i = 0; (int) i < argc; i++) {
       FUNCTION(gsl_matrix,set)(m, i, i, NUMCONV2(argv[i]));
     }
     break;
@@ -799,7 +799,7 @@ static VALUE FUNCTION(rb_gsl_matrix,set)(int argc, VALUE *argv, VALUE obj)
         FUNCTION(rb_gsl_vector,set_subvector)(2, row_set_argv, &vv.vector, other);
       } else {
         // m[...] = [[row0], [row1], ...] # multiple rows
-        if(n1 != RARRAY_LEN(other)) {
+        if((int) n1 != RARRAY_LEN(other)) {
           rb_raise(rb_eRangeError, "row counts do not match (%d != %d)",
 		   (int) n1, (int) RARRAY_LEN(other));
         }
@@ -908,7 +908,7 @@ static VALUE FUNCTION(rb_gsl_matrix,to_s)(VALUE obj)
       else
 	sprintf(buf, format2, x); 
       rb_str_cat(str, buf, strlen(buf));
-      if (j >= (55/dig)) {
+      if ((int) j >= (55/dig)) {
         strcpy(buf, "... ");
         rb_str_cat(str, buf, strlen(buf));
         break;
@@ -1062,7 +1062,7 @@ static VALUE FUNCTION(rb_gsl_matrix,set_row)(VALUE obj, VALUE i, VALUE vv)
   if (CLASS_OF(vv) == rb_cRange) vv = rb_gsl_range2ary(vv);
   if (TYPE(vv) == T_ARRAY) {
     v = FUNCTION(gsl_vector,alloc)(RARRAY_LEN(vv));
-    for (j = 0; j < RARRAY_LEN(vv); j++) {
+    for (j = 0; (int) j < RARRAY_LEN(vv); j++) {
       FUNCTION(gsl_vector,set)(v, j, NUMCONV2(rb_ary_entry(vv, j)));
     }
     flag = 1;
@@ -1086,7 +1086,7 @@ static VALUE FUNCTION(rb_gsl_matrix,set_col)(VALUE obj, VALUE j, VALUE vv)
   if (CLASS_OF(vv) == rb_cRange) vv = rb_gsl_range2ary(vv);
   if (TYPE(vv) == T_ARRAY) {
     v = FUNCTION(gsl_vector,alloc)(RARRAY_LEN(vv));
-    for (i = 0; i < RARRAY_LEN(vv); i++) {
+    for (i = 0; (int) i < RARRAY_LEN(vv); i++) {
       FUNCTION(gsl_vector,set)(v, i, NUMCONV2(rb_ary_entry(vv, i)));
     }
     flag = 1;

data/ext/monte.c

@@ -72,7 +72,7 @@ static VALUE rb_gsl_monte_function_set_f(int argc, VALUE *argv, VALUE obj)
     rb_ary_store(ary, 0, argv[0]);
     F->dim = FIX2INT(argv[1]);
     ary2 = rb_ary_new2(argc-2);
-    for (i = 2; i < argc; i++) rb_ary_store(ary2, i-2, argv[i]);
+    for (i = 2; i < (size_t) argc; i++) rb_ary_store(ary2, i-2, argv[i]);
     rb_ary_store(ary, 1, ary2);
     break;
   }
@@ -164,7 +164,7 @@ static VALUE rb_gsl_monte_function_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; i < (size_t) argc; i++) rb_ary_store(ary2, i, argv[i]);
     rb_ary_store(ary, 1, ary2);
   }
   return obj;

data/ext/multimin.c

@@ -99,7 +99,7 @@ static VALUE rb_gsl_multimin_function_new(int argc, VALUE *argv, VALUE klass)
     break;
   case 2:
   case 3:
-    for (i = 0; i < argc; i++) set_function(i, argv, F);
+    for (i = 0; (int) i < argc; i++) set_function(i, argv, F);
     break;
   default:
     rb_raise(rb_eArgError, "wrong number of arguments");
@@ -172,7 +172,7 @@ static VALUE rb_gsl_multimin_function_set_f(int argc, VALUE *argv, VALUE obj)
     break;
   case 2:
   case 3:
-    for (i = 0; i < argc; i++) set_function(i, argv, F);
+    for (i = 0; (int) i < argc; i++) set_function(i, argv, F);
     break;
   default:
     rb_raise(rb_eArgError, "wrong number of arguments");
@@ -197,7 +197,7 @@ static VALUE rb_gsl_multimin_function_set_params(int argc, VALUE *argv, VALUE ob
   if (argc == 1) 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;
@@ -358,7 +358,7 @@ static VALUE rb_gsl_multimin_function_fdf_set_params(int argc, VALUE *argv, VALU
 		rb_ary_store(ary, 3, 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, 3, ary2);
   }
   return obj;

data/ext/multiroots.c

@@ -80,7 +80,7 @@ static VALUE rb_gsl_multiroot_function_new(int argc, VALUE *argv, VALUE klass)
     break;
   case 2:
   case 3:
-    for (i = 0; i < argc; i++) set_function(i, argv, F);
+    for (i = 0; (int) i < argc; i++) set_function(i, argv, F);
     break;
   default:
     rb_raise(rb_eArgError, "wrong number of arguments");
@@ -99,7 +99,7 @@ static void gsl_multiroot_function_mark(gsl_multiroot_function *f)
   size_t i;
   rb_gc_mark((VALUE) f->params);
   //  for (i = 0; i < RARRAY(f->params)->len; i++) 
-  for (i = 0; i < RARRAY_LEN(f->params); i++) 
+  for (i = 0; (int) i < RARRAY_LEN(f->params); i++) 
     rb_gc_mark(rb_ary_entry((VALUE) f->params, i));
 }
 
@@ -161,7 +161,7 @@ static VALUE rb_gsl_multiroot_function_set_f(int argc, VALUE *argv, VALUE obj)
     break;
   case 2:
   case 3:
-    for (i = 0; i < argc; i++) set_function(i, argv, F);
+    for (i = 0; (int) i < argc; i++) set_function(i, argv, F);
     break;
   default:
     rb_raise(rb_eArgError, "wrong number of arguments");
@@ -187,7 +187,7 @@ static VALUE rb_gsl_multiroot_function_set_params(int argc, VALUE *argv, VALUE o
   if (argc == 1) 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;
@@ -236,7 +236,7 @@ static void gsl_multiroot_function_fdf_mark(gsl_multiroot_function_fdf *f)
   size_t i;
   rb_gc_mark((VALUE) f->params);
   //  for (i = 0; i < RARRAY(f->params)->len; i++) 
-  for (i = 0; i < RARRAY_LEN(f->params); i++) 
+  for (i = 0; (int) i < RARRAY_LEN(f->params); i++) 
     rb_gc_mark(rb_ary_entry((VALUE) f->params, i));
 }
 
@@ -322,7 +322,7 @@ static VALUE rb_gsl_multiroot_function_fdf_set_params(int argc, VALUE *argv, VAL
   if (argc == 1) rb_ary_store(ary, 3, 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, 3, ary2);
   }
   return obj;
@@ -813,7 +813,7 @@ static VALUE rb_gsl_multiroot_function_solve(int argc, VALUE *argv, VALUE obj)
   case 4:
   case 3:
   case 2:
-    for (i = 1; i < argc; i++) {
+    for (i = 1; (int) i < argc; i++) {
       switch (TYPE(argv[i])) {
       case T_STRING:
 	T = (gsl_multiroot_fsolver_type *) get_fsolver_type(argv[i]);
@@ -830,7 +830,7 @@ static VALUE rb_gsl_multiroot_function_solve(int argc, VALUE *argv, VALUE obj)
   case 1:
     if (TYPE(argv[0]) == T_ARRAY) {
       //      if (RARRAY(argv[0])->len != F->n)
-      if (RARRAY_LEN(argv[0]) != F->n)
+      if (RARRAY_LEN(argv[0]) != (int) F->n)
 	rb_raise(rb_eRangeError, "array size are different.");
       x0 = gsl_vector_alloc(F->n);
       for (i = 0; i < x0->size; i++)