gsl 1.15.3 → 1.16.0.6

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

@@ -9,9 +9,9 @@
   WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
-#include "rb_gsl_config.h"
-#include "rb_gsl_complex.h"
-#include "rb_gsl_array.h"
+
+#include "include/rb_gsl_complex.h"
+#include "include/rb_gsl_array.h"
 
 VALUE cgsl_complex;
 
@@ -31,11 +31,9 @@ gsl_complex rb_gsl_obj_to_gsl_complex(VALUE obj, gsl_complex *z)
     z = &tmp;
     GSL_SET_COMPLEX(z, 0.0, 0.0);
   }
-
   if(obj == Qnil) {
     return *z;
   }
-
   switch(TYPE(obj)) {
   case T_ARRAY:
     vre = rb_ary_entry(obj,0);
@@ -54,8 +52,8 @@ gsl_complex rb_gsl_obj_to_gsl_complex(VALUE obj, gsl_complex *z)
       *z = *zz;
     } else {
       rb_raise(rb_eTypeError,
-          "wrong type %s, (nil, Array, Float, Integer, or GSL::Complex expected)",
-          rb_class2name(CLASS_OF(obj)));
+               "wrong type %s, (nil, Array, Float, Integer, or GSL::Complex expected)",
+               rb_class2name(CLASS_OF(obj)));
     }
     break;
   }
@@ -66,7 +64,7 @@ static VALUE rb_gsl_complex_new(int argc, VALUE *argv, VALUE klass)
 {
   gsl_complex *c = NULL;
   VALUE obj;
-  obj = Data_Make_Struct(klass, gsl_complex, 0, free, c); 
+  obj = Data_Make_Struct(klass, gsl_complex, 0, free, c);
   switch (argc) {
   case 1:
     switch (TYPE(argv[0])) {
@@ -80,8 +78,8 @@ static VALUE rb_gsl_complex_new(int argc, VALUE *argv, VALUE klass)
       *c = gsl_complex_rect(NUM2DBL(argv[0]), 0.0);
       break;
     default:
-      rb_raise(rb_eTypeError, "wrong argument type %s", 
-	       rb_class2name(CLASS_OF(argv[0])));
+      rb_raise(rb_eTypeError, "wrong argument type %s",
+               rb_class2name(CLASS_OF(argv[0])));
     }
     break;
   case 2:
@@ -107,7 +105,7 @@ static VALUE rb_gsl_complex_polar(VALUE klass, VALUE r, VALUE theta)
   VALUE obj;
   gsl_complex *c = NULL;
   Need_Float(r); Need_Float(theta);
-  obj = Data_Make_Struct(klass, gsl_complex, 0, free, c); 
+  obj = Data_Make_Struct(klass, gsl_complex, 0, free, c);
   *c = gsl_complex_polar(NUM2DBL(r), NUM2DBL(theta));
   return obj;
 }
@@ -121,10 +119,10 @@ static VALUE rb_gsl_complex_get(VALUE obj, VALUE ii)
   i = FIX2INT(ii);
   switch (i) {
   case 0:
-    return  rb_float_new(GSL_REAL(*c));
+    return rb_float_new(GSL_REAL(*c));
     break;
   case 1:
-    return  rb_float_new(GSL_IMAG(*c));
+    return rb_float_new(GSL_IMAG(*c));
     break;
   default:
     rb_raise(rb_eArgError, "wrong argument (%d for 0 or 1)", i);
@@ -137,14 +135,14 @@ static VALUE rb_gsl_complex_real(VALUE obj)
 {
   gsl_complex *c = NULL;
   Data_Get_Struct(obj, gsl_complex, c);
-  return  rb_float_new(GSL_REAL(*c));
+  return rb_float_new(GSL_REAL(*c));
 }
 
 static VALUE rb_gsl_complex_imag(VALUE obj)
 {
   gsl_complex *c = NULL;
   Data_Get_Struct(obj, gsl_complex, c);
-  return  rb_float_new(GSL_IMAG(*c));
+  return rb_float_new(GSL_IMAG(*c));
 }
 
 static VALUE rb_gsl_complex_print(VALUE obj)
@@ -172,7 +170,7 @@ static VALUE rb_gsl_complex_return_double(double (*func)(gsl_complex), VALUE obj
 {
   gsl_complex *c = NULL;
   Data_Get_Struct(obj, gsl_complex, c);
-  return  rb_float_new((*func)(*c));
+  return rb_float_new((*func)(*c));
 }
 
 static VALUE rb_gsl_complex_arg(VALUE obj)
@@ -196,10 +194,10 @@ static VALUE rb_gsl_complex_logabs(VALUE obj)
 }
 
 static VALUE rb_gsl_complex_arithmetics2(gsl_complex (*func)(gsl_complex, double),
-					 VALUE obj, VALUE xx);
+                                         VALUE obj, VALUE xx);
 
 static VALUE rb_gsl_complex_arithmetics2(gsl_complex (*func)(gsl_complex, double),
-					VALUE obj, VALUE xx)
+                                         VALUE obj, VALUE xx)
 {
   gsl_complex *a = NULL, *c = NULL, tmp;
   VALUE obj2;
@@ -208,7 +206,7 @@ static VALUE rb_gsl_complex_arithmetics2(gsl_complex (*func)(gsl_complex, double
   Data_Get_Struct(obj, gsl_complex, a);
   x = NUM2DBL(xx);
   tmp = (*func)(*a, x);
-  obj2 = Data_Make_Struct(cgsl_complex, gsl_complex, 0, free, c); 
+  obj2 = Data_Make_Struct(cgsl_complex, gsl_complex, 0, free, c);
   *c = tmp;
   return obj2;
 }
@@ -262,7 +260,7 @@ static VALUE rb_gsl_complex_operate(gsl_complex (*func)(gsl_complex), VALUE obj)
   Data_Get_Struct(obj, gsl_complex, c);
   cnew = ALLOC(gsl_complex);
   *cnew = (*func)(*c);
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_operate2(gsl_complex (*func)(gsl_complex), int argc, VALUE *argv, VALUE obj)
@@ -279,33 +277,33 @@ static VALUE rb_gsl_complex_operate2(gsl_complex (*func)(gsl_complex), int argc,
     case 1:
       switch (TYPE(argv[0])) {
       case T_ARRAY:
-	tmp = ary2complex(argv[0]);
-	c = &tmp;
-	break;
+        tmp = ary2complex(argv[0]);
+        c = &tmp;
+        break;
       default:
-	if (VECTOR_COMPLEX_P(argv[0])) {
-	  Data_Get_Struct(argv[0], gsl_vector_complex, v);
-	  vnew = gsl_vector_complex_alloc(v->size);
-	  for (i = 0; i < v->size; i++) {
-	    c = GSL_COMPLEX_AT(v, i);
-	    gsl_vector_complex_set(vnew, i, (*func)(*c));
-	  }
-	  return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_free, vnew);
-	} else if (MATRIX_COMPLEX_P(obj)) {
-	  Data_Get_Struct(obj, gsl_matrix_complex, m);
-	  mnew = gsl_matrix_complex_alloc(m->size1, m->size2);
-	  for (i = 0; i < m->size1; i++) {
-	    for (j = 0; j < m->size2; j++) {
-	      gsl_matrix_complex_set(mnew, i, j, 
-				     (*func)(gsl_matrix_complex_get(m, i, j)));
-	    }
-	  }
-	  return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, mnew);
-	} else {
-	  CHECK_COMPLEX(argv[0]);
-	  Data_Get_Struct(argv[0], gsl_complex, c);
-	}
-	break;
+        if (VECTOR_COMPLEX_P(argv[0])) {
+          Data_Get_Struct(argv[0], gsl_vector_complex, v);
+          vnew = gsl_vector_complex_alloc(v->size);
+          for (i = 0; i < v->size; i++) {
+            c = GSL_COMPLEX_AT(v, i);
+            gsl_vector_complex_set(vnew, i, (*func)(*c));
+          }
+          return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_free, vnew);
+        } else if (MATRIX_COMPLEX_P(obj)) {
+          Data_Get_Struct(obj, gsl_matrix_complex, m);
+          mnew = gsl_matrix_complex_alloc(m->size1, m->size2);
+          for (i = 0; i < m->size1; i++) {
+            for (j = 0; j < m->size2; j++) {
+              gsl_matrix_complex_set(mnew, i, j,
+                                     (*func)(gsl_matrix_complex_get(m, i, j)));
+            }
+          }
+          return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, mnew);
+        } else {
+          CHECK_COMPLEX(argv[0]);
+          Data_Get_Struct(argv[0], gsl_complex, c);
+        }
+        break;
       }
       break;
     case 2:
@@ -324,7 +322,7 @@ static VALUE rb_gsl_complex_operate2(gsl_complex (*func)(gsl_complex), int argc,
   }
   cnew = ALLOC(gsl_complex);
   *cnew = (*func)(*c);
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_conjugate(VALUE obj)
@@ -350,7 +348,7 @@ static VALUE rb_gsl_complex_sqrt_real(VALUE obj, VALUE x)
   tmp = gsl_complex_sqrt_real(NUM2DBL(x));
   cnew = ALLOC(gsl_complex);
   *cnew = tmp;
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_sqrt(int argc, VALUE *argv, VALUE obj)
@@ -362,22 +360,22 @@ static VALUE rb_gsl_complex_sqrt(int argc, VALUE *argv, VALUE obj)
     case 1:
       switch (TYPE(argv[0])) {
       case T_FIXNUM: case T_FLOAT:
-	return rb_gsl_complex_sqrt_real(obj, argv[0]);
-	break;
+        return rb_gsl_complex_sqrt_real(obj, argv[0]);
+        break;
       case T_ARRAY:
-	tmp = ary2complex(argv[0]);
-	z = ALLOC(gsl_complex);
-	*z = gsl_complex_sqrt(tmp);
-	return Data_Wrap_Struct(cgsl_complex, 0, free, z);
-	break;
+        tmp = ary2complex(argv[0]);
+        z = ALLOC(gsl_complex);
+        *z = gsl_complex_sqrt(tmp);
+        return Data_Wrap_Struct(cgsl_complex, 0, free, z);
+        break;
       default:
-	CHECK_COMPLEX(argv[0]);
-	Data_Get_Struct(argv[0], gsl_complex, z);
-	tmp = *z;
-	z = ALLOC(gsl_complex);
-	*z = gsl_complex_sqrt(tmp);
-	return Data_Wrap_Struct(cgsl_complex, 0, free, z);      
-	break;
+        CHECK_COMPLEX(argv[0]);
+        Data_Get_Struct(argv[0], gsl_complex, z);
+        tmp = *z;
+        z = ALLOC(gsl_complex);
+        *z = gsl_complex_sqrt(tmp);
+        return Data_Wrap_Struct(cgsl_complex, 0, free, z);
+        break;
       }
       break;
     case 2:
@@ -385,7 +383,7 @@ static VALUE rb_gsl_complex_sqrt(int argc, VALUE *argv, VALUE obj)
       GSL_SET_REAL(&tmp, NUM2DBL(argv[0]));
       GSL_SET_IMAG(&tmp, NUM2DBL(argv[1]));
       *z = gsl_complex_sqrt(tmp);
-      return Data_Wrap_Struct(cgsl_complex, 0, free, z);     
+      return Data_Wrap_Struct(cgsl_complex, 0, free, z);
       break;
     default:
       rb_raise(rb_eArgError, "wrong number of arguments (%d for 1 or 2)", argc);
@@ -408,14 +406,14 @@ VALUE rb_gsl_complex_pow(int argc, VALUE *argv, VALUE obj)
   case T_MODULE:
   case T_CLASS:
   case T_OBJECT:
-    if (argc != 2) 
+    if (argc != 2)
       rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)", argc);
     switch (TYPE(argv[1])) {
     case T_ARRAY:
       tmpa = ary2complex(argv[1]);
       a = &tmpa;
       break;
-    case T_FLOAT: case T_FIXNUM: 
+    case T_FLOAT: case T_FIXNUM:
       return rb_gsl_complex_pow_real(argc, argv, obj);
       break;
     default:
@@ -431,24 +429,24 @@ VALUE rb_gsl_complex_pow(int argc, VALUE *argv, VALUE obj)
       break;
     default:
       if (VECTOR_COMPLEX_P(argv[0])) {
-	Data_Get_Struct(argv[0], gsl_vector_complex, v);
-	vnew = gsl_vector_complex_alloc(v->size);
-	for (i = 0; i < v->size; i++) {
-	  c = GSL_COMPLEX_AT(v, i);
-	  gsl_vector_complex_set(vnew, i, gsl_complex_pow(*c, *a));
-	}
-	return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vnew);
+        Data_Get_Struct(argv[0], gsl_vector_complex, v);
+        vnew = gsl_vector_complex_alloc(v->size);
+        for (i = 0; i < v->size; i++) {
+          c = GSL_COMPLEX_AT(v, i);
+          gsl_vector_complex_set(vnew, i, gsl_complex_pow(*c, *a));
+        }
+        return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vnew);
       }
       if (MATRIX_COMPLEX_P(argv[0])) {
-	Data_Get_Struct(argv[0], gsl_matrix_complex, m);
-	mnew = gsl_matrix_complex_alloc(m->size1, m->size2);
-	for (i = 0; i < m->size1; i++) {
-	  for (j = 0; j < m->size2; j++) {
-	    c = gsl_matrix_complex_ptr(m, i, j);
-	    gsl_matrix_complex_set(mnew, i, j, gsl_complex_pow(*c, *a));
-	  }
-	}
-	return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, mnew);
+        Data_Get_Struct(argv[0], gsl_matrix_complex, m);
+        mnew = gsl_matrix_complex_alloc(m->size1, m->size2);
+        for (i = 0; i < m->size1; i++) {
+          for (j = 0; j < m->size2; j++) {
+            c = gsl_matrix_complex_ptr(m, i, j);
+            gsl_matrix_complex_set(mnew, i, j, gsl_complex_pow(*c, *a));
+          }
+        }
+        return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, mnew);
       }
       CHECK_COMPLEX(argv[0]);
       Data_Get_Struct(argv[0], gsl_complex, c);
@@ -457,7 +455,7 @@ VALUE rb_gsl_complex_pow(int argc, VALUE *argv, VALUE obj)
 
     break;
   default:
-    if (argc != 1) 
+    if (argc != 1)
       rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)", argc);
     CHECK_COMPLEX(obj);
     Data_Get_Struct(obj, gsl_complex, c);
@@ -466,7 +464,7 @@ VALUE rb_gsl_complex_pow(int argc, VALUE *argv, VALUE obj)
       tmpa = ary2complex(argv[0]);
       a = &tmpa;
       break;
-    case T_FLOAT: case T_FIXNUM: 
+    case T_FLOAT: case T_FIXNUM:
       return rb_gsl_complex_pow_real(argc, argv, obj);
       break;
     default:
@@ -478,7 +476,7 @@ VALUE rb_gsl_complex_pow(int argc, VALUE *argv, VALUE obj)
   }
   cnew = ALLOC(gsl_complex);
   *cnew = gsl_complex_pow(*c, *a);
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 VALUE rb_gsl_complex_pow_real(int argc, VALUE *argv, VALUE obj)
@@ -492,7 +490,7 @@ VALUE rb_gsl_complex_pow_real(int argc, VALUE *argv, VALUE obj)
   case T_MODULE:
   case T_CLASS:
   case T_OBJECT:
-    if (argc != 2) 
+    if (argc != 2)
       rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)", argc);
     switch (TYPE(argv[0])) {
     case T_ARRAY:
@@ -502,26 +500,26 @@ VALUE rb_gsl_complex_pow_real(int argc, VALUE *argv, VALUE obj)
       break;
     default:
       if (VECTOR_COMPLEX_P(argv[0])) {
-	Data_Get_Struct(argv[0], gsl_vector_complex, v);
-	vnew = gsl_vector_complex_alloc(v->size);
-	a = NUM2DBL(argv[1]);
-	for (i = 0; i < v->size; i++) {
-	  c = GSL_COMPLEX_AT(v, i);
-	  tmpc = gsl_complex_pow_real(*c, a);
-	  gsl_vector_complex_set(vnew, i, tmpc);
-	}
-	return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_free, vnew);
+        Data_Get_Struct(argv[0], gsl_vector_complex, v);
+        vnew = gsl_vector_complex_alloc(v->size);
+        a = NUM2DBL(argv[1]);
+        for (i = 0; i < v->size; i++) {
+          c = GSL_COMPLEX_AT(v, i);
+          tmpc = gsl_complex_pow_real(*c, a);
+          gsl_vector_complex_set(vnew, i, tmpc);
+        }
+        return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_free, vnew);
       }
       if (MATRIX_COMPLEX_P(argv[0])) {
-	Data_Get_Struct(argv[0], gsl_matrix_complex, m);
-	mnew = gsl_matrix_complex_alloc(m->size1, m->size2);
-	for (i = 0; i < m->size1; i++) {
-	  for (j = 0; j < m->size2; j++) {
-	    tmpc = gsl_complex_pow_real(gsl_matrix_complex_get(m, i, j), a);
-	    gsl_matrix_complex_set(mnew, i, j, tmpc);
-	  }
-	}
-	return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, mnew);
+        Data_Get_Struct(argv[0], gsl_matrix_complex, m);
+        mnew = gsl_matrix_complex_alloc(m->size1, m->size2);
+        for (i = 0; i < m->size1; i++) {
+          for (j = 0; j < m->size2; j++) {
+            tmpc = gsl_complex_pow_real(gsl_matrix_complex_get(m, i, j), a);
+            gsl_matrix_complex_set(mnew, i, j, tmpc);
+          }
+        }
+        return Data_Wrap_Struct(cgsl_matrix_complex, 0, gsl_matrix_complex_free, mnew);
       }
       CHECK_COMPLEX(argv[0]);
       Data_Get_Struct(argv[0], gsl_complex, c);
@@ -531,7 +529,7 @@ VALUE rb_gsl_complex_pow_real(int argc, VALUE *argv, VALUE obj)
     a = NUM2DBL(argv[1]);
     break;
   default:
-    if (argc != 1) 
+    if (argc != 1)
       rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)", argc);
     CHECK_COMPLEX(obj);
     Need_Float(argv[0]);
@@ -541,7 +539,7 @@ VALUE rb_gsl_complex_pow_real(int argc, VALUE *argv, VALUE obj)
   }
   cnew = ALLOC(gsl_complex);
   *cnew = gsl_complex_pow_real(*c, a);
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_exp(int argc, VALUE *argv, VALUE obj)
@@ -566,7 +564,7 @@ static VALUE rb_gsl_complex_log_b(int argc, VALUE *argv, VALUE obj)
   case T_MODULE:
   case T_CLASS:
   case T_OBJECT:
-    if (argc != 2) 
+    if (argc != 2)
       rb_raise(rb_eArgError, "wrong number of arguments (%d for 2)", argc);
     switch (TYPE(argv[1])) {
     case T_ARRAY:
@@ -592,7 +590,7 @@ static VALUE rb_gsl_complex_log_b(int argc, VALUE *argv, VALUE obj)
 
     break;
   default:
-    if (argc != 1) 
+    if (argc != 1)
       rb_raise(rb_eArgError, "wrong number of arguments (%d for 1)", argc);
     CHECK_COMPLEX(obj);
     Data_Get_Struct(obj, gsl_complex, c);
@@ -610,7 +608,7 @@ static VALUE rb_gsl_complex_log_b(int argc, VALUE *argv, VALUE obj)
   }
   cnew = ALLOC(gsl_complex);
   *cnew = gsl_complex_log_b(*c, *a);
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_sin(int argc, VALUE *argv, VALUE obj)
@@ -658,7 +656,7 @@ static VALUE rb_gsl_complex_arcsin_real(VALUE obj, VALUE xx)
   tmp = gsl_complex_arcsin_real(x);
   cnew = ALLOC(gsl_complex);
   *cnew = tmp;
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_arccos(int argc, VALUE *argv, VALUE obj)
@@ -675,7 +673,7 @@ static VALUE rb_gsl_complex_arccos_real(VALUE obj, VALUE xx)
   tmp = gsl_complex_arccos_real(x);
   cnew = ALLOC(gsl_complex);
   *cnew = tmp;
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_arctan(int argc, VALUE *argv, VALUE obj)
@@ -697,7 +695,7 @@ static VALUE rb_gsl_complex_arcsec_real(VALUE obj, VALUE xx)
   tmp = gsl_complex_arcsec_real(x);
   cnew = ALLOC(gsl_complex);
   *cnew = tmp;
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_arccsc(int argc, VALUE *argv, VALUE obj)
@@ -714,7 +712,7 @@ static VALUE rb_gsl_complex_arccsc_real(VALUE obj, VALUE xx)
   tmp = gsl_complex_arccsc_real(x);
   cnew = ALLOC(gsl_complex);
   *cnew = tmp;
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_arccot(int argc, VALUE *argv, VALUE obj)
@@ -766,7 +764,7 @@ static VALUE rb_gsl_complex_arccosh_real(VALUE obj, VALUE xx)
   tmp = gsl_complex_arccosh_real(x);
   cnew = ALLOC(gsl_complex);
   *cnew = tmp;
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_arccosh(int argc, VALUE *argv, VALUE obj)
@@ -788,7 +786,7 @@ static VALUE rb_gsl_complex_arctanh_real(VALUE obj, VALUE xx)
   tmp = gsl_complex_arctanh_real(x);
   cnew = ALLOC(gsl_complex);
   *cnew = tmp;
-  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew); 
+  return Data_Wrap_Struct(cgsl_complex, 0, free, cnew);
 }
 
 static VALUE rb_gsl_complex_arcsech(int argc, VALUE *argv, VALUE obj)
@@ -807,10 +805,10 @@ static VALUE rb_gsl_complex_arccoth(int argc, VALUE *argv, VALUE obj)
 }
 
 int rbgsl_complex_equal(const gsl_complex *z1, const gsl_complex *z2,
-      double eps)
+                        double eps)
 {
   if (gsl_fcmp(z1->dat[0], z2->dat[0], eps) != 0) return 0;
-  if (gsl_fcmp(z1->dat[1], z2->dat[1], eps) != 0) return 0;  
+  if (gsl_fcmp(z1->dat[1], z2->dat[1], eps) != 0) return 0;
   return 1;
 }
 
@@ -825,11 +823,11 @@ static VALUE rb_gsl_complex_equal(int argc, VALUE *argv, VALUE obj)
 {
   gsl_complex *z1 = NULL, *z2 = NULL;
   double eps = 1e-8;
-  
+
   switch (argc) {
   case 1:
     CHECK_COMPLEX(argv[0]);
-    Data_Get_Struct(argv[0], gsl_complex, z2);      
+    Data_Get_Struct(argv[0], gsl_complex, z2);
     break;
   case 2:
     CHECK_COMPLEX(argv[0]);