gsl 1.15.3 → 1.16.0.6

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

@@ -10,7 +10,7 @@
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_dilog(VALUE obj, VALUE x)
 {

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

@@ -10,7 +10,7 @@
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_multiply_e(VALUE obj, VALUE x, VALUE y)
 {
@@ -25,7 +25,7 @@ static VALUE rb_gsl_sf_multiply_e(VALUE obj, VALUE x, VALUE y)
 }
 
 static VALUE rb_gsl_sf_multiply_err_e(VALUE obj, VALUE x, VALUE dx,
-				      VALUE y, VALUE dy)
+                                      VALUE y, VALUE dy)
 {
   gsl_sf_result *r;
   VALUE v;
@@ -35,7 +35,7 @@ static VALUE rb_gsl_sf_multiply_err_e(VALUE obj, VALUE x, VALUE dx,
   Need_Float(dx); Need_Float(dy);
   v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r);
   /*status =*/ gsl_sf_multiply_err_e(NUM2DBL(x), NUM2DBL(dx),
-				 NUM2DBL(y), NUM2DBL(dy), r);
+                                     NUM2DBL(y), NUM2DBL(dy), r);
   return v;
 }
 

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

@@ -10,7 +10,7 @@
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_ellint_Kcomp(int argc, VALUE *argv, VALUE obj)
 {
@@ -36,10 +36,10 @@ static VALUE rb_gsl_sf_ellint_Ecomp_e(VALUE obj, VALUE k, VALUE m)
 
 static VALUE rb_gsl_sf_ellint_F(int argc, VALUE *argv, VALUE obj)
 {
-  if (argc == 2) 
-    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_F, argv[0], argv[1], 
-				    INT2FIX(GSL_PREC_DOUBLE));
-  else 
+  if (argc == 2)
+    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_F, argv[0], argv[1],
+                                    INT2FIX(GSL_PREC_DOUBLE));
+  else
     return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_F, argv[0], argv[1], argv[2]);
 }
 
@@ -50,10 +50,10 @@ static VALUE rb_gsl_sf_ellint_F_e(VALUE obj, VALUE phi, VALUE k, VALUE m)
 
 static VALUE rb_gsl_sf_ellint_E(int argc, VALUE *argv, VALUE obj)
 {
-  if (argc == 2) 
-    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_E, argv[0], argv[1], 
-				    INT2FIX(GSL_PREC_DOUBLE));
-  else 
+  if (argc == 2)
+    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_E, argv[0], argv[1],
+                                    INT2FIX(GSL_PREC_DOUBLE));
+  else
     return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_E, argv[0], argv[1], argv[2]);
 }
 
@@ -65,42 +65,42 @@ static VALUE rb_gsl_sf_ellint_E_e(VALUE obj, VALUE phi, VALUE k, VALUE m)
 static VALUE rb_gsl_sf_ellint_P(int argc, VALUE *argv, VALUE obj)
 
 {
-  if (argc == 3) 
+  if (argc == 3)
     return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_P, argv[0], argv[1], argv[2],
-				    INT2FIX(GSL_PREC_DOUBLE));
-  else 
+                                    INT2FIX(GSL_PREC_DOUBLE));
+  else
     return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_P, argv[0], argv[1], argv[2],
-				    argv[3]);
+                                    argv[3]);
 }
 
-static VALUE rb_gsl_sf_ellint_P_e(VALUE obj, VALUE phi, VALUE k, 
-				  VALUE n, VALUE m)
+static VALUE rb_gsl_sf_ellint_P_e(VALUE obj, VALUE phi, VALUE k,
+                                  VALUE n, VALUE m)
 {
   return rb_gsl_sf_eval_e_double3_m(gsl_sf_ellint_P_e, phi, k, n, m);
 }
 
 static VALUE rb_gsl_sf_ellint_D(int argc, VALUE *argv, VALUE obj)
 {
-  if (argc == 3) 
+  if (argc == 3)
     return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_D, argv[0], argv[1], argv[2],
-				    INT2FIX(GSL_PREC_DOUBLE));
-  else 
+                                    INT2FIX(GSL_PREC_DOUBLE));
+  else
     return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_D, argv[0], argv[1], argv[2],
-				    argv[3]);
+                                    argv[3]);
 }
 
-static VALUE rb_gsl_sf_ellint_D_e(VALUE obj, VALUE phi, VALUE k, 
-				  VALUE n, VALUE m)
+static VALUE rb_gsl_sf_ellint_D_e(VALUE obj, VALUE phi, VALUE k,
+                                  VALUE n, VALUE m)
 {
   return rb_gsl_sf_eval_e_double3_m(gsl_sf_ellint_D_e, phi, k, n, m);
 }
 
 static VALUE rb_gsl_sf_ellint_RC(int argc, VALUE *argv, VALUE obj)
 {
-  if (argc == 2) 
-    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_RC, argv[0], argv[1], 
-				    INT2FIX(GSL_PREC_DOUBLE));
-  else 
+  if (argc == 2)
+    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_RC, argv[0], argv[1],
+                                    INT2FIX(GSL_PREC_DOUBLE));
+  else
     return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_RC, argv[0], argv[1], argv[2]);
 }
 
@@ -111,48 +111,48 @@ static VALUE rb_gsl_sf_ellint_RC_e(VALUE obj, VALUE x, VALUE y, VALUE m)
 
 static VALUE rb_gsl_sf_ellint_RD(int argc, VALUE *argv, VALUE obj)
 {
-  if (argc == 3) 
+  if (argc == 3)
     return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_RD, argv[0], argv[1], argv[2],
-				    INT2FIX(GSL_PREC_DOUBLE));
-  else 
+                                    INT2FIX(GSL_PREC_DOUBLE));
+  else
     return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_RD, argv[0], argv[1], argv[2],
-				    argv[3]);
+                                    argv[3]);
 }
 
-static VALUE rb_gsl_sf_ellint_RD_e(VALUE obj, VALUE x, VALUE y, 
-				   VALUE z, VALUE m)
+static VALUE rb_gsl_sf_ellint_RD_e(VALUE obj, VALUE x, VALUE y,
+                                   VALUE z, VALUE m)
 {
   return rb_gsl_sf_eval_e_double3_m(gsl_sf_ellint_RD_e, x, y, z, m);
 }
 
 static VALUE rb_gsl_sf_ellint_RF(int argc, VALUE *argv, VALUE obj)
 {
-  if (argc == 3) 
+  if (argc == 3)
     return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_RF, argv[0], argv[1], argv[2],
-				    INT2FIX(GSL_PREC_DOUBLE));
-  else 
+                                    INT2FIX(GSL_PREC_DOUBLE));
+  else
     return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_RF, argv[0], argv[1], argv[2],
-				    argv[3]);
+                                    argv[3]);
 }
 
-static VALUE rb_gsl_sf_ellint_RF_e(VALUE obj, VALUE x, VALUE y, 
-				   VALUE z, VALUE m)
+static VALUE rb_gsl_sf_ellint_RF_e(VALUE obj, VALUE x, VALUE y,
+                                   VALUE z, VALUE m)
 {
   return rb_gsl_sf_eval_e_double3_m(gsl_sf_ellint_RF_e, x, y, z, m);
 }
 
 static VALUE rb_gsl_sf_ellint_RJ(int argc, VALUE *argv, VALUE obj)
 {
-  if (argc == 4) 
+  if (argc == 4)
     return rb_gsl_sf_eval_double4_m(gsl_sf_ellint_RJ, argv[0], argv[1], argv[2],
-				    argv[3], INT2FIX(GSL_PREC_DOUBLE));
-  else 
+                                    argv[3], INT2FIX(GSL_PREC_DOUBLE));
+  else
     return rb_gsl_sf_eval_double4_m(gsl_sf_ellint_RJ, argv[0], argv[1], argv[2],
-				    argv[3], argv[4]);
+                                    argv[3], argv[4]);
 }
 
-static VALUE rb_gsl_sf_ellint_RJ_e(VALUE obj, VALUE x, VALUE y, 
-				   VALUE z, VALUE p, VALUE m)
+static VALUE rb_gsl_sf_ellint_RJ_e(VALUE obj, VALUE x, VALUE y,
+                                   VALUE z, VALUE p, VALUE m)
 {
   return rb_gsl_sf_eval_e_double4_m(gsl_sf_ellint_RJ_e, x, y, z, p, m);
 }

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

@@ -10,7 +10,7 @@
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_elljac_e(VALUE obj, VALUE n, VALUE m)
 {
@@ -19,8 +19,8 @@ static VALUE rb_gsl_sf_elljac_e(VALUE obj, VALUE n, VALUE m)
   //int status;
   Need_Float(n); Need_Float(m);
   /*status =*/ gsl_sf_elljac_e(NUM2DBL(n), NUM2DBL(m), &sn, &cn, &dn);
-  return rb_ary_new3(3, rb_float_new(sn), 
-		     rb_float_new(cn), rb_float_new(dn));
+  return rb_ary_new3(3, rb_float_new(sn),
+                     rb_float_new(cn), rb_float_new(dn));
 }
 
 void Init_gsl_sf_elljac(VALUE module)

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

@@ -10,7 +10,7 @@
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_erf(VALUE obj, VALUE x)
 {
@@ -62,7 +62,6 @@ static VALUE rb_gsl_sf_erf_Q_e(VALUE obj, VALUE x)
   return rb_gsl_sf_eval_e(gsl_sf_erf_Q_e, x);
 }
 
-#ifdef GSL_1_4_LATER
 static VALUE rb_gsl_sf_hazard(VALUE obj, VALUE x)
 {
   return rb_gsl_sf_eval1(gsl_sf_hazard, x);
@@ -72,7 +71,6 @@ static VALUE rb_gsl_sf_hazard_e(VALUE obj, VALUE x)
 {
   return rb_gsl_sf_eval_e(gsl_sf_hazard_e, x);
 }
-#endif
 
 void Init_gsl_sf_erfc(VALUE module)
 {
@@ -86,8 +84,6 @@ void Init_gsl_sf_erfc(VALUE module)
   rb_define_module_function(module, "erf_Z_e",  rb_gsl_sf_erf_Z_e, 1);
   rb_define_module_function(module, "erf_Q",  rb_gsl_sf_erf_Q, 1);
   rb_define_module_function(module, "erf_Q_e",  rb_gsl_sf_erf_Q_e, 1);
-#ifdef GSL_1_4_LATER
   rb_define_module_function(module, "hazard",  rb_gsl_sf_hazard, 1);
   rb_define_module_function(module, "hazard_e",  rb_gsl_sf_hazard_e, 1);
-#endif
 }

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

@@ -10,7 +10,7 @@
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_exp(VALUE obj, VALUE x)
 {
@@ -119,7 +119,7 @@ static VALUE rb_gsl_sf_exp_err_e10_e(VALUE obj, VALUE x, VALUE dx)
 }
 
 static VALUE rb_gsl_sf_exp_mult_err_e(VALUE obj, VALUE x, VALUE dx,
-				      VALUE y, VALUE dy)
+                                      VALUE y, VALUE dy)
 {
   gsl_sf_result *rslt = NULL;
   VALUE v;
@@ -133,7 +133,7 @@ static VALUE rb_gsl_sf_exp_mult_err_e(VALUE obj, VALUE x, VALUE dx,
 }
 
 static VALUE rb_gsl_sf_exp_mult_err_e10_e(VALUE obj, VALUE x, VALUE dx,
-					  VALUE y, VALUE dy)
+                                          VALUE y, VALUE dy)
 {
   gsl_sf_result_e10 *rslt = NULL;
   VALUE v;

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

@@ -10,7 +10,7 @@
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_expint_E1(VALUE obj, VALUE x)
 {
@@ -102,7 +102,6 @@ static VALUE rb_gsl_sf_atanint_e(VALUE obj, VALUE x)
   return rb_gsl_sf_eval_e(gsl_sf_atanint_e, x);
 }
 
-#ifdef GSL_1_3_LATER
 static VALUE rb_gsl_sf_expint_E1_scaled(VALUE obj, VALUE x)
 {
   return rb_gsl_sf_eval1(gsl_sf_expint_E1_scaled, x);
@@ -132,9 +131,7 @@ static VALUE rb_gsl_sf_expint_Ei_scaled_e(VALUE obj, VALUE x)
 {
   return rb_gsl_sf_eval_e(gsl_sf_expint_Ei_scaled_e, x);
 }
-#endif
 
-#ifdef GSL_1_10_LATER
 static VALUE rb_gsl_sf_expint_En(VALUE obj, VALUE n, VALUE x)
 {
   return rb_gsl_sf_eval_int_double(gsl_sf_expint_En, n, x);
@@ -148,8 +145,6 @@ static VALUE rb_gsl_sf_expint_En_e(VALUE obj, VALUE n, VALUE x)
   return val;
 }
 
-#endif
-
 void Init_gsl_sf_expint(VALUE module)
 {
   VALUE mgsl_sf_expint;
@@ -185,7 +180,6 @@ void Init_gsl_sf_expint(VALUE module)
   rb_define_module_function(mgsl_sf_expint, "three",  rb_gsl_sf_expint_3, 1);
   rb_define_module_function(mgsl_sf_expint, "three_e",  rb_gsl_sf_expint_3_e, 1);
 
-#ifdef GSL_1_3_LATER
   rb_define_module_function(module, "expint_E1_scaled",  rb_gsl_sf_expint_E1_scaled, 1);
   rb_define_module_function(module, "expint_E1_scaled_e",  rb_gsl_sf_expint_E1_scaled_e, 1);
   rb_define_module_function(module, "expint_E2_scaled",  rb_gsl_sf_expint_E2_scaled, 1);
@@ -199,13 +193,9 @@ void Init_gsl_sf_expint(VALUE module)
   rb_define_module_function(mgsl_sf_expint, "E2_scaled_e",  rb_gsl_sf_expint_E2_scaled_e, 1);
   rb_define_module_function(mgsl_sf_expint, "Ei_scaled",  rb_gsl_sf_expint_Ei_scaled, 1);
   rb_define_module_function(mgsl_sf_expint, "Ei_scaled_e",  rb_gsl_sf_expint_Ei_scaled_e, 1);
-#endif
 
-#ifdef GSL_1_10_LATER
   rb_define_module_function(module, "expint_En",  rb_gsl_sf_expint_En, 2);
   rb_define_module_function(mgsl_sf_expint, "En",  rb_gsl_sf_expint_En, 2);
   rb_define_module_function(module, "expint_En_e",  rb_gsl_sf_expint_En_e, 2);
-  rb_define_module_function(mgsl_sf_expint, "En_e",  rb_gsl_sf_expint_En_e, 2);  
-#endif
-
+  rb_define_module_function(mgsl_sf_expint, "En_e",  rb_gsl_sf_expint_En_e, 2);
 }

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

@@ -10,7 +10,7 @@
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_fermi_dirac_m1(VALUE obj, VALUE x)
 {

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

@@ -9,7 +9,7 @@
   WITHOUT ANY WARRANTY.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_gamma(VALUE obj, VALUE x)
 {
@@ -242,12 +242,10 @@ static VALUE rb_gsl_sf_gamma_inc(VALUE obj, VALUE a, VALUE x)
   return rb_gsl_sf_eval_double_double(gsl_sf_gamma_inc_P, a, x);
 }
 
-#ifdef GSL_1_4_LATER
 static VALUE rb_gsl_sf_gamma_inc_e(VALUE obj, VALUE a, VALUE x)
 {
   return rb_gsl_sf_eval_e_double2(gsl_sf_gamma_inc_e, a, x);
 }
-#endif
 
 static VALUE rb_gsl_sf_beta(VALUE obj, VALUE a, VALUE b)
 {
@@ -284,7 +282,7 @@ static VALUE rb_gsl_sf_beta_inc_e(VALUE obj, VALUE a, VALUE b, VALUE x)
 double mygsl_binomial_coef(unsigned int n, unsigned int k);
 double mygsl_binomial_coef(unsigned int n, unsigned int k)
 {
-  return floor(0.5 + exp(gsl_sf_lnfact(n) - gsl_sf_lnfact(k) - gsl_sf_lnfact(n-k))); 
+  return floor(0.5 + exp(gsl_sf_lnfact(n) - gsl_sf_lnfact(k) - gsl_sf_lnfact(n-k)));
 }
 
 static VALUE rb_gsl_sf_bincoef(VALUE obj, VALUE n, VALUE k)
@@ -333,9 +331,7 @@ void Init_gsl_sf_gamma(VALUE module)
   rb_define_module_function(module, "gamma_inc_Q_e",  rb_gsl_sf_gamma_inc_Q_e, 2);
   rb_define_module_function(module, "gamma_inc",  rb_gsl_sf_gamma_inc, 2);
 
-#ifdef GSL_1_4_LATER
   rb_define_module_function(module, "gamma_inc_e",  rb_gsl_sf_gamma_inc_e, 2);
-#endif
   rb_define_module_function(module, "beta",  rb_gsl_sf_beta, 2);
   rb_define_module_function(module, "beta_e",  rb_gsl_sf_beta_e, 2);
   rb_define_module_function(module, "lnbeta",  rb_gsl_sf_lnbeta, 2);

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

@@ -9,7 +9,7 @@
   WITHOUT ANY WARRANTY.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_gegenpoly_1(VALUE obj, VALUE lambda, VALUE x)
 {

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

@@ -9,7 +9,7 @@
   WITHOUT ANY WARRANTY.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 /* Checked with Mathematica. */
 static VALUE rb_gsl_sf_hyperg_0F1(VALUE obj, VALUE c, VALUE x)

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

@@ -9,10 +9,10 @@
   WITHOUT ANY WARRANTY.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_laguerre_X(int argc, VALUE *argv, VALUE obj,
-				  double (*f)(double, double))
+                                  double (*f)(double, double))
 {
   switch (argc) {
   case 2:
@@ -62,11 +62,11 @@ static VALUE rb_gsl_sf_laguerre_n(int argc, VALUE *argv, VALUE obj)
   switch (argc) {
   case 3:
     return rb_gsl_sf_eval_int_double_double(gsl_sf_laguerre_n, argv[0],
-					    argv[1], argv[2]);
+                                            argv[1], argv[2]);
     break;
   case 2:
     return rb_gsl_sf_eval_int_double_double(gsl_sf_laguerre_n, argv[0],
-					    INT2FIX(0), argv[1]);
+                                            INT2FIX(0), argv[1]);
     break;
   default:
     rb_raise(rb_eArgError, "wrong number of arguments (%d for 1 or 2)", argc);

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

@@ -9,7 +9,7 @@
   WITHOUT ANY WARRANTY.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 
 static VALUE rb_gsl_sf_lambert_W0(VALUE obj, VALUE x)
 {

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

@@ -9,7 +9,7 @@
   WITHOUT ANY WARRANTY.
 */
 
-#include "rb_gsl_sf.h"
+#include "include/rb_gsl_sf.h"
 EXTERN VALUE cgsl_vector;
 
 static VALUE rb_gsl_sf_legendre_P1(VALUE obj, VALUE x)