gsl 1.15.3 → 1.16.0.6

data/ext/gsl_native/bspline.c

@@ -0,0 +1,122 @@
+#include "include/rb_gsl.h"
+#include "gsl/gsl_bspline.h"
+
+static VALUE cBSWS;
+
+static VALUE rb_gsl_bspline_alloc(VALUE klass, VALUE k, VALUE n)
+{
+  gsl_bspline_workspace *w;
+  w = gsl_bspline_alloc(FIX2INT(k), FIX2INT(n));
+  return Data_Wrap_Struct(klass, 0, gsl_bspline_free, w);
+}
+
+static VALUE rb_gsl_bspline_ncoeffs(VALUE obj)
+{
+  gsl_bspline_workspace *w;
+  Data_Get_Struct(obj, gsl_bspline_workspace, w);
+  return INT2FIX((int)gsl_bspline_ncoeffs(w));
+}
+static VALUE rb_gsl_bspline_order(VALUE obj)
+{
+  gsl_bspline_workspace *w;
+  Data_Get_Struct(obj, gsl_bspline_workspace, w);
+  return INT2FIX((int)gsl_bspline_order(w));
+}
+static VALUE rb_gsl_bspline_nbreak(VALUE obj)
+{
+  gsl_bspline_workspace *w;
+  Data_Get_Struct(obj, gsl_bspline_workspace, w);
+  return INT2FIX((int)gsl_bspline_nbreak(w));
+}
+static VALUE rb_gsl_bspline_breakpoint(VALUE obj, VALUE i)
+{
+  gsl_bspline_workspace *w;
+  Data_Get_Struct(obj, gsl_bspline_workspace, w);
+  return rb_float_new(gsl_bspline_breakpoint(FIX2INT(i), w));
+}
+static VALUE rb_gsl_bspline_knots(VALUE obj, VALUE b)
+{
+  gsl_bspline_workspace *w;
+  gsl_vector *bpts;
+  CHECK_VECTOR(b);
+  Data_Get_Struct(obj, gsl_bspline_workspace, w);
+  Data_Get_Struct(b, gsl_vector, bpts);
+  gsl_bspline_knots(bpts, w);
+  return Data_Wrap_Struct(cgsl_vector_view_ro, 0, NULL, w->knots);
+}
+static VALUE rb_gsl_bspline_knots_uniform(int argc, VALUE *argv, VALUE obj)
+{
+  gsl_bspline_workspace *w;
+  int argc2;
+  switch (TYPE(obj)) {
+  case T_MODULE:
+  case T_CLASS:
+  case T_OBJECT:
+    if (!rb_obj_is_kind_of(argv[argc-1], cBSWS)) {
+      rb_raise(rb_eTypeError, "Wrong argument type %s (GSL::BSpline expected)",
+               rb_class2name(CLASS_OF(argv[argc-1])));
+    }
+    Data_Get_Struct(argv[argc-1], gsl_bspline_workspace, w);
+    argc2 = argc-1;
+    break;
+  default:
+    Data_Get_Struct(obj, gsl_bspline_workspace, w);
+    argc2 = argc;
+  }
+  if (argc2 != 2) rb_raise(rb_eArgError, "Wrong number of arguments.");
+  gsl_bspline_knots_uniform(NUM2DBL(argv[0]), NUM2DBL(argv[1]), w);
+  return Data_Wrap_Struct(cgsl_vector_view_ro, 0, NULL, w->knots);
+}
+static VALUE rb_gsl_bspline_eval(int argc, VALUE *argv, VALUE obj)
+{
+  gsl_bspline_workspace *w;
+  double x;
+  gsl_vector *B;
+  VALUE vB;
+
+  Data_Get_Struct(obj, gsl_bspline_workspace, w);
+
+  switch (argc) {
+  case 2:
+    CHECK_VECTOR(argv[1]);
+    Data_Get_Struct(argv[1], gsl_vector, B);
+    vB = argv[1];
+    x = NUM2DBL(argv[0]);
+    break;
+  case 1:
+    x = NUM2DBL(argv[0]);
+    B = gsl_vector_alloc(w->nbreak+w->k-2);
+    vB = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, B);
+    break;
+  default:
+    rb_raise(rb_eArgError, "Wrong number of arguments (%d for 1 or 2)", argc);
+  }
+
+  gsl_bspline_eval(x, B, w);
+
+  return vB;
+}
+
+static VALUE rb_gsl_bspline_greville_abscissa(VALUE obj, VALUE i)
+{
+  gsl_bspline_workspace *w;
+  Data_Get_Struct(obj, gsl_bspline_workspace, w);
+  return rb_float_new(gsl_bspline_greville_abscissa(i, w));
+}
+
+void Init_bspline(VALUE module)
+{
+  cBSWS = rb_define_class_under(module, "BSpline", cGSL_Object);
+
+  rb_define_singleton_method(cBSWS, "alloc", rb_gsl_bspline_alloc, 2);
+
+  rb_define_method(cBSWS, "ncoeffs", rb_gsl_bspline_ncoeffs, 0);
+  rb_define_method(cBSWS, "order", rb_gsl_bspline_order, 0);
+  rb_define_method(cBSWS, "nbreak", rb_gsl_bspline_nbreak, 0);
+  rb_define_method(cBSWS, "breakpoint", rb_gsl_bspline_breakpoint, 1);
+  rb_define_method(cBSWS, "knots", rb_gsl_bspline_knots, 1);
+  rb_define_method(cBSWS, "knots_uniform", rb_gsl_bspline_knots_uniform, -1);
+  rb_define_singleton_method(cBSWS, "knots_uniform", rb_gsl_bspline_knots_uniform, -1);
+  rb_define_method(cBSWS, "eval", rb_gsl_bspline_eval, -1);
+  rb_define_method(cBSWS, "greville_abscissa", rb_gsl_bspline_greville_abscissa, 1);
+}

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

@@ -9,21 +9,18 @@
   WITHOUT ANY WARRANTY.
 */
 
-#include "rb_gsl_config.h"
-
-#ifdef GSL_1_4_LATER
-#include "rb_gsl_array.h"
-#include "rb_gsl_common.h"
-#include "rb_gsl_rng.h"
+#include "include/rb_gsl_array.h"
+#include "include/rb_gsl_common.h"
+#include "include/rb_gsl_rng.h"
 #include <gsl/gsl_randist.h>
 #include <gsl/gsl_cdf.h>
 
 VALUE rb_gsl_eval_pdf_cdf(VALUE xx, double (*f)(double));
 VALUE rb_gsl_eval_pdf_cdf2(VALUE xx, VALUE aa, double (*f)(double, double));
-VALUE rb_gsl_eval_pdf_cdf3(VALUE xx, VALUE aa, VALUE bb, 
-			      double (*f)(double, double, double));
-VALUE rb_gsl_eval_pdf_cdf2_uint(VALUE xx, VALUE aa, 
-				       double (*f)(unsigned int, double));
+VALUE rb_gsl_eval_pdf_cdf3(VALUE xx, VALUE aa, VALUE bb,
+                           double (*f)(double, double, double));
+VALUE rb_gsl_eval_pdf_cdf2_uint(VALUE xx, VALUE aa,
+                                double (*f)(unsigned int, double));
 
 
 static VALUE rb_gsl_cdf_gaussian_P(int argc, VALUE *argv, VALUE obj)
@@ -72,316 +69,313 @@ static VALUE rb_gsl_cdf_gaussian_Qinv(int argc, VALUE *argv, VALUE obj)
 static VALUE rb_gsl_cdf_exponential_P(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_exponential_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_exponential_Q(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_exponential_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_exponential_Pinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_exponential_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_exponential_Qinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_exponential_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_laplace_P(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_laplace_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_laplace_Q(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_laplace_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_laplace_Pinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_laplace_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_laplace_Qinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_laplace_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_cauchy_P(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_cauchy_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_cauchy_Q(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_cauchy_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_cauchy_Pinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_cauchy_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_cauchy_Qinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_cauchy_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_rayleigh_P(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_rayleigh_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_rayleigh_Q(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_rayleigh_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_rayleigh_Pinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_rayleigh_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_rayleigh_Qinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_rayleigh_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_gamma_P(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gamma_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_gamma_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_gamma_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_gamma_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_gamma_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_gamma_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_gamma_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_flat_P(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_flat_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_flat_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_flat_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_flat_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_flat_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_flat_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_flat_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_lognormal_P(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_lognormal_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_lognormal_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_lognormal_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_lognormal_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_lognormal_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_lognormal_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b,  gsl_cdf_lognormal_Qinv);
-}  
+}
 
-#ifdef GSL_1_6_LATER
 static VALUE rb_gsl_cdf_exppow_P(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_exppow_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_exppow_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_exppow_Q);
-}  
-#endif
+}
 
 static VALUE rb_gsl_cdf_chisq_P(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_chisq_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_chisq_Q(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_chisq_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_chisq_Pinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_chisq_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_chisq_Qinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_chisq_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_fdist_P(VALUE obj, VALUE x, VALUE nu1, VALUE nu2)
 {
   return rb_gsl_eval_pdf_cdf3(x, nu1, nu2, gsl_cdf_fdist_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_fdist_Q(VALUE obj, VALUE x, VALUE nu1, VALUE nu2)
 {
   return rb_gsl_eval_pdf_cdf3(x, nu1, nu2, gsl_cdf_fdist_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_tdist_P(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_tdist_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_tdist_Q(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_tdist_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_tdist_Pinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_tdist_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_tdist_Qinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_tdist_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_beta_P(VALUE obj, VALUE x, VALUE nu1, VALUE nu2)
 {
   return rb_gsl_eval_pdf_cdf3(x, nu1, nu2, gsl_cdf_beta_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_beta_Q(VALUE obj, VALUE x, VALUE nu1, VALUE nu2)
 {
   return rb_gsl_eval_pdf_cdf3(x, nu1, nu2, gsl_cdf_beta_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_logistic_P(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_logistic_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_logistic_Q(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_logistic_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_logistic_Pinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_logistic_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_logistic_Qinv(VALUE obj, VALUE x, VALUE mu)
 {
   return rb_gsl_eval_pdf_cdf2(x, mu, gsl_cdf_logistic_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_pareto_P(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_pareto_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_pareto_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_pareto_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_pareto_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_pareto_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_pareto_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_pareto_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_weibull_P(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_weibull_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_weibull_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_weibull_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_weibull_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_weibull_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_weibull_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_weibull_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_gumbel1_P(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel1_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_gumbel1_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel1_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_gumbel1_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel1_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_gumbel1_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel1_Qinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_gumbel2_P(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel2_P);
-}  
+}
 
 static VALUE rb_gsl_cdf_gumbel2_Q(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel2_Q);
-}  
+}
 
 static VALUE rb_gsl_cdf_gumbel2_Pinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel2_Pinv);
-}  
+}
 
 static VALUE rb_gsl_cdf_gumbel2_Qinv(VALUE obj, VALUE x, VALUE a, VALUE b)
 {
   return rb_gsl_eval_pdf_cdf3(x, a, b, gsl_cdf_gumbel2_Qinv);
-}  
+}
 
-#ifdef GSL_1_8_LATER
 static VALUE rb_gsl_cdf_binomial_P(VALUE obj, VALUE kk, VALUE pp, VALUE nn)
 {
   unsigned int k, n;
@@ -518,7 +512,6 @@ static VALUE rb_gsl_cdf_fdist_Qinv(VALUE obj, VALUE qq, VALUE aa, VALUE bb)
   b = NUM2DBL(bb);
   return rb_float_new(gsl_cdf_fdist_Qinv(Q, a, b));
 }
-#endif
 
 void Init_gsl_cdf(VALUE module)
 {
@@ -537,7 +530,7 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(module, "cdf_ugaussian_Qinv", rb_gsl_cdf_gaussian_Qinv, -1);
 
   rb_define_module_function(module, "cdf_exponential_P", rb_gsl_cdf_exponential_P, 2);
- rb_define_module_function(module, "cdf_exponential_Q", rb_gsl_cdf_exponential_Q, 2);
+  rb_define_module_function(module, "cdf_exponential_Q", rb_gsl_cdf_exponential_Q, 2);
   rb_define_module_function(module, "cdf_exponential_Pinv", rb_gsl_cdf_exponential_Pinv, 2);
   rb_define_module_function(module, "cdf_exponential_Qinv", rb_gsl_cdf_exponential_Qinv, 2);
 
@@ -575,7 +568,7 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(module, "cdf_chisq_Q", rb_gsl_cdf_chisq_Q, 2);
   rb_define_module_function(module, "cdf_chisq_Pinv", rb_gsl_cdf_chisq_Pinv, 2);
   rb_define_module_function(module, "cdf_chisq_Qinv", rb_gsl_cdf_chisq_Qinv, 2);
-  
+
   rb_define_module_function(module, "cdf_fdist_P", rb_gsl_cdf_fdist_P, 3);
   rb_define_module_function(module, "cdf_fdist_Q", rb_gsl_cdf_fdist_Q, 3);
 
@@ -583,7 +576,7 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(module, "cdf_tdist_Q", rb_gsl_cdf_tdist_Q, 2);
   rb_define_module_function(module, "cdf_tdist_Pinv", rb_gsl_cdf_tdist_Pinv, 2);
   rb_define_module_function(module, "cdf_tdist_Qinv", rb_gsl_cdf_tdist_Qinv, 2);
-  
+
   rb_define_module_function(module, "cdf_beta_P", rb_gsl_cdf_beta_P, 3);
   rb_define_module_function(module, "cdf_beta_Q", rb_gsl_cdf_beta_Q, 3);
 
@@ -591,7 +584,7 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(module, "cdf_logistic_Q", rb_gsl_cdf_logistic_Q, 2);
   rb_define_module_function(module, "cdf_logistic_Pinv", rb_gsl_cdf_logistic_Pinv, 2);
   rb_define_module_function(module, "cdf_logistic_Qinv", rb_gsl_cdf_logistic_Qinv, 2);
-  
+
   rb_define_module_function(module, "cdf_pareto_P", rb_gsl_cdf_pareto_P, 3);
   rb_define_module_function(module, "cdf_pareto_Q", rb_gsl_cdf_pareto_Q, 3);
   rb_define_module_function(module, "cdf_pareto_Pinv", rb_gsl_cdf_pareto_Pinv, 3);
@@ -624,12 +617,12 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(mgsl_cdf, "ugaussian_Qinv", rb_gsl_cdf_gaussian_Qinv, -1);
 
   rb_define_module_function(mgsl_cdf, "exponential_P", rb_gsl_cdf_exponential_P, 2);
- rb_define_module_function(mgsl_cdf, "exponential_Q", rb_gsl_cdf_exponential_Q, 2);
+  rb_define_module_function(mgsl_cdf, "exponential_Q", rb_gsl_cdf_exponential_Q, 2);
   rb_define_module_function(mgsl_cdf, "exponential_Pinv", rb_gsl_cdf_exponential_Pinv, 2);
   rb_define_module_function(mgsl_cdf, "exponential_Qinv", rb_gsl_cdf_exponential_Qinv, 2);
 
- rb_define_module_function(mgsl_cdf, "laplace_P", rb_gsl_cdf_laplace_P, 2);
- rb_define_module_function(mgsl_cdf, "laplace_Q", rb_gsl_cdf_laplace_Q, 2);
+  rb_define_module_function(mgsl_cdf, "laplace_P", rb_gsl_cdf_laplace_P, 2);
+  rb_define_module_function(mgsl_cdf, "laplace_Q", rb_gsl_cdf_laplace_Q, 2);
   rb_define_module_function(mgsl_cdf, "laplace_Pinv", rb_gsl_cdf_laplace_Pinv, 2);
   rb_define_module_function(mgsl_cdf, "laplace_Qinv", rb_gsl_cdf_laplace_Qinv, 2);
 
@@ -662,7 +655,7 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(mgsl_cdf, "chisq_Q", rb_gsl_cdf_chisq_Q, 2);
   rb_define_module_function(mgsl_cdf, "chisq_Pinv", rb_gsl_cdf_chisq_Pinv, 2);
   rb_define_module_function(mgsl_cdf, "chisq_Qinv", rb_gsl_cdf_chisq_Qinv, 2);
-  
+
   rb_define_module_function(mgsl_cdf, "fdist_P", rb_gsl_cdf_fdist_P, 3);
   rb_define_module_function(mgsl_cdf, "fdist_Q", rb_gsl_cdf_fdist_Q, 3);
 
@@ -670,7 +663,7 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(mgsl_cdf, "tdist_Q", rb_gsl_cdf_tdist_Q, 2);
   rb_define_module_function(mgsl_cdf, "tdist_Pinv", rb_gsl_cdf_tdist_Pinv, 2);
   rb_define_module_function(mgsl_cdf, "tdist_Qinv", rb_gsl_cdf_tdist_Qinv, 2);
-  
+
   rb_define_module_function(mgsl_cdf, "beta_P", rb_gsl_cdf_beta_P, 3);
   rb_define_module_function(mgsl_cdf, "beta_Q", rb_gsl_cdf_beta_Q, 3);
 
@@ -678,7 +671,7 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(mgsl_cdf, "logistic_Q", rb_gsl_cdf_logistic_Q, 2);
   rb_define_module_function(mgsl_cdf, "logistic_Pinv", rb_gsl_cdf_logistic_Pinv, 2);
   rb_define_module_function(mgsl_cdf, "logistic_Qinv", rb_gsl_cdf_logistic_Qinv, 2);
-  
+
   rb_define_module_function(mgsl_cdf, "pareto_P", rb_gsl_cdf_pareto_P, 3);
   rb_define_module_function(mgsl_cdf, "pareto_Q", rb_gsl_cdf_pareto_Q, 3);
   rb_define_module_function(mgsl_cdf, "pareto_Pinv", rb_gsl_cdf_pareto_Pinv, 3);
@@ -699,16 +692,12 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(mgsl_cdf, "gumbel2_Pinv", rb_gsl_cdf_gumbel2_Pinv, 3);
   rb_define_module_function(mgsl_cdf, "gumbel2_Qinv", rb_gsl_cdf_gumbel2_Qinv, 3);
 
-
-#ifdef GSL_1_6_LATER
   rb_define_module_function(module, "cdf_exppow_P", rb_gsl_cdf_exppow_P, 3);
   rb_define_module_function(module, "cdf_exppow_Q", rb_gsl_cdf_exppow_Q, 3);
 
   rb_define_module_function(mgsl_cdf, "exppow_P", rb_gsl_cdf_exppow_P, 3);
   rb_define_module_function(mgsl_cdf, "exppow_Q", rb_gsl_cdf_exppow_Q, 3);
-#endif
 
-#ifdef GSL_1_8_LATER
   rb_define_module_function(module, "cdf_binomial_P", rb_gsl_cdf_binomial_P, 3);
   rb_define_module_function(module, "cdf_binomial_Q", rb_gsl_cdf_binomial_Q, 3);
   rb_define_module_function(mgsl_cdf, "binomial_P", rb_gsl_cdf_binomial_P, 3);
@@ -748,7 +737,4 @@ void Init_gsl_cdf(VALUE module)
   rb_define_module_function(mgsl_cdf, "fdist_Pinv", rb_gsl_cdf_fdist_Pinv, 3);
   rb_define_module_function(module, "cdf_fdist_Qinv", rb_gsl_cdf_fdist_Qinv, 3);
   rb_define_module_function(mgsl_cdf, "fdist_Qinv", rb_gsl_cdf_fdist_Qinv, 3);
-#endif
-
 }
-#endif