rb-gsl 1.16.0.5 → 1.16.0.6

data/ext/gsl_native/sf_clausen.c

@@ -1,28 +0,0 @@
-/*
-  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 "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/gsl_native/sf_coulomb.c

@@ -1,206 +0,0 @@
-/*
-  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 "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;
-  // local variable "status" declared and set, but never used
-  //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;
-  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);
-
-  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;
-  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);
-
-  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;
-  // local variable "status" declared and set, but never used
-  //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/gsl_native/sf_coupling.c

@@ -1,121 +0,0 @@
-/*
-  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 "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;
-  // local variable "status" declared and set, but never used
-  //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;
-  // local variable "status" declared and set, but never used
-  //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;
-  // local variable "status" declared and set, but never used
-  //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/gsl_native/sf_dawson.c

@@ -1,29 +0,0 @@
-/*
-  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 "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/gsl_native/sf_debye.c

@@ -1,148 +0,0 @@
-/*
-  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 "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;
-  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;
-  default:
-    rb_raise(rb_eRuntimeError, "n must be 1, 2, 3, or 4");
-    break;
-  }
-}
-
-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);
-}
-
-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);
-  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);
-  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);
-  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);
-  rb_define_module_function(mgsl_sf_debye, "n",  rb_gsl_sf_debye_n, -1);
-}