gsl 1.12.108

data/ext/sf_clausen.c

@@ -0,0 +1,28 @@
+/*
+  sf_clausen.c
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2006 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY.
+*/
+
+#include "rb_gsl_sf.h"
+
+static VALUE rb_gsl_sf_clausen(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_clausen, x);
+}
+
+static VALUE rb_gsl_sf_clausen_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_clausen_e, x);
+}
+
+void Init_gsl_sf_clausen(VALUE module)
+{
+  rb_define_module_function(module, "clausen",  rb_gsl_sf_clausen, 1);
+  rb_define_module_function(module, "clausen_e",  rb_gsl_sf_clausen_e, 1);
+}

data/ext/sf_coulomb.c

@@ -0,0 +1,206 @@
+/*
+  sf_coulomb.c
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2006 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#include "rb_gsl_sf.h"
+EXTERN VALUE cgsl_vector;
+
+static VALUE rb_gsl_sf_hydrogenicR_1(VALUE obj, VALUE Z, VALUE r)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_hydrogenicR_1, Z, r);
+}
+
+static VALUE rb_gsl_sf_hydrogenicR_1_e(VALUE obj,  VALUE Z, VALUE r)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_hydrogenicR_1_e, Z, r);
+}
+
+static VALUE rb_gsl_sf_hydrogenicR(VALUE obj, VALUE n, VALUE l, 
+				   VALUE Z, VALUE r)
+{
+  return rb_float_new(gsl_sf_hydrogenicR(FIX2INT(n), FIX2INT(l),
+					   NUM2DBL(Z), NUM2DBL(r)));
+}
+
+static VALUE rb_gsl_sf_hydrogenicR_e(VALUE obj, VALUE n, VALUE l, 
+				     VALUE Z, VALUE r)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  int status;
+  CHECK_FIXNUM(n); CHECK_FIXNUM(l);
+  Need_Float(Z); Need_Float(r);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  status = gsl_sf_hydrogenicR_e(FIX2INT(n), FIX2INT(l),
+			    NUM2DBL(Z), NUM2DBL(r), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_coulomb_wave_FG_e(VALUE obj, VALUE eta, VALUE x,
+					 VALUE L_F, VALUE k)
+
+{
+  gsl_sf_result *F, *Fp, *G, *Gp;
+  VALUE vF, vFp, vG, vGp;
+  double exp_G, exp_F;
+  int status;
+  Need_Float(eta); Need_Float(x); Need_Float(L_F);
+  CHECK_FIXNUM(k);
+  vF = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, F);
+  vFp = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, Fp);
+  vG = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, G);
+  vGp = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, Gp);
+  status = gsl_sf_coulomb_wave_FG_e(NUM2DBL(eta), NUM2DBL(x), NUM2DBL(L_F),
+				    FIX2INT(k), F, Fp, G, Gp, &exp_F, &exp_G);
+  return rb_ary_new3(7, vF, vFp, vG, vGp,
+		     rb_float_new(exp_F), rb_float_new(exp_G), INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_coulomb_wave_F_array(VALUE obj, VALUE Lmin, VALUE kmax,
+					 VALUE eta, VALUE x)
+{
+  double F_exponent;
+  int status;
+  size_t size;
+  gsl_vector *v = NULL;
+  CHECK_FIXNUM(kmax);
+  Need_Float(Lmin); Need_Float(eta); Need_Float(x);
+  size = FIX2INT(kmax);
+  v = gsl_vector_alloc(size);
+  status = gsl_sf_coulomb_wave_F_array(NUM2DBL(Lmin), size, NUM2DBL(eta), 
+				   NUM2DBL(x), v->data, &F_exponent);
+  
+  return rb_ary_new3(3, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v), 
+		     rb_float_new(F_exponent), INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_coulomb_wave_FG_array(VALUE obj, VALUE Lmin, VALUE kmax,
+					 VALUE eta, VALUE x)
+{
+  double F_exponent, G_exponent;
+  int status;
+  size_t size, stride;
+  gsl_vector *vf = NULL, *vg = NULL;
+  VALUE fary, gary;
+  CHECK_FIXNUM(kmax);
+  Need_Float(Lmin); Need_Float(eta); Need_Float(x);
+  size = FIX2INT(kmax);
+  vf = gsl_vector_alloc(size);
+  vg = gsl_vector_alloc(size);
+  stride = vf->stride;
+
+  status = gsl_sf_coulomb_wave_FG_array(NUM2DBL(Lmin), size, NUM2DBL(eta), 
+				   NUM2DBL(x), vf->data, vg->data,
+				    &F_exponent, &G_exponent);
+  fary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vf);
+  gary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vf);
+  return rb_ary_new3(5, fary, gary,
+		     rb_float_new(F_exponent), rb_float_new(G_exponent),
+		     INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_coulomb_wave_FGp_array(VALUE obj, VALUE Lmin, VALUE kmax,
+					 VALUE eta, VALUE x)
+{
+  double F_exponent, G_exponent;
+  int status;
+  size_t size, stride;
+  gsl_vector *vf = NULL, *vg = NULL, *vfp = NULL, *vgp = NULL;
+  VALUE fary, gary, fpary, gpary;
+  CHECK_FIXNUM(kmax);
+  Need_Float(Lmin); Need_Float(eta); Need_Float(x);
+  size = FIX2INT(kmax);
+  vf = gsl_vector_alloc(size);
+  vfp = gsl_vector_alloc(size);
+  vg = gsl_vector_alloc(size);
+  vgp = gsl_vector_alloc(size);
+  stride = vf->stride;
+
+  status = gsl_sf_coulomb_wave_FGp_array(NUM2DBL(Lmin), size, NUM2DBL(eta), 
+				     NUM2DBL(x), vf->data, vfp->data, 
+				     vg->data, vgp->data,
+				     &F_exponent, &G_exponent);
+  fary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vf);
+  fpary =Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vfp);
+  gary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vg);
+  gpary =Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vgp);
+  return rb_ary_new3(7, fary, fpary, gary, gpary,
+		     rb_float_new(F_exponent), rb_float_new(G_exponent),
+		     INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_coulomb_wave_sphF_array(VALUE obj, VALUE Lmin, VALUE kmax,
+					 VALUE eta, VALUE x)
+{
+  int status;
+  size_t size;
+  gsl_vector *v = NULL, *v2 = NULL;
+  CHECK_FIXNUM(kmax);
+  Need_Float(Lmin); Need_Float(eta); Need_Float(x);
+  size = FIX2INT(kmax);
+  v = gsl_vector_alloc(size);
+  v2 = gsl_vector_alloc(size);
+  status =  gsl_sf_coulomb_wave_sphF_array(NUM2DBL(Lmin), size, NUM2DBL(eta), 
+				       NUM2DBL(x), v->data, v2->data);
+  return rb_ary_new3(3, Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v), 
+		     Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v2), 
+		     INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_coulomb_CL_e(VALUE obj, VALUE L, VALUE eta)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_coulomb_CL_e, L, eta);
+}
+
+static VALUE rb_gsl_sf_coulomb_CL_array(VALUE obj, VALUE Lmin, VALUE kmax,
+					VALUE eta)
+{
+  gsl_vector *v = NULL;
+  size_t size;
+  int status;
+  CHECK_FIXNUM(kmax);
+  Need_Float(Lmin); Need_Float(eta); 
+  size = FIX2INT(kmax);
+  v = gsl_vector_alloc(size);
+  status = gsl_sf_coulomb_CL_array(NUM2DBL(Lmin), size, NUM2DBL(eta), v->data);
+  return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
+}
+
+void Init_gsl_sf_coulomb(VALUE module)
+{
+  VALUE mgsl_sf_coulomb;
+
+  rb_define_module_function(module, "hydrogenicR_1",  rb_gsl_sf_hydrogenicR_1, 2);
+  rb_define_module_function(module, "hydrogenicR_1_e",  rb_gsl_sf_hydrogenicR_1_e, 2);
+  rb_define_module_function(module, "hydrogenicR",  rb_gsl_sf_hydrogenicR, 4);
+  rb_define_module_function(module, "hydrogenicR_e",  rb_gsl_sf_hydrogenicR_e, 4);
+  rb_define_module_function(module, "coulomb_wave_FG_e",  rb_gsl_sf_coulomb_wave_FG_e, 4);
+  rb_define_module_function(module, "coulomb_wave_F_array",  rb_gsl_sf_coulomb_wave_F_array, 4);
+  rb_define_module_function(module, "coulomb_wave_FG_array",  rb_gsl_sf_coulomb_wave_FG_array, 4);
+  rb_define_module_function(module, "coulomb_wave_FGp_array",  rb_gsl_sf_coulomb_wave_FGp_array, 4);
+  rb_define_module_function(module, "coulomb_wave_sphF_array ",  rb_gsl_sf_coulomb_wave_sphF_array, 4);
+  rb_define_module_function(module, "coulomb_CL_e",  rb_gsl_sf_coulomb_CL_e, 2);
+  rb_define_module_function(module, "coulomb_CL_array",  rb_gsl_sf_coulomb_CL_array, 3);
+
+  mgsl_sf_coulomb = rb_define_module_under(module, "Coulomb");
+  
+  rb_define_module_function(mgsl_sf_coulomb, "hydrogenicR_1",  rb_gsl_sf_hydrogenicR_1, 2);
+  rb_define_module_function(mgsl_sf_coulomb, "hydrogenicR_1_e",  rb_gsl_sf_hydrogenicR_1_e, 2);
+  rb_define_module_function(mgsl_sf_coulomb, "hydrogenicR",  rb_gsl_sf_hydrogenicR, 4);
+  rb_define_module_function(mgsl_sf_coulomb, "hydrogenicR_e",  rb_gsl_sf_hydrogenicR_e, 4);
+  rb_define_module_function(mgsl_sf_coulomb, "wave_FG_e",  rb_gsl_sf_coulomb_wave_FG_e, 4);
+  rb_define_module_function(mgsl_sf_coulomb, "wave_F_array",  rb_gsl_sf_coulomb_wave_F_array, 4);
+  rb_define_module_function(mgsl_sf_coulomb, "wave_FG_array",  rb_gsl_sf_coulomb_wave_FG_array, 4);
+  rb_define_module_function(mgsl_sf_coulomb, "wave_FGp_array",  rb_gsl_sf_coulomb_wave_FGp_array, 4);
+  rb_define_module_function(mgsl_sf_coulomb, "wave_sphF_array ",  rb_gsl_sf_coulomb_wave_sphF_array, 4);
+  rb_define_module_function(mgsl_sf_coulomb, "CL_e",  rb_gsl_sf_coulomb_CL_e, 2);
+  rb_define_module_function(mgsl_sf_coulomb, "CL_array",  rb_gsl_sf_coulomb_CL_array, 3);
+}

data/ext/sf_coupling.c

@@ -0,0 +1,118 @@
+/*
+  sf_coupling.c
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2006 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#include "rb_gsl_sf.h"
+
+static VALUE rb_gsl_sf_coupling_3j(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_ma, VALUE two_mb, VALUE two_mc) 
+{
+  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
+  CHECK_FIXNUM(two_ma); CHECK_FIXNUM(two_mb); CHECK_FIXNUM(two_mc);
+  return rb_float_new(gsl_sf_coupling_3j(FIX2INT(two_ja), FIX2INT(two_jb),
+					 FIX2INT(two_jc), FIX2INT(two_ma),
+					 FIX2INT(two_mb), FIX2INT(two_mc)));
+}
+
+static VALUE rb_gsl_sf_coupling_3j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_ma, VALUE two_mb, VALUE two_mc) 
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  int status;
+  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
+  CHECK_FIXNUM(two_ma); CHECK_FIXNUM(two_mb); CHECK_FIXNUM(two_mc);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  status = gsl_sf_coupling_3j_e(FIX2INT(two_ja), FIX2INT(two_jb),
+			  FIX2INT(two_jc), FIX2INT(two_ma),
+			  FIX2INT(two_mb), FIX2INT(two_mc),
+			  rslt);
+
+  return v;
+}
+
+static VALUE rb_gsl_sf_coupling_6j(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_jd, VALUE two_je, VALUE two_jf) 
+{
+  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
+  CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
+  return rb_float_new(gsl_sf_coupling_6j(FIX2INT(two_ja), FIX2INT(two_jb),
+					 FIX2INT(two_jc), FIX2INT(two_jd),
+					 FIX2INT(two_je), FIX2INT(two_jf)));
+}
+
+static VALUE rb_gsl_sf_coupling_6j_e(VALUE obj, VALUE two_ja, VALUE two_jb, VALUE two_jc, VALUE two_jd, VALUE two_je, VALUE two_jf) 
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  int status;
+  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
+  CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  status = gsl_sf_coupling_6j_e(FIX2INT(two_ja), FIX2INT(two_jb),
+			  FIX2INT(two_jc), FIX2INT(two_jd),
+			  FIX2INT(two_je), FIX2INT(two_jf),
+			  rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_coupling_9j(VALUE obj, VALUE two_ja, VALUE two_jb,
+				   VALUE two_jc, VALUE two_jd, VALUE two_je,
+				   VALUE two_jf, VALUE two_jg, VALUE two_jh,
+				   VALUE two_ji) 
+{
+  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
+  CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
+  CHECK_FIXNUM(two_jg); CHECK_FIXNUM(two_jh); CHECK_FIXNUM(two_ji);
+  return rb_float_new(gsl_sf_coupling_9j(FIX2INT(two_ja), FIX2INT(two_jb),
+					FIX2INT(two_jc), FIX2INT(two_jd),
+					FIX2INT(two_je), FIX2INT(two_jf),
+					FIX2INT(two_jg), FIX2INT(two_jh),
+					FIX2INT(two_ji)));
+}
+
+static VALUE rb_gsl_sf_coupling_9j_e(VALUE obj, VALUE two_ja, VALUE two_jb,
+				   VALUE two_jc, VALUE two_jd, VALUE two_je,
+				   VALUE two_jf, VALUE two_jg, VALUE two_jh,
+				   VALUE two_ji) 
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  int status;
+  CHECK_FIXNUM(two_ja); CHECK_FIXNUM(two_jb); CHECK_FIXNUM(two_jc);
+  CHECK_FIXNUM(two_jd); CHECK_FIXNUM(two_je); CHECK_FIXNUM(two_jf);
+  CHECK_FIXNUM(two_jg); CHECK_FIXNUM(two_jh); CHECK_FIXNUM(two_ji);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  status = gsl_sf_coupling_9j_e(FIX2INT(two_ja), FIX2INT(two_jb),
+			   FIX2INT(two_jc), FIX2INT(two_jd),
+			   FIX2INT(two_je), FIX2INT(two_jf),
+			   FIX2INT(two_jg), FIX2INT(two_jh),
+			   FIX2INT(two_ji), rslt);
+  return v;
+}
+
+void Init_gsl_sf_coupling(VALUE module)
+{
+  VALUE mgsl_sf_coupling;
+
+  rb_define_module_function(module, "coupling_3j",  rb_gsl_sf_coupling_3j, 6);
+  rb_define_module_function(module, "coupling_3j_e",  rb_gsl_sf_coupling_3j_e, 6);
+  rb_define_module_function(module, "coupling_6j",  rb_gsl_sf_coupling_6j, 6);
+  rb_define_module_function(module, "coupling_6j_e",  rb_gsl_sf_coupling_6j_e, 6);
+  rb_define_module_function(module, "coupling_9j",  rb_gsl_sf_coupling_9j, 9);
+  rb_define_module_function(module, "coupling_9j_e",  rb_gsl_sf_coupling_9j_e, 9);
+
+  mgsl_sf_coupling = rb_define_module_under(module, "Coupling");
+  
+  rb_define_module_function(mgsl_sf_coupling, "3j",  rb_gsl_sf_coupling_3j, 6);
+  rb_define_module_function(mgsl_sf_coupling, "3j_e",  rb_gsl_sf_coupling_3j_e, 6);
+  rb_define_module_function(mgsl_sf_coupling, "6j",  rb_gsl_sf_coupling_6j, 6);
+  rb_define_module_function(mgsl_sf_coupling, "6j_e",  rb_gsl_sf_coupling_6j_e, 6);
+  rb_define_module_function(mgsl_sf_coupling, "9j",  rb_gsl_sf_coupling_9j, 9);
+  rb_define_module_function(mgsl_sf_coupling, "9j_e",  rb_gsl_sf_coupling_9j_e, 9);
+}

data/ext/sf_dawson.c

@@ -0,0 +1,29 @@
+/*
+  sf_dawson.c
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2006 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#include "rb_gsl_sf.h"
+
+static VALUE rb_gsl_sf_dawson(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_dawson, x);
+}
+
+static VALUE rb_gsl_sf_dawson_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_dawson_e, x);
+}
+
+void Init_gsl_sf_dawson(VALUE module)
+{
+  rb_define_module_function(module, "dawson",  rb_gsl_sf_dawson, 1);
+  rb_define_module_function(module, "dawson_e",  rb_gsl_sf_dawson_e, 1);
+}

data/ext/sf_debye.c

@@ -0,0 +1,157 @@
+/*
+  sf_debye.c
+  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
+    (C) Copyright 2001-2006 by Yoshiki Tsunesada
+
+  Ruby/GSL is free software: you can redistribute it and/or modify it
+  under the terms of the GNU General Public License.
+  This library is distributed in the hope that it will be useful, but
+  WITHOUT ANY WARRANTY; without even the implied warranty of
+  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+*/
+
+#include "rb_gsl_sf.h"
+
+static VALUE rb_gsl_sf_debye_1(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_debye_1, x);
+}
+
+static VALUE rb_gsl_sf_debye_1_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_debye_1_e, x);
+}
+
+static VALUE rb_gsl_sf_debye_2(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_debye_2, x);
+}
+
+static VALUE rb_gsl_sf_debye_2_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_debye_2_e, x);
+}
+
+static VALUE rb_gsl_sf_debye_3(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_debye_3, x);
+}
+
+static VALUE rb_gsl_sf_debye_3_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_debye_3_e, x);
+}
+
+static VALUE rb_gsl_sf_debye_4(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_debye_4, x);
+}
+
+static VALUE rb_gsl_sf_debye_4_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_debye_4_e, x);
+}
+
+static VALUE rb_gsl_sf_debye_n(int argc, VALUE *argv, VALUE obj)
+{
+  int n;
+  VALUE x;
+  switch (argc) {
+  case 1:
+    n = 1;
+    x = argv[0];
+    break;
+  case 2:
+    CHECK_FIXNUM(argv[0]);
+    n = FIX2INT(argv[0]);
+    x = argv[1];
+    break;
+  default:
+    rb_raise(rb_eArgError, "wrong number of arguments (%d for 1 or 2)", argc);
+    break;
+  }
+  switch (n) {
+  case 1:
+    return rb_gsl_sf_eval1(gsl_sf_debye_1, x);
+    break;
+  case 2:
+    return rb_gsl_sf_eval1(gsl_sf_debye_2, x);
+    break;
+  case 3:
+    return rb_gsl_sf_eval1(gsl_sf_debye_3, x);
+    break;
+  case 4:
+    return rb_gsl_sf_eval1(gsl_sf_debye_4, x);
+    break;
+#ifdef GSL_1_8_LATER
+  case 5:
+    return rb_gsl_sf_eval1(gsl_sf_debye_5, x);
+    break;
+  case 6:
+    return rb_gsl_sf_eval1(gsl_sf_debye_6, x);
+    break;
+#endif
+  default:
+    rb_raise(rb_eRuntimeError, "n must be 1, 2, 3, or 4");
+    break;
+  }
+}
+
+#ifdef GSL_1_8_LATER
+static VALUE rb_gsl_sf_debye_5(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_debye_5, x);
+}
+
+static VALUE rb_gsl_sf_debye_5_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_debye_5_e, x);
+}
+static VALUE rb_gsl_sf_debye_6(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_debye_6, x);
+}
+
+static VALUE rb_gsl_sf_debye_6_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_debye_6_e, x);
+}
+
+#endif
+
+void Init_gsl_sf_debye(VALUE module)
+{
+  VALUE mgsl_sf_debye;
+  rb_define_module_function(module, "debye_1",  rb_gsl_sf_debye_1, 1);
+  rb_define_module_function(module, "debye_1_e",  rb_gsl_sf_debye_1_e, 1);
+  rb_define_module_function(module, "debye_2",  rb_gsl_sf_debye_2, 1);
+  rb_define_module_function(module, "debye_2_e",  rb_gsl_sf_debye_2_e, 1);
+  rb_define_module_function(module, "debye_3",  rb_gsl_sf_debye_3, 1);
+  rb_define_module_function(module, "debye_3_e",  rb_gsl_sf_debye_3_e, 1);
+  rb_define_module_function(module, "debye_4",  rb_gsl_sf_debye_4, 1);
+  rb_define_module_function(module, "debye_4_e",  rb_gsl_sf_debye_4_e, 1);
+#ifdef GSL_1_8_LATER
+  rb_define_module_function(module, "debye_5",  rb_gsl_sf_debye_5, 1);
+  rb_define_module_function(module, "debye_5_e",  rb_gsl_sf_debye_5_e, 1);
+  rb_define_module_function(module, "debye_6",  rb_gsl_sf_debye_6, 1);
+  rb_define_module_function(module, "debye_6_e",  rb_gsl_sf_debye_6_e, 1);
+#endif
+  rb_define_module_function(module, "debye_n",  rb_gsl_sf_debye_n, -1);
+
+  mgsl_sf_debye = rb_define_module_under(module, "Debye");
+  rb_define_module_function(mgsl_sf_debye, "one",  rb_gsl_sf_debye_1, 1);
+  rb_define_module_function(mgsl_sf_debye, "one_e",  rb_gsl_sf_debye_1_e, 1);
+  rb_define_module_function(mgsl_sf_debye, "two",  rb_gsl_sf_debye_2, 1);
+  rb_define_module_function(mgsl_sf_debye, "two_e",  rb_gsl_sf_debye_2_e, 1);
+  rb_define_module_function(mgsl_sf_debye, "three",  rb_gsl_sf_debye_3, 1);
+  rb_define_module_function(mgsl_sf_debye, "three_e",  rb_gsl_sf_debye_3_e, 1);
+  rb_define_module_function(mgsl_sf_debye, "four",  rb_gsl_sf_debye_4, 1);
+  rb_define_module_function(mgsl_sf_debye, "four_e",  rb_gsl_sf_debye_4_e, 1);
+#ifdef GSL_1_8_LATER
+  rb_define_module_function(mgsl_sf_debye, "five",  rb_gsl_sf_debye_5, 1);
+  rb_define_module_function(mgsl_sf_debye, "five_e",  rb_gsl_sf_debye_5_e, 1);
+  rb_define_module_function(mgsl_sf_debye, "six",  rb_gsl_sf_debye_6, 1);
+  rb_define_module_function(mgsl_sf_debye, "six_e",  rb_gsl_sf_debye_6_e, 1);
+#endif
+  rb_define_module_function(mgsl_sf_debye, "n",  rb_gsl_sf_debye_n, -1);
+}