gsl 1.15.3 → 1.16.0.6

data/ext/{matrix.c → gsl_native/matrix.c}

@@ -14,26 +14,23 @@
   which are provided by the GSL source.
  */
 
-#include "rb_gsl_config.h"
-#include "rb_gsl_array.h"
-#include "rb_gsl_histogram.h"
-#include "rb_gsl_complex.h"
-#include "rb_gsl_poly.h"
-#ifdef HAVE_NARRAY_H
-#include "rb_gsl_with_narray.h"
-#endif
+#include "include/rb_gsl_array.h"
+#include "include/rb_gsl_histogram.h"
+#include "include/rb_gsl_complex.h"
+#include "include/rb_gsl_poly.h"
+#include "include/rb_gsl_with_narray.h"
 
-int gsl_linalg_matmult_int(const gsl_matrix_int *A, 
-			   const gsl_matrix_int *B, gsl_matrix_int *C);
+int gsl_linalg_matmult_int(const gsl_matrix_int *A,
+                           const gsl_matrix_int *B, gsl_matrix_int *C);
 
 #define BASE_DOUBLE
-#include "templates_on.h"
-#include "matrix_source.c"
-#include "templates_off.h"
+#include "include/templates_on.h"
+#include "matrix_source.h"
+#include "include/templates_off.h"
 #undef  BASE_DOUBLE
 
 #define BASE_INT
-#include "templates_on.h"
-#include "matrix_source.c"
-#include "templates_off.h"
+#include "include/templates_on.h"
+#include "matrix_source.h"
+#include "include/templates_off.h"
 #undef  BASE_INT

data/ext/{matrix_complex.c → gsl_native/matrix_complex.c}

@@ -9,9 +9,8 @@
   WITHOUT ANY WARRANTY.
 */
 
-#include "rb_gsl_config.h"
-#include "rb_gsl_array.h"
-#include "rb_gsl_complex.h"
+#include "include/rb_gsl_array.h"
+#include "include/rb_gsl_complex.h"
 
 enum {
   GSL_MATRIX_COMPLEX_ADD,
@@ -28,7 +27,7 @@ void get_range_beg_en_n(VALUE range, double *beg, double *en, size_t *n, int *st
 
 // From ext/matrix_source.c
 void parse_submatrix_args(int argc, VALUE *argv, size_t size1, size_t size2,
-    size_t *i, size_t *j, size_t *n1, size_t *n2);
+                          size_t *i, size_t *j, size_t *n1, size_t *n2);
 
 static VALUE rb_gsl_matrix_complex_arithmetics(int flag, VALUE obj, VALUE bb)
 {
@@ -74,28 +73,28 @@ static VALUE rb_gsl_matrix_complex_arithmetics(int flag, VALUE obj, VALUE bb)
       Data_Get_Struct(bb, gsl_matrix_complex, cmb);
       switch (flag) {
       case GSL_MATRIX_COMPLEX_ADD:
-	cmnew = make_matrix_complex_clone(cm);
-	gsl_matrix_complex_add(cmnew, cmb);
-	return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
-	break;
+        cmnew = make_matrix_complex_clone(cm);
+        gsl_matrix_complex_add(cmnew, cmb);
+        return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
+        break;
       case GSL_MATRIX_COMPLEX_SUB:
-	cmnew = make_matrix_complex_clone(cm);
-	gsl_matrix_complex_sub(cmnew,cmb);
-	return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
-	break;
+        cmnew = make_matrix_complex_clone(cm);
+        gsl_matrix_complex_sub(cmnew,cmb);
+        return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
+        break;
       case GSL_MATRIX_COMPLEX_MUL:
-	cmnew = make_matrix_complex_clone(cm);
-	gsl_matrix_complex_mul_elements(cmnew, cmb);
-	return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
-	break;
+        cmnew = make_matrix_complex_clone(cm);
+        gsl_matrix_complex_mul_elements(cmnew, cmb);
+        return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
+        break;
       case GSL_MATRIX_COMPLEX_DIV:
-	cmnew = make_matrix_complex_clone(cm);
-	gsl_matrix_complex_div_elements(cmnew, cmb);
-	return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
-	break;
+        cmnew = make_matrix_complex_clone(cm);
+        gsl_matrix_complex_div_elements(cmnew, cmb);
+        return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
+        break;
       default:
-	rb_raise(rb_eRuntimeError, "operation not defined");
-	break;
+        rb_raise(rb_eRuntimeError, "operation not defined");
+        break;
       }
     } else if (rb_obj_is_kind_of(bb, cgsl_matrix)) {
       Data_Get_Struct(bb, gsl_matrix, m);
@@ -103,80 +102,80 @@ static VALUE rb_gsl_matrix_complex_arithmetics(int flag, VALUE obj, VALUE bb)
       cmnew = make_matrix_complex_clone(cm);
       switch (flag) {
       case GSL_MATRIX_COMPLEX_ADD:
-	gsl_matrix_complex_add(cmnew, cmb);
-	break;
+        gsl_matrix_complex_add(cmnew, cmb);
+        break;
       case GSL_MATRIX_COMPLEX_SUB:
-	gsl_matrix_complex_sub(cmnew,cmb);
-	break;
+        gsl_matrix_complex_sub(cmnew,cmb);
+        break;
       case GSL_MATRIX_COMPLEX_MUL:
-	gsl_matrix_complex_mul_elements(cmnew, cmb);
-	break;
+        gsl_matrix_complex_mul_elements(cmnew, cmb);
+        break;
       case GSL_MATRIX_COMPLEX_DIV:
-	gsl_matrix_complex_div_elements(cmnew, cmb);
-	break;
+        gsl_matrix_complex_div_elements(cmnew, cmb);
+        break;
       default:
-	rb_raise(rb_eRuntimeError, "operation not defined");
-	break;
+        rb_raise(rb_eRuntimeError, "operation not defined");
+        break;
       }
       return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
     } else if (rb_obj_is_kind_of(bb, cgsl_complex)) {
       Data_Get_Struct(bb, gsl_complex, c);
       switch (flag) {
       case GSL_MATRIX_COMPLEX_ADD:
-	cmnew = make_matrix_complex_clone(cm);
-	gsl_matrix_complex_add_constant(cmnew, *c);
-	break;
+        cmnew = make_matrix_complex_clone(cm);
+        gsl_matrix_complex_add_constant(cmnew, *c);
+        break;
       case GSL_MATRIX_COMPLEX_SUB:
-	cmnew = make_matrix_complex_clone(cm);
-	gsl_matrix_complex_add_constant(cmnew, gsl_complex_negative(*c));
-	break;
+        cmnew = make_matrix_complex_clone(cm);
+        gsl_matrix_complex_add_constant(cmnew, gsl_complex_negative(*c));
+        break;
       case GSL_MATRIX_COMPLEX_MUL:
-	cmnew = make_matrix_complex_clone(cm);
-	gsl_matrix_complex_scale(cmnew, *c);
-	break;
+        cmnew = make_matrix_complex_clone(cm);
+        gsl_matrix_complex_scale(cmnew, *c);
+        break;
       case GSL_MATRIX_COMPLEX_DIV:
-	cmnew = make_matrix_complex_clone(cm);
-	gsl_matrix_complex_scale(cmnew, gsl_complex_inverse(*c));
-	break;
+        cmnew = make_matrix_complex_clone(cm);
+        gsl_matrix_complex_scale(cmnew, gsl_complex_inverse(*c));
+        break;
       default:
-	rb_raise(rb_eRuntimeError, "operation not defined");
-	break;
+        rb_raise(rb_eRuntimeError, "operation not defined");
+        break;
       }
       return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmnew);
     } else if (rb_obj_is_kind_of(bb, cgsl_vector)) {
       Data_Get_Struct(bb, gsl_vector, v);
       switch (flag) {
       case GSL_MATRIX_COMPLEX_MUL:
-	cvb = vector_to_complex(v);
-	cvnew = gsl_vector_complex_alloc(v->size);
-	if (cvnew == NULL) rb_raise(rb_eNoMemError, "gsl_vector_complex_alloc failed");
-	gsl_matrix_complex_mul_vector(cvnew, cm, cvb);
-	gsl_vector_complex_free(cvb);
-	return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, cvnew);
-	break;
+        cvb = vector_to_complex(v);
+        cvnew = gsl_vector_complex_alloc(v->size);
+        if (cvnew == NULL) rb_raise(rb_eNoMemError, "gsl_vector_complex_alloc failed");
+        gsl_matrix_complex_mul_vector(cvnew, cm, cvb);
+        gsl_vector_complex_free(cvb);
+        return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, cvnew);
+        break;
       default:
-	rb_raise(rb_eRuntimeError, 
-		 "operation is not defined %s and Matrix_Complex",
-		 rb_class2name(CLASS_OF(bb)));
-	break;
+        rb_raise(rb_eRuntimeError,
+                 "operation is not defined %s and Matrix_Complex",
+                 rb_class2name(CLASS_OF(bb)));
+        break;
       }
     } else if (rb_obj_is_kind_of(bb, cgsl_vector_complex)) {
       if (!VECTOR_COMPLEX_COL_P(bb))
-	rb_raise(rb_eTypeError, 
-		 "Operation is not defined with %s (Vector::Complex::Col expected)", 
-		 rb_class2name(CLASS_OF(bb)));
+        rb_raise(rb_eTypeError,
+                 "Operation is not defined with %s (Vector::Complex::Col expected)",
+                 rb_class2name(CLASS_OF(bb)));
       Data_Get_Struct(bb, gsl_vector_complex, cv);
       switch (flag) {
       case GSL_MATRIX_COMPLEX_MUL:
-	cvnew = gsl_vector_complex_alloc(cv->size);
-	if (cvnew == NULL) rb_raise(rb_eNoMemError, "gsl_vector_complex_alloc failed");
-	gsl_matrix_complex_mul_vector(cvnew, cm, cv);
-	return Data_Wrap_Struct(cgsl_vector_complex_col, 0, gsl_vector_complex_free, cvnew);;
+        cvnew = gsl_vector_complex_alloc(cv->size);
+        if (cvnew == NULL) rb_raise(rb_eNoMemError, "gsl_vector_complex_alloc failed");
+        gsl_matrix_complex_mul_vector(cvnew, cm, cv);
+        return Data_Wrap_Struct(cgsl_vector_complex_col, 0, gsl_vector_complex_free, cvnew);
       default:
-	rb_raise(rb_eRuntimeError, 
-		 "operation is not defined %s and Matrix_Complex",
-		 rb_class2name(CLASS_OF(bb)));
-	break;
+        rb_raise(rb_eRuntimeError,
+                 "operation is not defined %s and Matrix_Complex",
+                 rb_class2name(CLASS_OF(bb)));
+        break;
       }
     } else {
       rb_raise(rb_eTypeError, "wrong argument type %s", rb_class2name(CLASS_OF(bb)));
@@ -219,15 +218,15 @@ static VALUE rb_gsl_matrix_complex_eye(int argc, VALUE *argv, VALUE klass)
       //      if (RARRAY(argv[1])->len < 2) rb_raise(rb_eArgError, "wrong argument");
       if (RARRAY_LEN(argv[1]) < 2) rb_raise(rb_eArgError, "wrong argument");
       z = gsl_complex_rect(NUM2DBL(rb_ary_entry(argv[1], 0)),
-			   NUM2DBL(rb_ary_entry(argv[1], 1)));
+                           NUM2DBL(rb_ary_entry(argv[1], 1)));
       break;
     default:
       if (rb_obj_is_kind_of(argv[1], cgsl_complex)) {
-	Data_Get_Struct(argv[1], gsl_complex, pz);
+        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])));
+        rb_raise(rb_eTypeError,
+                 "wrong argument type %s", rb_class2name(CLASS_OF(argv[1])));
       }
       break;
     }
@@ -245,7 +244,7 @@ static VALUE rb_gsl_matrix_complex_eye(int argc, VALUE *argv, VALUE klass)
   m = gsl_matrix_complex_calloc(n, n);
   if (m == NULL) rb_raise(rb_eNoMemError, "gsl_matrix_complex_alloc failed");
   for (i = 0; i < n; i++) gsl_matrix_complex_set(m, i, i, z);
-  return  Data_Wrap_Struct(klass, 0, gsl_matrix_complex_free, m);
+  return Data_Wrap_Struct(klass, 0, gsl_matrix_complex_free, m);
 }
 
 static VALUE rb_gsl_matrix_complex_identity(VALUE klass, VALUE nn)
@@ -259,7 +258,7 @@ static VALUE rb_gsl_matrix_complex_identity(VALUE klass, VALUE nn)
   if (m == NULL) rb_raise(rb_eNoMemError, "gsl_matrix_complex_calloc failed");
   z = gsl_complex_rect(1.0, 0.0);
   for (i = 0; i < n; i++) gsl_matrix_complex_set(m, i, i, z);
-  return  Data_Wrap_Struct(klass, 0, gsl_matrix_complex_free, m);
+  return Data_Wrap_Struct(klass, 0, gsl_matrix_complex_free, m);
 }
 
 static VALUE rb_gsl_matrix_complex_set_zero(VALUE obj)
@@ -301,8 +300,8 @@ static VALUE rb_gsl_matrix_complex_set_all(VALUE obj, VALUE s)
       Data_Get_Struct(s, gsl_complex, z);
       gsl_matrix_complex_set_all(m, *z);
     } else {
-      rb_raise(rb_eTypeError, 
-	       "wrong argument type %s", rb_class2name(CLASS_OF(s)));
+      rb_raise(rb_eTypeError,
+               "wrong argument type %s", rb_class2name(CLASS_OF(s)));
     }
     break;
   }
@@ -324,7 +323,6 @@ static VALUE rb_gsl_matrix_complex_set(int argc, VALUE *argv, VALUE obj)
   if(argc < 1 || argc > 5) {
     rb_raise(rb_eArgError, "wrong number of arguments (%d for 1-5)", argc);
   }
-
   Data_Get_Struct(obj, gsl_matrix_complex, m);
   other = argv[argc-1];
 
@@ -343,7 +341,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]);
     }
@@ -358,7 +356,7 @@ static VALUE rb_gsl_matrix_complex_set(int argc, VALUE *argv, VALUE obj)
       Data_Get_Struct(other, gsl_matrix_complex, mother);
       if(n1 * n2 != mother->size1 * mother->size2) {
         rb_raise(rb_eRangeError, "sizes do not match (%d x %d != %d x %d)",
-		 (int) n1, (int) n2, (int) mother->size1, (int) mother->size2);
+                 (int) n1, (int) n2, (int) mother->size1, (int) mother->size2);
       }
       // TODO Change to gsl_matrix_memmove if/when GSL has such a function
       // because gsl_matrix_memcpy does not handle overlapping regions (e.g.
@@ -374,9 +372,9 @@ 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));
+                   (int) n1, (int) RARRAY_LEN(other));
         }
         for(k = 0; k < n1; k++) {
           vv = gsl_matrix_complex_row(&mv.matrix, k);
@@ -400,7 +398,6 @@ static VALUE rb_gsl_matrix_complex_set(int argc, VALUE *argv, VALUE obj)
       gsl_matrix_complex_set_all(&mv.matrix, rb_gsl_obj_to_gsl_complex(other, NULL));
     }
   }
-
   return obj;
 }
 
@@ -410,11 +407,11 @@ static VALUE rb_gsl_matrix_complex_set_row(int argc, VALUE *argv, VALUE obj)
   size_t i, k;
   gsl_complex z, *pz = &z;
   if (argc < 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 2)",
-			 argc);
+                         argc);
   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:
@@ -437,12 +434,12 @@ static VALUE rb_gsl_matrix_complex_set_col(int argc, VALUE *argv, VALUE obj)
   int j, k;
   gsl_complex z, *pz = &z;
   if (argc < 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 2)",
-			 argc);
+                         argc);
   CHECK_FIXNUM(argv[0]);
   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]);
@@ -474,7 +471,7 @@ static VALUE rb_gsl_matrix_complex_get(int argc, VALUE *argv, VALUE obj)
     if(ii < 0) ii += m->size1;
     if(ij < 0) ij += m->size2;
     c = ALLOC(gsl_complex);
-    *c = gsl_matrix_complex_get(m, (size_t)ii, (size_t)ij); 
+    *c = gsl_matrix_complex_get(m, (size_t)ii, (size_t)ij);
     retval = Data_Wrap_Struct(cgsl_complex, 0, free, c);
   } else if(argc == 1 && TYPE(argv[0]) == T_FIXNUM) {
     // m[i]
@@ -493,7 +490,7 @@ static VALUE rb_gsl_matrix_complex_get(int argc, VALUE *argv, VALUE obj)
       if(ii < 0) ii += m->size1;
       if(ij < 0) ij += m->size2;
       c = ALLOC(gsl_complex);
-      *c = gsl_matrix_complex_get(m, (size_t)ii, (size_t)ij); 
+      *c = gsl_matrix_complex_get(m, (size_t)ii, (size_t)ij);
       retval = Data_Wrap_Struct(cgsl_complex, 0, free, c);
     } else {
       rb_raise(rb_eArgError, "Array index must have length 2, not %d", (int) RARRAY_LEN(argv[0]));
@@ -501,7 +498,6 @@ static VALUE rb_gsl_matrix_complex_get(int argc, VALUE *argv, VALUE obj)
   } else {
     retval = rb_gsl_matrix_complex_submatrix(argc, argv, obj);
   }
-
   return retval;
 }
 
@@ -547,12 +543,12 @@ static VALUE rb_gsl_matrix_complex_to_a(VALUE obj)
   size_t i, j;
   Data_Get_Struct(obj, gsl_matrix_complex, m);
   ma = rb_ary_new2(m->size1);
-  for(i=0; i < m->size1; i++) {
+  for(i = 0; i < m->size1; i++) {
     ra = rb_ary_new2(m->size2);
     rb_ary_store(ma, i, ra);
-    for(j=0; j < m->size2; j++) {
+    for(j = 0; j < m->size2; j++) {
       c = ALLOC(gsl_complex);
-      *c = gsl_matrix_complex_get(m, i, j); 
+      *c = gsl_matrix_complex_get(m, i, j);
       rb_ary_store(ra, j, Data_Wrap_Struct(cgsl_complex, 0, free, c));
     }
   }
@@ -611,11 +607,11 @@ static VALUE rb_gsl_matrix_complex_to_s(int argc, VALUE *argv, VALUE obj)
   int max_rows = 4;
   int max_cols = 4;
 
-  switch(argc){
-    case 2: max_cols = NUM2INT(argv[1]);
-    case 1: max_rows = NUM2INT(argv[0]);
-    case 0: break;
-    default:
+  switch(argc) {
+  case 2: max_cols = NUM2INT(argv[1]);
+  case 1: max_rows = NUM2INT(argv[0]);
+  case 0: break;
+  default:
     rb_raise(rb_eArgError, "wrong number of arguments (%d for 0, 1, or 2)", argc);
   }
 
@@ -629,16 +625,16 @@ static VALUE rb_gsl_matrix_complex_to_s(int argc, VALUE *argv, VALUE obj)
     for (j = 0; j < m->size2; j++) {
       z = gsl_matrix_complex_get(m, i, j);
       sprintf(buf,
-          "%s[ %4.3e %4.3e ]", (j==0) ? "" : " ", GSL_REAL(z), GSL_IMAG(z));
+              "%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;
     }
@@ -927,8 +923,8 @@ static VALUE rb_gsl_matrix_complex_add_diagonal(VALUE obj, VALUE s)
       Data_Get_Struct(s, gsl_complex, z);
       gsl_matrix_complex_add_diagonal(m, *z);
     } else {
-      rb_raise(rb_eTypeError, 
-	       "wrong argument type %s", rb_class2name(CLASS_OF(s)));
+      rb_raise(rb_eTypeError,
+               "wrong argument type %s", rb_class2name(CLASS_OF(s)));
     }
     break;
   }
@@ -1002,7 +998,7 @@ static VALUE rb_gsl_matrix_complex_set_diagonal(VALUE obj, VALUE diag)
     for (i = 0; i < m->size1; i++) gsl_matrix_complex_set(m, i, i, gsl_vector_complex_get(v, i));
   } else {
     rb_raise(rb_eTypeError, "wrong argument type %s (GSL::Vector_Complex or Array expected)",
-	     rb_class2name(CLASS_OF(diag)));
+             rb_class2name(CLASS_OF(diag)));
   }
   return obj;
 }
@@ -1056,7 +1052,7 @@ static VALUE rb_gsl_matrix_complex_coerce(VALUE obj, VALUE other)
       return rb_ary_new3(2, vcm, obj);
     } else {
       rb_raise(rb_eTypeError, "cannot coerce %s to GSL::Matrix::Complex",
-	       rb_class2name(CLASS_OF(other)));
+               rb_class2name(CLASS_OF(other)));
     }
     break;
   }
@@ -1274,8 +1270,8 @@ static VALUE rb_gsl_matrix_complex_logabs(VALUE obj)
   return rb_gsl_matrix_complex_XXX(obj, gsl_complex_logabs);
 }
 
-static VALUE rb_gsl_matrix_complex_XXX_complex(VALUE obj, 
-					       gsl_complex (*f)(gsl_complex))
+static VALUE rb_gsl_matrix_complex_XXX_complex(VALUE obj,
+                                               gsl_complex (*f)(gsl_complex))
 {
   gsl_matrix_complex *m, *mnew;
   size_t i, j;
@@ -1526,14 +1522,14 @@ static VALUE rb_gsl_matrix_complex_indgen_singleton(int argc, VALUE *argv, VALUE
 // 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)
+                                        const gsl_matrix_complex *m2, double eps)
 {
   gsl_complex z1, z2;
   size_t i, j;
   if (m1->size1 != m2->size1) return 0;
-  if (m1->size2 != m2->size2) return 0;  
+  if (m1->size2 != m2->size2) return 0;
   for (i = 0; i < m1->size1; i++) {
-    for (j = 0; j < m1->size2; j++) {    
+    for (j = 0; j < m1->size2; j++) {
       z1 = gsl_matrix_complex_get(m1, i, j);
       z2 = gsl_matrix_complex_get(m2, i, j);
       if (!rbgsl_complex_equal(&z1, &z2, eps)) return 0;
@@ -1603,7 +1599,7 @@ void Init_gsl_matrix_complex(VALUE module)
   rb_define_method(cgsl_matrix_complex, "fwrite", rb_gsl_matrix_complex_fwrite, 1);
   rb_define_method(cgsl_matrix_complex, "fread", rb_gsl_matrix_complex_fread, 1);
   rb_define_method(cgsl_matrix_complex, "fscanf", rb_gsl_matrix_complex_fscanf, 1);
-  
+
   rb_define_singleton_method(cgsl_matrix_complex, "memcpy", rb_gsl_matrix_complex_memcpy, 2);
   rb_define_method(cgsl_matrix_complex, "clone", rb_gsl_matrix_complex_clone, 0);
   rb_define_alias(cgsl_matrix_complex, "duplicate", "clone");
@@ -1611,10 +1607,10 @@ void Init_gsl_matrix_complex(VALUE module)
   rb_define_method(cgsl_matrix_complex, "swap_rows", rb_gsl_matrix_complex_swap_rows, 2);
   rb_define_method(cgsl_matrix_complex, "swap_columns", rb_gsl_matrix_complex_swap_columns, 2);
   rb_define_method(cgsl_matrix_complex, "swap_rowcol", rb_gsl_matrix_complex_swap_rowcol, 2);
-  
+
   rb_define_method(cgsl_matrix_complex, "transpose", rb_gsl_matrix_complex_transpose, 0);
   rb_define_method(cgsl_matrix_complex, "isnull", rb_gsl_matrix_complex_isnull, 0);
-  
+
   rb_define_method(cgsl_matrix_complex, "add", rb_gsl_matrix_complex_add, 1);
   rb_define_alias(cgsl_matrix_complex, "add_constant", "add");
   rb_define_alias(cgsl_matrix_complex, "+", "add");
@@ -1625,13 +1621,13 @@ void Init_gsl_matrix_complex(VALUE module)
   rb_define_alias(cgsl_matrix_complex, "/", "div_elements");
   rb_define_method(cgsl_matrix_complex, "scale", rb_gsl_matrix_complex_scale, 1);
   rb_define_method(cgsl_matrix_complex, "scale!", rb_gsl_matrix_complex_scale_bang, 1);
-  
+
   rb_define_method(cgsl_matrix_complex, "add_diagonal", rb_gsl_matrix_complex_add_diagonal, 1);
-  
+
   rb_define_method(cgsl_matrix_complex, "set_zero", rb_gsl_matrix_complex_set_zero, 0);
   rb_define_method(cgsl_matrix_complex, "set_identity", rb_gsl_matrix_complex_set_identity, 0);
   rb_define_method(cgsl_matrix_complex, "set_all", rb_gsl_matrix_complex_set_all, 1);
-  
+
   rb_define_method(cgsl_matrix_complex, "submatrix", rb_gsl_matrix_complex_submatrix, -1);
   rb_define_alias(cgsl_matrix_complex, "view", "submatrix");
   rb_define_method(cgsl_matrix_complex, "row", rb_gsl_matrix_complex_row, 1);
@@ -1643,13 +1639,13 @@ void Init_gsl_matrix_complex(VALUE module)
   rb_define_method(cgsl_matrix_complex, "set_diagonal", rb_gsl_matrix_complex_set_diagonal, 1);
   rb_define_method(cgsl_matrix_complex, "subdiagonal", rb_gsl_matrix_complex_subdiagonal, 1);
   rb_define_method(cgsl_matrix_complex, "superdiagonal", rb_gsl_matrix_complex_superdiagonal, 1);
-  
+
   rb_define_method(cgsl_matrix_complex, "coerce", rb_gsl_matrix_complex_coerce, 1);
-  
+
   rb_define_method(cgsl_matrix_complex, "mul", rb_gsl_matrix_complex_mul, 1);
   rb_define_alias(cgsl_matrix_complex, "*", "mul");
   rb_define_method(cgsl_matrix_complex, "mul!", rb_gsl_matrix_complex_mul2, 1);
-  
+
   rb_define_method(cgsl_matrix_complex, "real", rb_gsl_matrix_complex_real, 0);
   rb_define_alias(cgsl_matrix_complex, "to_real", "real");
   rb_define_alias(cgsl_matrix_complex, "re", "real");
@@ -1663,7 +1659,7 @@ void Init_gsl_matrix_complex(VALUE module)
   rb_define_alias(cgsl_matrix_complex, "conj", "conjugate");
   rb_define_method(cgsl_matrix_complex, "dagger!", rb_gsl_matrix_complex_dagger, 0);
   rb_define_method(cgsl_matrix_complex, "dagger", rb_gsl_matrix_complex_dagger2, 0);
-  
+
   rb_define_method(cgsl_matrix_complex, "trace", rb_gsl_matrix_complex_trace, 0);
   rb_define_method(cgsl_matrix_complex, "each_row", rb_gsl_matrix_complex_each_row, 0);
   rb_define_method(cgsl_matrix_complex, "each_col", rb_gsl_matrix_complex_each_col, 0);
@@ -1672,14 +1668,14 @@ void Init_gsl_matrix_complex(VALUE module)
   rb_define_method(cgsl_matrix_complex, "collect!", rb_gsl_matrix_complex_collect_bang, 0);
   rb_define_alias(cgsl_matrix_complex, "map", "collect");
   rb_define_alias(cgsl_matrix_complex, "map!", "collect!");
-  
+
   rb_define_method(cgsl_matrix_complex, "to_a", rb_gsl_matrix_complex_to_a, 0);
 
   rb_define_method(cgsl_matrix_complex, "size1", rb_gsl_matrix_complex_size1, 0);
   rb_define_method(cgsl_matrix_complex, "size2", rb_gsl_matrix_complex_size2, 0);
   rb_define_method(cgsl_matrix_complex, "shape", rb_gsl_matrix_complex_shape, 0);
   rb_define_alias(cgsl_matrix_complex, "size", "shape");
-  
+
   /*****/
   rb_define_method(cgsl_matrix_complex, "-@", rb_gsl_matrix_complex_uminus, 0);
   rb_define_method(cgsl_matrix_complex, "+@", rb_gsl_matrix_complex_uplus, 0);
@@ -1729,9 +1725,9 @@ void Init_gsl_matrix_complex(VALUE module)
   rb_define_method(cgsl_matrix_complex, "indgen", rb_gsl_matrix_complex_indgen, -1);
   rb_define_method(cgsl_matrix_complex, "indgen!", rb_gsl_matrix_complex_indgen_bang, -1);
   rb_define_singleton_method(cgsl_matrix_complex, "indgen", rb_gsl_matrix_complex_indgen_singleton, -1);
-  
+
   rb_define_method(cgsl_matrix_complex, "equal?", rb_gsl_matrix_complex_equal, -1);
   rb_define_alias(cgsl_matrix_complex, "==", "equal?");
   rb_define_method(cgsl_matrix_complex, "not_equal?", rb_gsl_matrix_complex_not_equal, -1);
-  rb_define_alias(cgsl_matrix_complex, "!=", "not_equal?");    
+  rb_define_alias(cgsl_matrix_complex, "!=", "not_equal?");
 }