rb-gsl 1.16.0.4 → 1.16.0.5

data/ext/{gsl → gsl_native}/array_complex.c

@@ -74,7 +74,7 @@ static VALUE rb_gsl_complex_arithmetics5(int flag, VALUE obj, VALUE bb)
     case GSL_COMPLEX_DIV:
       c = ALLOC(gsl_complex);
       *c = tmp;
-      return Data_Wrap_Struct(cgsl_complex, 0, free, c); 
+      return Data_Wrap_Struct(cgsl_complex, 0, free, c);
       break;
     }
     break;
@@ -87,83 +87,81 @@ static VALUE rb_gsl_complex_arithmetics5(int flag, VALUE obj, VALUE bb)
       case GSL_COMPLEX_SUB:
       case GSL_COMPLEX_MUL:
       case GSL_COMPLEX_DIV:
-	c = ALLOC(gsl_complex);
-	*c = tmp;
-	return Data_Wrap_Struct(cgsl_complex, 0, free, c); 
-	break;
+        c = ALLOC(gsl_complex);
+        *c = tmp;
+        return Data_Wrap_Struct(cgsl_complex, 0, free, c);
+        break;
       }
     } else {
       if (VECTOR_P(bb)) {
-	Data_Get_Struct(bb, gsl_vector, v);
-	cv = vector_to_complex(v);
-	cvnew = gsl_vector_complex_alloc(v->size);
-	if (cvnew == NULL) rb_raise(rb_eNoMemError, "gsl_vector_complex_alloc failed");
-	gsl_vector_complex_set_all(cvnew, *a);
-	switch (flag) {
-	case GSL_COMPLEX_ADD:
-	  gsl_vector_complex_add(cvnew, cv);
-	  break;
-	case GSL_COMPLEX_SUB:
-	  gsl_vector_complex_sub(cvnew, cv);
-	  break;
-	case GSL_COMPLEX_MUL:
-	  gsl_vector_complex_mul(cvnew, cv);
-	  break;
-	case GSL_COMPLEX_DIV:
-	  gsl_vector_complex_add(cvnew, cv);
-	  break;
-	}
-	gsl_vector_complex_free(cv);
-	return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, cvnew);
+        Data_Get_Struct(bb, gsl_vector, v);
+        cv = vector_to_complex(v);
+        cvnew = gsl_vector_complex_alloc(v->size);
+        if (cvnew == NULL) rb_raise(rb_eNoMemError, "gsl_vector_complex_alloc failed");
+        gsl_vector_complex_set_all(cvnew, *a);
+        switch (flag) {
+        case GSL_COMPLEX_ADD:
+          gsl_vector_complex_add(cvnew, cv);
+          break;
+        case GSL_COMPLEX_SUB:
+          gsl_vector_complex_sub(cvnew, cv);
+          break;
+        case GSL_COMPLEX_MUL:
+          gsl_vector_complex_mul(cvnew, cv);
+          break;
+        case GSL_COMPLEX_DIV:
+          gsl_vector_complex_add(cvnew, cv);
+          break;
+        }
+        gsl_vector_complex_free(cv);
+        return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, cvnew);
       }
       if (VECTOR_COMPLEX_P(bb)) {
-	Data_Get_Struct(bb, gsl_vector_complex, cv);
-	cvnew = gsl_vector_complex_alloc(v->size);
-	if (cvnew == NULL) rb_raise(rb_eNoMemError, "gsl_vector_complex_alloc failed");
-	gsl_vector_complex_set_all(cvnew, *a);
-	switch (flag) {
-	case GSL_COMPLEX_ADD:
-	  gsl_vector_complex_add(cvnew, cv);
-	  break;
-	case GSL_COMPLEX_SUB:
-	  gsl_vector_complex_sub(cvnew, cv);
-	  break;
-	case GSL_COMPLEX_MUL:
-	  gsl_vector_complex_mul(cvnew, cv);
-	  break;
-	case GSL_COMPLEX_DIV:
-	  gsl_vector_complex_add(cvnew, cv);
-	  break;
-	}
-	return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, cvnew);
+        Data_Get_Struct(bb, gsl_vector_complex, cv);
+        cvnew = gsl_vector_complex_alloc(v->size);
+        if (cvnew == NULL) rb_raise(rb_eNoMemError, "gsl_vector_complex_alloc failed");
+        gsl_vector_complex_set_all(cvnew, *a);
+        switch (flag) {
+        case GSL_COMPLEX_ADD:
+          gsl_vector_complex_add(cvnew, cv);
+          break;
+        case GSL_COMPLEX_SUB:
+          gsl_vector_complex_sub(cvnew, cv);
+          break;
+        case GSL_COMPLEX_MUL:
+          gsl_vector_complex_mul(cvnew, cv);
+          break;
+        case GSL_COMPLEX_DIV:
+          gsl_vector_complex_add(cvnew, cv);
+          break;
+        }
+        return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, cvnew);
       }
-
-
       if (MATRIX_P(bb)) {
-	Data_Get_Struct(bb, gsl_matrix, m);
-	cm = matrix_to_complex(m);
-	flagcm = 1;
-      }	else if (MATRIX_COMPLEX_P(bb)) {
-	Data_Get_Struct(bb, gsl_matrix_complex, cm);
+        Data_Get_Struct(bb, gsl_matrix, m);
+        cm = matrix_to_complex(m);
+        flagcm = 1;
+      } else if (MATRIX_COMPLEX_P(bb)) {
+        Data_Get_Struct(bb, gsl_matrix_complex, cm);
       } else {
-	rb_raise(rb_eTypeError, "wrong argument type %s", rb_class2name(CLASS_OF(bb)));
+        rb_raise(rb_eTypeError, "wrong argument type %s", rb_class2name(CLASS_OF(bb)));
       }
       cmself = gsl_matrix_complex_alloc(m->size1, m->size2);
       if (cmself == NULL) rb_raise(rb_eNoMemError, "gsl_matrix_complex_alloc failed");
       gsl_matrix_complex_set_all(cmself, *a);
       switch (flag) {
       case GSL_COMPLEX_ADD:
-	gsl_matrix_complex_add(cmself, cm);
-	break;
+        gsl_matrix_complex_add(cmself, cm);
+        break;
       case GSL_COMPLEX_SUB:
-	gsl_matrix_complex_sub(cmself, cm);
-	break;
+        gsl_matrix_complex_sub(cmself, cm);
+        break;
       case GSL_COMPLEX_MUL:
-	gsl_matrix_complex_mul_elements(cmself, cm);
-	break;
+        gsl_matrix_complex_mul_elements(cmself, cm);
+        break;
       case GSL_COMPLEX_DIV:
-	gsl_matrix_complex_div_elements(cmself, cm);
-	break;
+        gsl_matrix_complex_div_elements(cmself, cm);
+        break;
       }
       if (flagcm == 1) gsl_matrix_complex_free(cm);
       return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmself);
@@ -206,7 +204,7 @@ static VALUE rb_gsl_complex_coerce(VALUE obj, VALUE other)
     c = ALLOC(gsl_complex);
     *c = gsl_complex_rect(x, 0.0);
     return rb_ary_new3(2, Data_Wrap_Struct(cgsl_complex, 0, free, c),
-		       obj);
+                       obj);
     break;
   default:
     if (MATRIX_P(other)) {
@@ -219,7 +217,7 @@ static VALUE rb_gsl_complex_coerce(VALUE obj, VALUE other)
       gsl_matrix_complex_set_all(cmself, *c);
       vcmself = Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, cmself);
       return rb_ary_new3(2, vcmself, vcmnew);
-    } 
+    }
     if (MATRIX_COMPLEX_P(other)) {
       Data_Get_Struct(other, gsl_matrix_complex, cmnew);
       cmself = gsl_matrix_complex_alloc(cmnew->size1, cmnew->size2);

data/ext/{gsl → gsl_native}/blas1.c

@@ -16,10 +16,10 @@
 
 static int get_vector1(int argc, VALUE *argv, VALUE obj, gsl_vector **x);
 static int get_vector_complex1(int argc, VALUE *argv, VALUE obj, gsl_vector_complex **x);
-static int get_vector2(int argc, VALUE *argv, VALUE obj, 
-		       gsl_vector **x, gsl_vector **y);
-static int get_vector_complex2(int argc, VALUE *argv, VALUE obj, 
-		       gsl_vector_complex **x, gsl_vector_complex **y);
+static int get_vector2(int argc, VALUE *argv, VALUE obj,
+                       gsl_vector **x, gsl_vector **y);
+static int get_vector_complex2(int argc, VALUE *argv, VALUE obj,
+                               gsl_vector_complex **x, gsl_vector_complex **y);
 
 static int get_vector1(int argc, VALUE *argv, VALUE obj, gsl_vector **x)
 {
@@ -29,7 +29,7 @@ static int get_vector1(int argc, VALUE *argv, VALUE obj, gsl_vector **x)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     Data_Get_Vector(argv[0], (*x));
     break;
   default:
@@ -48,7 +48,7 @@ static int get_vector_complex1(int argc, VALUE *argv, VALUE obj, gsl_vector_comp
   case T_CLASS:
   case T_OBJECT:
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[0]);
     Data_Get_Struct(argv[0], gsl_vector_complex, (*x));
     break;
@@ -60,8 +60,8 @@ static int get_vector_complex1(int argc, VALUE *argv, VALUE obj, gsl_vector_comp
   return flag;
 }
 
-static int get_vector2(int argc, VALUE *argv, VALUE obj, 
-		       gsl_vector **x, gsl_vector **y)
+static int get_vector2(int argc, VALUE *argv, VALUE obj,
+                       gsl_vector **x, gsl_vector **y)
 {
   int flag = 0;
   switch (TYPE(obj)) {
@@ -69,13 +69,13 @@ static int get_vector2(int argc, VALUE *argv, VALUE obj,
   case T_CLASS:
   case T_OBJECT:
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     Data_Get_Vector(argv[0], (*x));
     Data_Get_Vector(argv[1], (*y));
     break;
   default:
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     Data_Get_Vector(obj, (*x));
     Data_Get_Vector(argv[0], (*y));
     flag = 1;
@@ -85,8 +85,8 @@ static int get_vector2(int argc, VALUE *argv, VALUE obj,
 }
 
 
-static int get_vector_complex2(int argc, VALUE *argv, VALUE obj, 
-		       gsl_vector_complex **x, gsl_vector_complex **y)
+static int get_vector_complex2(int argc, VALUE *argv, VALUE obj,
+                               gsl_vector_complex **x, gsl_vector_complex **y)
 {
   int flag = 0;
   switch (TYPE(obj)) {
@@ -94,7 +94,7 @@ static int get_vector_complex2(int argc, VALUE *argv, VALUE obj,
   case T_CLASS:
   case T_OBJECT:
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[0]);
     CHECK_VECTOR_COMPLEX(argv[1]);
     Data_Get_Struct(argv[0], gsl_vector_complex, (*x));
@@ -102,7 +102,7 @@ static int get_vector_complex2(int argc, VALUE *argv, VALUE obj,
     break;
   default:
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[0]);
     Data_Get_Struct(obj, gsl_vector_complex, (*x));
     Data_Get_Struct(argv[0], gsl_vector_complex, (*y));
@@ -242,7 +242,7 @@ static VALUE rb_gsl_blas_daxpy(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_vector, x);
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     //    a = RFLOAT(argv[0])->value;
     a = NUM2DBL(argv[0]);
@@ -269,7 +269,7 @@ static VALUE rb_gsl_blas_daxpy2(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_vector, x);
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     CHECK_VECTOR(argv[1]);
     //    a = RFLOAT(argv[0])->value;
@@ -298,7 +298,7 @@ static VALUE rb_gsl_blas_zaxpy(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_vector_complex, x);
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     CHECK_COMPLEX(argv[0]);
     CHECK_VECTOR_COMPLEX(argv[1]);
     Data_Get_Struct(argv[0], gsl_complex, a);
@@ -325,7 +325,7 @@ static VALUE rb_gsl_blas_zaxpy2(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_vector_complex, x);
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     CHECK_COMPLEX(argv[0]);
     CHECK_VECTOR_COMPLEX(argv[1]);
     Data_Get_Struct(argv[0], gsl_complex, a);
@@ -347,7 +347,7 @@ static VALUE rb_gsl_blas_dscal(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     CHECK_VECTOR(argv[1]);
     //    a = RFLOAT(argv[0])->value;
@@ -358,7 +358,7 @@ static VALUE rb_gsl_blas_dscal(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     //    a = RFLOAT(argv[0])->value;
     a = NUM2DBL(argv[0]);
@@ -379,7 +379,7 @@ static VALUE rb_gsl_blas_dscal2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     CHECK_VECTOR(argv[1]);
     a = NUM2DBL(argv[0]);
@@ -388,7 +388,7 @@ static VALUE rb_gsl_blas_dscal2(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_vector, x);
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     a = NUM2DBL(argv[0]);
     break;
@@ -408,7 +408,7 @@ static VALUE rb_gsl_blas_zdscal(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     CHECK_VECTOR_COMPLEX(argv[1]);
     //    a = RFLOAT(argv[0])->value;
@@ -420,7 +420,7 @@ static VALUE rb_gsl_blas_zdscal(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_vector_complex, x);
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     a = NUM2DBL(argv[0]);
     gsl_blas_zdscal(a, x);
@@ -438,7 +438,7 @@ static VALUE rb_gsl_blas_zdscal2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     CHECK_VECTOR_COMPLEX(argv[1]);
     a = NUM2DBL(argv[0]);
@@ -447,7 +447,7 @@ static VALUE rb_gsl_blas_zdscal2(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_vector_complex, x);
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     Need_Float(argv[0]);
     a = NUM2DBL(argv[0]);
     break;
@@ -468,7 +468,7 @@ static VALUE rb_gsl_blas_zscal(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[1]);
     Data_Get_Struct(argv[0], gsl_complex, a);
     Data_Get_Struct(argv[1], gsl_vector_complex, x);
@@ -477,7 +477,7 @@ static VALUE rb_gsl_blas_zscal(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     Data_Get_Struct(obj, gsl_vector_complex, x);
     Data_Get_Struct(argv[0], gsl_complex, a);
     gsl_blas_zscal(*a, x);
@@ -496,14 +496,14 @@ static VALUE rb_gsl_blas_zscal2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 2) rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[1]);
     Data_Get_Struct(argv[0], gsl_complex, a);
     Data_Get_Struct(argv[1], gsl_vector_complex, x);
     break;
   default:
     if (argc != 1) rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)",
-			    argc);
+                            argc);
     Data_Get_Struct(obj, gsl_vector_complex, x);
     Data_Get_Struct(argv[0], gsl_complex, a);
     break;
@@ -548,7 +548,7 @@ static VALUE rb_gsl_blas_drot2(VALUE obj, VALUE xx, VALUE yy, VALUE cc, VALUE ss
   gsl_vector_memcpy(ynew, y);
   gsl_blas_drot(xnew, ynew, c, s);
   return rb_ary_new3(2, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, xnew),
-		     Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, ynew));
+                     Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, ynew));
 }
 
 static VALUE rb_gsl_blas_drotm(VALUE obj, VALUE xx, VALUE yy, VALUE PP)
@@ -601,7 +601,7 @@ static VALUE rb_gsl_blas_drotm2(VALUE obj, VALUE xx, VALUE yy, VALUE PP)
   gsl_blas_drotm(xnew, ynew, p->data);
   if (flag == 1) gsl_vector_free(p);
   return rb_ary_new3(2, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, xnew),
-		     Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, ynew));
+                     Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, ynew));
 }
 
 void Init_gsl_blas1(VALUE module)

data/ext/{gsl → gsl_native}/blas2.c

@@ -57,7 +57,7 @@ static VALUE rb_gsl_blas_dgemv(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc < 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 4)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     Need_Float(argv[1]);
     CHECK_MATRIX(argv[2]);
@@ -71,7 +71,7 @@ static VALUE rb_gsl_blas_dgemv(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_matrix, A);
     if (argc < 3) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 3)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     Need_Float(argv[1]);
     CHECK_VECTOR(argv[2]);
@@ -100,7 +100,6 @@ static VALUE rb_gsl_blas_dgemv(int argc, VALUE *argv, VALUE obj)
   gsl_blas_dgemv(type, a, A, x, b, y);
   if (flag == 1) return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, y);
   else return argv[argc-1];
-  
 }
 
 static VALUE rb_gsl_blas_dgemv2(int argc, VALUE *argv, VALUE obj)
@@ -114,7 +113,7 @@ static VALUE rb_gsl_blas_dgemv2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc < 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 4)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     Need_Float(argv[1]);
     CHECK_MATRIX(argv[2]);
@@ -128,7 +127,7 @@ static VALUE rb_gsl_blas_dgemv2(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_matrix, A);
     if (argc < 3) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 3)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     Need_Float(argv[1]);
     CHECK_VECTOR(argv[2]);
@@ -172,7 +171,7 @@ static VALUE rb_gsl_blas_zgemv(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc < 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 4)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     CHECK_COMPLEX(argv[1]);
     CHECK_MATRIX_COMPLEX(argv[2]);
@@ -186,7 +185,7 @@ static VALUE rb_gsl_blas_zgemv(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_matrix_complex, A);
     if (argc < 3) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 3)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     CHECK_COMPLEX(argv[1]);
     CHECK_VECTOR_COMPLEX(argv[2]);
@@ -216,7 +215,6 @@ static VALUE rb_gsl_blas_zgemv(int argc, VALUE *argv, VALUE obj)
   gsl_blas_zgemv(type, *a, A, x, *b, y);
   if (flag == 1) return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, y);
   else return argv[argc-1];
-  
 }
 
 
@@ -231,7 +229,7 @@ static VALUE rb_gsl_blas_zgemv2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc < 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 4)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     CHECK_COMPLEX(argv[1]);
     CHECK_MATRIX_COMPLEX(argv[2]);
@@ -245,7 +243,7 @@ static VALUE rb_gsl_blas_zgemv2(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_matrix_complex, A);
     if (argc < 3) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 3)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     CHECK_COMPLEX(argv[1]);
     CHECK_VECTOR_COMPLEX(argv[2]);
@@ -288,7 +286,7 @@ static VALUE rb_gsl_blas_dtrmv(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 5) rb_raise(rb_eArgError, "wrong number of arguments (%d for 5)",
-			    argc);
+                            argc);
     CHECK_MATRIX(argv[3]);
     CHECK_VECTOR(argv[4]);
     Data_Get_Struct(argv[3], gsl_matrix, A);
@@ -296,7 +294,7 @@ static VALUE rb_gsl_blas_dtrmv(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for 4)",
-			    argc);
+                            argc);
     CHECK_VECTOR(argv[3]);
     Data_Get_Struct(obj, gsl_matrix, A);
     Data_Get_Struct(argv[3], gsl_vector, x);
@@ -304,7 +302,7 @@ static VALUE rb_gsl_blas_dtrmv(int argc, VALUE *argv, VALUE obj)
   }
   CHECK_FIXNUM(argv[0]); CHECK_FIXNUM(argv[1]); CHECK_FIXNUM(argv[2]);
   gsl_blas_dtrmv(FIX2INT(argv[0]), FIX2INT(argv[1]), FIX2INT(argv[2]),
-		 A, x);
+                 A, x);
   return argv[argc-1];
 }
 
@@ -318,7 +316,7 @@ static VALUE rb_gsl_blas_dtrmv2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 5) rb_raise(rb_eArgError, "wrong number of arguments (%d for 5)",
-			    argc);
+                            argc);
     CHECK_MATRIX(argv[3]);
     CHECK_VECTOR(argv[4]);
     Data_Get_Struct(argv[3], gsl_matrix, A);
@@ -326,7 +324,7 @@ static VALUE rb_gsl_blas_dtrmv2(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for 4)",
-			    argc);
+                            argc);
     CHECK_VECTOR(argv[3]);
     Data_Get_Struct(obj, gsl_matrix, A);
     Data_Get_Struct(argv[3], gsl_vector, x);
@@ -336,7 +334,7 @@ static VALUE rb_gsl_blas_dtrmv2(int argc, VALUE *argv, VALUE obj)
   xnew = gsl_vector_alloc(x->size);
   gsl_vector_memcpy(xnew, x);
   gsl_blas_dtrmv(FIX2INT(argv[0]), FIX2INT(argv[1]), FIX2INT(argv[2]),
-		 A, xnew);
+                 A, xnew);
   return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, xnew);
 }
 
@@ -350,7 +348,7 @@ static VALUE rb_gsl_blas_ztrmv(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 5) rb_raise(rb_eArgError, "wrong number of arguments (%d for 5)",
-			    argc);
+                            argc);
     CHECK_MATRIX_COMPLEX(argv[3]);
     CHECK_VECTOR_COMPLEX(argv[4]);
     Data_Get_Struct(argv[3], gsl_matrix_complex, A);
@@ -358,7 +356,7 @@ static VALUE rb_gsl_blas_ztrmv(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for 4)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[3]);
     Data_Get_Struct(obj, gsl_matrix_complex, A);
     Data_Get_Struct(argv[3], gsl_vector_complex, x);
@@ -366,7 +364,7 @@ static VALUE rb_gsl_blas_ztrmv(int argc, VALUE *argv, VALUE obj)
   }
   CHECK_FIXNUM(argv[0]); CHECK_FIXNUM(argv[1]); CHECK_FIXNUM(argv[2]);
   gsl_blas_ztrmv(FIX2INT(argv[0]), FIX2INT(argv[1]), FIX2INT(argv[2]),
-		 A, x);
+                 A, x);
   return argv[argc-1];
 }
 
@@ -379,7 +377,7 @@ static VALUE rb_gsl_blas_ztrmv2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 5) rb_raise(rb_eArgError, "wrong number of arguments (%d for 5)",
-			    argc);
+                            argc);
     CHECK_MATRIX_COMPLEX(argv[3]);
     CHECK_VECTOR_COMPLEX(argv[4]);
     Data_Get_Struct(argv[3], gsl_matrix_complex, A);
@@ -387,7 +385,7 @@ static VALUE rb_gsl_blas_ztrmv2(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for 4)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[3]);
     Data_Get_Struct(obj, gsl_matrix_complex, A);
     Data_Get_Struct(argv[3], gsl_vector_complex, x);
@@ -397,7 +395,7 @@ static VALUE rb_gsl_blas_ztrmv2(int argc, VALUE *argv, VALUE obj)
   xnew = gsl_vector_complex_alloc(x->size);
   gsl_vector_complex_memcpy(xnew, x);
   gsl_blas_ztrmv(FIX2INT(argv[0]), FIX2INT(argv[1]), FIX2INT(argv[2]),
-		 A, xnew);
+                 A, xnew);
   return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, xnew);
 }
 
@@ -410,7 +408,7 @@ static VALUE rb_gsl_blas_dtrsv(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 5) rb_raise(rb_eArgError, "wrong number of arguments (%d for 5)",
-			    argc);
+                            argc);
     CHECK_MATRIX(argv[3]);
     CHECK_VECTOR(argv[4]);
     Data_Get_Struct(argv[3], gsl_matrix, A);
@@ -418,7 +416,7 @@ static VALUE rb_gsl_blas_dtrsv(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for 4)",
-			    argc);
+                            argc);
     CHECK_VECTOR(argv[3]);
     Data_Get_Struct(obj, gsl_matrix, A);
     Data_Get_Struct(argv[3], gsl_vector, x);
@@ -426,7 +424,7 @@ static VALUE rb_gsl_blas_dtrsv(int argc, VALUE *argv, VALUE obj)
   }
   CHECK_FIXNUM(argv[0]); CHECK_FIXNUM(argv[1]); CHECK_FIXNUM(argv[2]);
   gsl_blas_dtrsv(FIX2INT(argv[0]), FIX2INT(argv[1]), FIX2INT(argv[2]),
-		 A, x);
+                 A, x);
   return argv[argc-1];
 }
 
@@ -439,7 +437,7 @@ static VALUE rb_gsl_blas_dtrsv2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 5) rb_raise(rb_eArgError, "wrong number of arguments (%d for 5)",
-			    argc);
+                            argc);
     CHECK_MATRIX(argv[3]);
     CHECK_VECTOR(argv[4]);
     Data_Get_Struct(argv[3], gsl_matrix, A);
@@ -447,7 +445,7 @@ static VALUE rb_gsl_blas_dtrsv2(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for 4)",
-			    argc);
+                            argc);
     CHECK_VECTOR(argv[3]);
     Data_Get_Struct(obj, gsl_matrix, A);
     Data_Get_Struct(argv[3], gsl_vector, x);
@@ -457,7 +455,7 @@ static VALUE rb_gsl_blas_dtrsv2(int argc, VALUE *argv, VALUE obj)
   xnew = gsl_vector_alloc(x->size);
   gsl_vector_memcpy(xnew, x);
   gsl_blas_dtrsv(FIX2INT(argv[0]), FIX2INT(argv[1]), FIX2INT(argv[2]),
-		 A, xnew);
+                 A, xnew);
   return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, xnew);
 }
 
@@ -470,7 +468,7 @@ static VALUE rb_gsl_blas_ztrsv(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 5) rb_raise(rb_eArgError, "wrong number of arguments (%d for 5)",
-			    argc);
+                            argc);
     CHECK_MATRIX_COMPLEX(argv[3]);
     CHECK_VECTOR_COMPLEX(argv[4]);
     Data_Get_Struct(argv[3], gsl_matrix_complex, A);
@@ -478,7 +476,7 @@ static VALUE rb_gsl_blas_ztrsv(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for 4)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[3]);
     Data_Get_Struct(obj, gsl_matrix_complex, A);
     Data_Get_Struct(argv[3], gsl_vector_complex, x);
@@ -486,7 +484,7 @@ static VALUE rb_gsl_blas_ztrsv(int argc, VALUE *argv, VALUE obj)
   }
   CHECK_FIXNUM(argv[0]); CHECK_FIXNUM(argv[1]); CHECK_FIXNUM(argv[2]);
   gsl_blas_ztrsv(FIX2INT(argv[0]), FIX2INT(argv[1]), FIX2INT(argv[2]),
-		 A, x);
+                 A, x);
   return argv[argc-1];
 }
 
@@ -499,7 +497,7 @@ static VALUE rb_gsl_blas_ztrsv2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc != 5) rb_raise(rb_eArgError, "wrong number of arguments (%d for 5)",
-			    argc);
+                            argc);
     CHECK_MATRIX_COMPLEX(argv[3]);
     CHECK_VECTOR_COMPLEX(argv[4]);
     Data_Get_Struct(argv[3], gsl_matrix_complex, A);
@@ -507,7 +505,7 @@ static VALUE rb_gsl_blas_ztrsv2(int argc, VALUE *argv, VALUE obj)
     break;
   default:
     if (argc != 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for 4)",
-			    argc);
+                            argc);
     CHECK_VECTOR_COMPLEX(argv[3]);
     Data_Get_Struct(obj, gsl_matrix_complex, A);
     Data_Get_Struct(argv[3], gsl_vector_complex, x);
@@ -517,7 +515,7 @@ static VALUE rb_gsl_blas_ztrsv2(int argc, VALUE *argv, VALUE obj)
   xnew = gsl_vector_complex_alloc(x->size);
   gsl_vector_complex_memcpy(xnew, x);
   gsl_blas_ztrsv(FIX2INT(argv[0]), FIX2INT(argv[1]), FIX2INT(argv[2]),
-		 A, xnew);
+                 A, xnew);
   return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, xnew);
 }
 
@@ -532,7 +530,7 @@ static VALUE rb_gsl_blas_dsymv(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc < 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 4)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     Need_Float(argv[1]);
     CHECK_MATRIX(argv[2]);
@@ -546,7 +544,7 @@ static VALUE rb_gsl_blas_dsymv(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_matrix, A);
     if (argc < 3) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 3)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     Need_Float(argv[1]);
     CHECK_VECTOR(argv[2]);
@@ -575,7 +573,6 @@ static VALUE rb_gsl_blas_dsymv(int argc, VALUE *argv, VALUE obj)
   gsl_blas_dsymv(type, a, A, x, b, y);
   if (flag == 1) return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, y);
   else return argv[argc-1];
-  
 }
 
 static VALUE rb_gsl_blas_dsymv2(int argc, VALUE *argv, VALUE obj)
@@ -589,7 +586,7 @@ static VALUE rb_gsl_blas_dsymv2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc < 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 4)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     Need_Float(argv[1]);
     CHECK_MATRIX(argv[2]);
@@ -603,7 +600,7 @@ static VALUE rb_gsl_blas_dsymv2(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_matrix, A);
     if (argc < 3) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 3)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     Need_Float(argv[1]);
     CHECK_VECTOR(argv[2]);
@@ -647,7 +644,7 @@ static VALUE rb_gsl_blas_zhemv(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc < 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 4)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     CHECK_COMPLEX(argv[1]);
     CHECK_MATRIX_COMPLEX(argv[2]);
@@ -661,7 +658,7 @@ static VALUE rb_gsl_blas_zhemv(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_matrix_complex, A);
     if (argc < 3) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 3)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     CHECK_COMPLEX(argv[1]);
     CHECK_VECTOR_COMPLEX(argv[2]);
@@ -691,7 +688,6 @@ static VALUE rb_gsl_blas_zhemv(int argc, VALUE *argv, VALUE obj)
   gsl_blas_zhemv(type, *a, A, x, *b, y);
   if (flag == 1) return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, y);
   else return argv[argc-1];
-  
 }
 
 static VALUE rb_gsl_blas_zhemv2(int argc, VALUE *argv, VALUE obj)
@@ -705,7 +701,7 @@ static VALUE rb_gsl_blas_zhemv2(int argc, VALUE *argv, VALUE obj)
   case T_CLASS:
   case T_OBJECT:
     if (argc < 4) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 4)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     CHECK_COMPLEX(argv[1]);
     CHECK_MATRIX_COMPLEX(argv[2]);
@@ -719,7 +715,7 @@ static VALUE rb_gsl_blas_zhemv2(int argc, VALUE *argv, VALUE obj)
   default:
     Data_Get_Struct(obj, gsl_matrix_complex, A);
     if (argc < 3) rb_raise(rb_eArgError, "wrong number of arguments (%d for >= 3)",
-			   argc);
+                           argc);
     CHECK_FIXNUM(argv[0]);
     CHECK_COMPLEX(argv[1]);
     CHECK_VECTOR_COMPLEX(argv[2]);
@@ -759,7 +755,7 @@ static VALUE rb_gsl_blas_dger(VALUE obj, VALUE aa, VALUE xx, VALUE yy, VALUE AA)
   gsl_vector *x = NULL, *y = NULL;
   double a;
   Need_Float(aa);
-  CHECK_VECTOR(xx);  CHECK_VECTOR(yy);  
+  CHECK_VECTOR(xx);  CHECK_VECTOR(yy);
   CHECK_MATRIX(AA);
   a = NUM2DBL(aa);
   Data_Get_Struct(xx, gsl_vector, x);
@@ -775,7 +771,7 @@ static VALUE rb_gsl_blas_dger2(VALUE obj, VALUE aa, VALUE xx, VALUE yy, VALUE AA
   gsl_vector *x = NULL, *y = NULL;
   double a;
   Need_Float(aa);
-  CHECK_VECTOR(xx);  CHECK_VECTOR(yy);  
+  CHECK_VECTOR(xx);  CHECK_VECTOR(yy);
   CHECK_MATRIX(AA);
   a = NUM2DBL(aa);
   Data_Get_Struct(xx, gsl_vector, x);
@@ -793,7 +789,7 @@ static VALUE rb_gsl_blas_zgeru(VALUE obj, VALUE aa, VALUE xx, VALUE yy, VALUE AA
   gsl_vector_complex *x = NULL, *y = NULL;
   gsl_complex *a;
   CHECK_COMPLEX(aa);
-  CHECK_VECTOR_COMPLEX(xx);CHECK_VECTOR_COMPLEX(yy); 
+  CHECK_VECTOR_COMPLEX(xx); CHECK_VECTOR_COMPLEX(yy);
   CHECK_MATRIX_COMPLEX(AA);
   Data_Get_Struct(aa, gsl_complex, a);
   Data_Get_Struct(xx, gsl_vector_complex, x);
@@ -809,7 +805,7 @@ static VALUE rb_gsl_blas_zgeru2(VALUE obj, VALUE aa, VALUE xx, VALUE yy, VALUE A
   gsl_vector_complex *x = NULL, *y = NULL;
   gsl_complex *a;
   CHECK_COMPLEX(aa);
-  CHECK_VECTOR_COMPLEX(xx);CHECK_VECTOR_COMPLEX(yy); 
+  CHECK_VECTOR_COMPLEX(xx); CHECK_VECTOR_COMPLEX(yy);
   CHECK_MATRIX_COMPLEX(AA);
   Data_Get_Struct(aa, gsl_complex, a);
   Data_Get_Struct(xx, gsl_vector_complex, x);
@@ -827,7 +823,7 @@ static VALUE rb_gsl_blas_zgerc(VALUE obj, VALUE aa, VALUE xx, VALUE yy, VALUE AA
   gsl_vector_complex *x = NULL, *y = NULL;
   gsl_complex *a;
   CHECK_COMPLEX(aa);
-  CHECK_VECTOR_COMPLEX(xx);CHECK_VECTOR_COMPLEX(yy); 
+  CHECK_VECTOR_COMPLEX(xx); CHECK_VECTOR_COMPLEX(yy);
   CHECK_MATRIX_COMPLEX(AA);
   Data_Get_Struct(aa, gsl_complex, a);
   Data_Get_Struct(xx, gsl_vector_complex, x);
@@ -843,7 +839,7 @@ static VALUE rb_gsl_blas_zgerc2(VALUE obj, VALUE aa, VALUE xx, VALUE yy, VALUE A
   gsl_vector_complex *x = NULL, *y = NULL;
   gsl_complex *a;
   CHECK_COMPLEX(aa);
-  CHECK_VECTOR_COMPLEX(xx);CHECK_VECTOR_COMPLEX(yy); 
+  CHECK_VECTOR_COMPLEX(xx); CHECK_VECTOR_COMPLEX(yy);
   CHECK_MATRIX_COMPLEX(AA);
   Data_Get_Struct(aa, gsl_complex, a);
   Data_Get_Struct(xx, gsl_vector_complex, x);
@@ -919,8 +915,8 @@ static VALUE rb_gsl_blas_zher_a(VALUE obj, VALUE tt, VALUE aa, VALUE xx, VALUE A
   return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, Anew);
 }
 
-static VALUE rb_gsl_blas_dsyr2(VALUE obj, VALUE tt, VALUE aa, VALUE xx, 
-			      VALUE yy, VALUE AA)
+static VALUE rb_gsl_blas_dsyr2(VALUE obj, VALUE tt, VALUE aa, VALUE xx,
+                               VALUE yy, VALUE AA)
 {
   gsl_matrix *A = NULL;
   gsl_vector *x = NULL, *y = NULL;
@@ -936,8 +932,8 @@ static VALUE rb_gsl_blas_dsyr2(VALUE obj, VALUE tt, VALUE aa, VALUE xx,
   return AA;
 }
 
-static VALUE rb_gsl_blas_dsyr2_a(VALUE obj, VALUE tt, VALUE aa, VALUE xx, 
-				 VALUE yy, VALUE AA)
+static VALUE rb_gsl_blas_dsyr2_a(VALUE obj, VALUE tt, VALUE aa, VALUE xx,
+                                 VALUE yy, VALUE AA)
 {
   gsl_matrix *A = NULL, *Anew = NULL;
   gsl_vector *x = NULL, *y = NULL;
@@ -955,8 +951,8 @@ static VALUE rb_gsl_blas_dsyr2_a(VALUE obj, VALUE tt, VALUE aa, VALUE xx,
   return Data_Wrap_Struct(cgsl_matrix, 0, gsl_matrix_free, Anew);
 }
 
-static VALUE rb_gsl_blas_zher2(VALUE obj, VALUE tt, VALUE aa, VALUE xx, 
-			      VALUE yy, VALUE AA)
+static VALUE rb_gsl_blas_zher2(VALUE obj, VALUE tt, VALUE aa, VALUE xx,
+                               VALUE yy, VALUE AA)
 {
   gsl_matrix_complex *A = NULL;
   gsl_vector_complex *x = NULL, *y = NULL;
@@ -972,8 +968,8 @@ static VALUE rb_gsl_blas_zher2(VALUE obj, VALUE tt, VALUE aa, VALUE xx,
   return AA;
 }
 
-static VALUE rb_gsl_blas_zher2_a(VALUE obj, VALUE tt, VALUE aa, VALUE xx, 
-				 VALUE yy, VALUE AA)
+static VALUE rb_gsl_blas_zher2_a(VALUE obj, VALUE tt, VALUE aa, VALUE xx,
+                                 VALUE yy, VALUE AA)
 {
   gsl_matrix_complex *A = NULL, *Anew = NULL;
   gsl_vector_complex *x = NULL, *y = NULL;