rb-gsl 1.15.3.1

data/ext/sf_gegenbauer.c

@@ -0,0 +1,97 @@
+/*
+  sf_gegenbauer.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_gegenpoly_1(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_gegenpoly_1, lambda, x);
+}
+
+static VALUE rb_gsl_sf_gegenpoly_1_e(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_gegenpoly_1_e, lambda, x);
+}
+
+static VALUE rb_gsl_sf_gegenpoly_2(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_gegenpoly_2, lambda, x);
+}
+
+static VALUE rb_gsl_sf_gegenpoly_2_e(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_gegenpoly_2_e, lambda, x);
+}
+
+static VALUE rb_gsl_sf_gegenpoly_3(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_gegenpoly_3, lambda, x);
+}
+
+static VALUE rb_gsl_sf_gegenpoly_3_e(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_gegenpoly_3_e, lambda, x);
+}
+
+static VALUE rb_gsl_sf_gegenpoly_n(VALUE obj, VALUE n, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_int_double_double(gsl_sf_gegenpoly_n, n, lambda, x);
+}
+
+static VALUE rb_gsl_sf_gegenpoly_n_e(VALUE obj, VALUE n, VALUE lambda, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  CHECK_FIXNUM(n);
+  Need_Float(lambda); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  /*status =*/ gsl_sf_gegenpoly_n_e(FIX2INT(n), NUM2DBL(lambda), NUM2DBL(x), rslt);
+
+  return v;
+}
+
+static VALUE rb_gsl_sf_gegenpoly_array(VALUE obj, VALUE nmax, VALUE lambda, VALUE x)
+{
+  gsl_vector *v = NULL;
+  CHECK_FIXNUM(nmax);
+  Need_Float(lambda); Need_Float(x);
+  v = gsl_vector_alloc(nmax);
+  gsl_sf_gegenpoly_array(FIX2INT(nmax), NUM2DBL(lambda), NUM2DBL(x), v->data);
+  return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
+}
+
+void Init_gsl_sf_gegenbauer(VALUE module)
+{
+  VALUE mgsl_sf_gegenpoly;
+  rb_define_module_function(module, "gegenpoly_1",  rb_gsl_sf_gegenpoly_1, 2);
+  rb_define_module_function(module, "gegenpoly_1_e",  rb_gsl_sf_gegenpoly_1_e, 2);
+  rb_define_module_function(module, "gegenpoly_2",  rb_gsl_sf_gegenpoly_2, 2);
+  rb_define_module_function(module, "gegenpoly_2_e",  rb_gsl_sf_gegenpoly_2_e, 2);
+  rb_define_module_function(module, "gegenpoly_3",  rb_gsl_sf_gegenpoly_3, 2);
+  rb_define_module_function(module, "gegenpoly_3_e",  rb_gsl_sf_gegenpoly_3_e, 2);
+  rb_define_module_function(module, "gegenpoly_n",  rb_gsl_sf_gegenpoly_n, 3);
+  rb_define_module_function(module, "gegenpoly_n_e",  rb_gsl_sf_gegenpoly_n_e, 3);
+  rb_define_module_function(module, "gegenpoly_array",  rb_gsl_sf_gegenpoly_array, 3);
+
+  mgsl_sf_gegenpoly = rb_define_module_under(module, "Gegenpoly");
+  rb_define_module_function(mgsl_sf_gegenpoly, "one",  rb_gsl_sf_gegenpoly_1, 2);
+  rb_define_module_function(mgsl_sf_gegenpoly, "one_e",  rb_gsl_sf_gegenpoly_1_e, 2);
+  rb_define_module_function(mgsl_sf_gegenpoly, "two",  rb_gsl_sf_gegenpoly_2, 2);
+  rb_define_module_function(mgsl_sf_gegenpoly, "two_e",  rb_gsl_sf_gegenpoly_2_e, 2);
+  rb_define_module_function(mgsl_sf_gegenpoly, "three",  rb_gsl_sf_gegenpoly_3, 2);
+  rb_define_module_function(mgsl_sf_gegenpoly, "three_e",  rb_gsl_sf_gegenpoly_3_e, 2);
+  rb_define_module_function(mgsl_sf_gegenpoly, "n",  rb_gsl_sf_gegenpoly_n, 3);
+  rb_define_module_function(mgsl_sf_gegenpoly, "n_e",  rb_gsl_sf_gegenpoly_n_e, 3);
+  rb_define_module_function(mgsl_sf_gegenpoly, "array",  rb_gsl_sf_gegenpoly_array, 3);
+
+}

data/ext/sf_hyperg.c

@@ -0,0 +1,203 @@
+/*
+  sf_hyperg.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"
+
+/* Checked with Mathematica. */
+static VALUE rb_gsl_sf_hyperg_0F1(VALUE obj, VALUE c, VALUE x)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_hyperg_0F1, c, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_0F1_e(VALUE obj, VALUE c, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_hyperg_0F1_e, c, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_1F1_int(VALUE obj, VALUE m, VALUE n, VALUE x)
+{
+  return rb_gsl_sf_eval_int_int_double(gsl_sf_hyperg_1F1_int, m, n, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_1F1_int_e(VALUE obj, VALUE m, VALUE n, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  CHECK_FIXNUM(m); CHECK_FIXNUM(n); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  /*status =*/ gsl_sf_hyperg_1F1_int_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_hyperg_1F1(VALUE obj, VALUE a, VALUE b, VALUE x)
+{
+  return rb_gsl_sf_eval_double3(gsl_sf_hyperg_1F1, a, b, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_1F1_e(VALUE obj, VALUE a, VALUE b, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double3(gsl_sf_hyperg_1F1_e, a, b, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_U_int(VALUE obj, VALUE m, VALUE n, VALUE x)
+{
+  return rb_gsl_sf_eval_int_int_double(gsl_sf_hyperg_U_int, m, n, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_U_int_e(VALUE obj, VALUE m, VALUE n, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  CHECK_FIXNUM(m); CHECK_FIXNUM(n); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  /*status =*/ gsl_sf_hyperg_U_int_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_hyperg_U_int_e10_e(VALUE obj, VALUE m, VALUE n, VALUE x)
+{
+  gsl_sf_result_e10 *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  CHECK_FIXNUM(m); CHECK_FIXNUM(n); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
+  /*status =*/ gsl_sf_hyperg_U_int_e10_e(FIX2INT(m), FIX2INT(n), NUM2DBL(x), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_hyperg_U(VALUE obj, VALUE a, VALUE b, VALUE x)
+{
+  return rb_gsl_sf_eval_double3(gsl_sf_hyperg_U, a, b, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_U_e(VALUE obj, VALUE a, VALUE b, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double3(gsl_sf_hyperg_U_e, a, b, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_U_e10_e(VALUE obj, VALUE a, VALUE b, VALUE x)
+{
+  gsl_sf_result_e10 *rslt = NULL;
+  VALUE v;
+  int status;
+  Need_Float(a); Need_Float(b); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
+  status = gsl_sf_hyperg_U_e10_e(NUM2DBL(a), NUM2DBL(b), NUM2DBL(x), rslt);
+  return rb_ary_new3(2, v, INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_hyperg_2F1(VALUE obj, VALUE a, VALUE b, VALUE c, VALUE x)
+{
+  return rb_gsl_sf_eval_double4(gsl_sf_hyperg_2F1, a, b, c, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_2F1_e(VALUE obj, VALUE a, VALUE b, VALUE c, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  int status;
+  Need_Float(a); Need_Float(b); Need_Float(c); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  status = gsl_sf_hyperg_2F1_e(NUM2DBL(a), NUM2DBL(b), NUM2DBL(c), NUM2DBL(x), rslt);
+  return rb_ary_new3(2, v, INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_hyperg_2F1_conj(VALUE obj, VALUE aR, VALUE aI, VALUE c, VALUE x)
+{
+  return rb_gsl_sf_eval_double4(gsl_sf_hyperg_2F1_conj, aR, aI, c, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_2F1_conj_e(VALUE obj, VALUE aR, VALUE aI, VALUE c, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  Need_Float(aR); Need_Float(aI);  Need_Float(c); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  /*status =*/ gsl_sf_hyperg_2F1_conj_e(NUM2DBL(aR), NUM2DBL(aI), NUM2DBL(c), NUM2DBL(x), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_hyperg_2F1_renorm(VALUE obj, VALUE a, VALUE b, VALUE c, VALUE x)
+{
+  return rb_gsl_sf_eval_double4(gsl_sf_hyperg_2F1_renorm, a, b, c, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_2F1_renorm_e(VALUE obj, VALUE a, VALUE b, VALUE c, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  Need_Float(a); Need_Float(b);  Need_Float(c); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  /*status =*/ gsl_sf_hyperg_2F1_renorm_e(NUM2DBL(a), NUM2DBL(b), NUM2DBL(c), NUM2DBL(x), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_hyperg_2F1_conj_renorm(VALUE obj, VALUE aR, VALUE aI, VALUE c, VALUE x)
+{
+  return rb_gsl_sf_eval_double4(gsl_sf_hyperg_2F1_conj_renorm, aR, aI, c, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_2F1_conj_renorm_e(VALUE obj, VALUE aR, VALUE aI, VALUE c, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  Need_Float(aR); Need_Float(aI);  Need_Float(c); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  /*status =*/ gsl_sf_hyperg_2F1_conj_renorm_e(NUM2DBL(aR), NUM2DBL(aI), NUM2DBL(c), NUM2DBL(x), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_hyperg_2F0(VALUE obj, VALUE a, VALUE b,  VALUE x)
+{
+  return rb_gsl_sf_eval_double3(gsl_sf_hyperg_2F0, a, b, x);
+}
+
+static VALUE rb_gsl_sf_hyperg_2F0_e(VALUE obj, VALUE a, VALUE b, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double3(gsl_sf_hyperg_2F0_e, a, b, x);
+}
+
+void Init_gsl_sf_hyperg(VALUE module)
+{
+  rb_define_module_function(module, "hyperg_0F1",  rb_gsl_sf_hyperg_0F1, 2);
+  rb_define_module_function(module, "hyperg_0F1_e",  rb_gsl_sf_hyperg_0F1_e, 2);
+  rb_define_module_function(module, "hyperg_1F1_int",  rb_gsl_sf_hyperg_1F1_int, 3);
+  rb_define_module_function(module, "hyperg_1F1_int_e",  rb_gsl_sf_hyperg_1F1_int_e, 3);
+  rb_define_module_function(module, "hyperg_1F1",  rb_gsl_sf_hyperg_1F1, 3);
+  rb_define_module_function(module, "hyperg_1F1_e",  rb_gsl_sf_hyperg_1F1_e, 3);
+  rb_define_module_function(module, "hyperg_U_int",  rb_gsl_sf_hyperg_U_int, 3);
+  rb_define_module_function(module, "hyperg_U_int_e",  rb_gsl_sf_hyperg_U_int_e, 3);
+  rb_define_module_function(module, "hyperg_U_int_e10_e",  rb_gsl_sf_hyperg_U_int_e10_e, 3);
+  rb_define_module_function(module, "hyperg_U",  rb_gsl_sf_hyperg_U, 3);
+  rb_define_module_function(module, "hyperg_U_e",  rb_gsl_sf_hyperg_U_e, 3);
+  rb_define_module_function(module, "hyperg_U_e10_e",  rb_gsl_sf_hyperg_U_e10_e, 3);
+  rb_define_module_function(module, "hyperg_2F1",  rb_gsl_sf_hyperg_2F1, 4);
+  rb_define_module_function(module, "hyperg_2F1_e",  rb_gsl_sf_hyperg_2F1_e, 4);
+  rb_define_module_function(module, "hyperg_2F1_conj",  rb_gsl_sf_hyperg_2F1_conj, 4);
+  rb_define_module_function(module, "hyperg_2F1_conj_e",  rb_gsl_sf_hyperg_2F1_conj_e, 4);
+  rb_define_module_function(module, "hyperg_2F1_renorm",  rb_gsl_sf_hyperg_2F1_renorm, 4);
+  rb_define_module_function(module, "hyperg_2F1_renorm_e",  rb_gsl_sf_hyperg_2F1_renorm_e, 4);
+  rb_define_module_function(module, "hyperg_2F1_conj_renorm",  rb_gsl_sf_hyperg_2F1_conj_renorm, 4);
+  rb_define_module_function(module, "hyperg_2F1_conj_renorm_e",  rb_gsl_sf_hyperg_2F1_conj_renorm_e, 4);
+  rb_define_module_function(module, "hyperg_2F0",  rb_gsl_sf_hyperg_2F0, 3);
+  rb_define_module_function(module, "hyperg_2F0_e",  rb_gsl_sf_hyperg_2F0_e, 3);
+}

data/ext/sf_laguerre.c

@@ -0,0 +1,113 @@
+/*
+  sf_laguerre.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_laguerre_X(int argc, VALUE *argv, VALUE obj,
+				  double (*f)(double, double))
+{
+  switch (argc) {
+  case 2:
+    return rb_gsl_sf_eval_double_double(f, argv[0], argv[1]);
+    break;
+  case 1:
+    return rb_gsl_sf_eval_double_double(f, INT2FIX(0), argv[0]);
+    break;
+  default:
+    rb_raise(rb_eArgError, "wrong number of arguments (%d for 1 or 2)", argc);
+    break;
+  }
+}
+
+static VALUE rb_gsl_sf_laguerre_1(int argc, VALUE *argv, VALUE obj)
+{
+  return rb_gsl_sf_laguerre_X(argc, argv, obj, gsl_sf_laguerre_1);
+}
+
+static VALUE rb_gsl_sf_laguerre_1_e(VALUE obj, VALUE a, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_laguerre_1_e, a, x);
+}
+
+static VALUE rb_gsl_sf_laguerre_2(int argc, VALUE *argv, VALUE obj)
+{
+  return rb_gsl_sf_laguerre_X(argc, argv, obj, gsl_sf_laguerre_2);
+}
+
+static VALUE rb_gsl_sf_laguerre_2_e(VALUE obj, VALUE a, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_laguerre_2_e, a, x);
+}
+
+static VALUE rb_gsl_sf_laguerre_3(int argc, VALUE *argv, VALUE obj)
+{
+  return rb_gsl_sf_laguerre_X(argc, argv, obj, gsl_sf_laguerre_3);
+}
+
+static VALUE rb_gsl_sf_laguerre_3_e(VALUE obj, VALUE a, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_laguerre_3_e, a, x);
+}
+
+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]);
+    break;
+  case 2:
+    return rb_gsl_sf_eval_int_double_double(gsl_sf_laguerre_n, argv[0],
+					    INT2FIX(0), argv[1]);
+    break;
+  default:
+    rb_raise(rb_eArgError, "wrong number of arguments (%d for 1 or 2)", argc);
+    break;
+  }
+}
+
+static VALUE rb_gsl_sf_laguerre_n_e(VALUE obj, VALUE n, VALUE a, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  CHECK_FIXNUM(n);
+  Need_Float(a); Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  /*status =*/ gsl_sf_laguerre_n_e(FIX2INT(n), NUM2DBL(a), NUM2DBL(x), rslt);
+  return v;
+}
+
+void Init_gsl_sf_laguerre(VALUE module)
+{
+  VALUE mgsl_sf_laguerre;
+
+  rb_define_module_function(module, "laguerre_1",  rb_gsl_sf_laguerre_1, -1);
+  rb_define_module_function(module, "laguerre_1_e",  rb_gsl_sf_laguerre_1_e, 2);
+  rb_define_module_function(module, "laguerre_2",  rb_gsl_sf_laguerre_2, -1);
+  rb_define_module_function(module, "laguerre_2_e",  rb_gsl_sf_laguerre_2_e, 2);
+  rb_define_module_function(module, "laguerre_3",  rb_gsl_sf_laguerre_3, -1);
+  rb_define_module_function(module, "laguerre_3_e",  rb_gsl_sf_laguerre_3_e, 2);
+  rb_define_module_function(module, "laguerre_n",  rb_gsl_sf_laguerre_n, -1);
+  rb_define_module_function(module, "laguerre_n_e",  rb_gsl_sf_laguerre_n_e, 3);
+
+  mgsl_sf_laguerre = rb_define_module_under(module, "Laguerre");
+  rb_define_module_function(mgsl_sf_laguerre, "one",  rb_gsl_sf_laguerre_1, -1);
+  rb_define_module_function(mgsl_sf_laguerre, "one_e",  rb_gsl_sf_laguerre_1_e, 2);
+  rb_define_module_function(mgsl_sf_laguerre, "two",  rb_gsl_sf_laguerre_2, -1);
+  rb_define_module_function(mgsl_sf_laguerre, "two_e",  rb_gsl_sf_laguerre_2_e, 2);
+  rb_define_module_function(mgsl_sf_laguerre, "three_3",  rb_gsl_sf_laguerre_3, -1);
+  rb_define_module_function(mgsl_sf_laguerre, "three_e",  rb_gsl_sf_laguerre_3_e, 2);
+  rb_define_module_function(mgsl_sf_laguerre, "n",  rb_gsl_sf_laguerre_n, -1);
+  rb_define_module_function(mgsl_sf_laguerre, "n_e",  rb_gsl_sf_laguerre_n_e, 3);
+
+}

data/ext/sf_lambert.c

@@ -0,0 +1,47 @@
+/*
+  sf_lambert.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_lambert_W0(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_lambert_W0, x);
+}
+
+static VALUE rb_gsl_sf_lambert_W0_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_lambert_W0_e, x);
+}
+
+static VALUE rb_gsl_sf_lambert_Wm1(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_lambert_Wm1, x);
+}
+
+static VALUE rb_gsl_sf_lambert_Wm1_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_lambert_Wm1_e, x);
+}
+
+void Init_gsl_sf_lambert(VALUE module)
+{
+  VALUE mgsl_sf_lambert;
+  rb_define_module_function(module, "lambert_W0",  rb_gsl_sf_lambert_W0, 1);
+  rb_define_module_function(module, "lambert_W0_e",  rb_gsl_sf_lambert_W0_e, 1);
+  rb_define_module_function(module, "lambert_Wm1",  rb_gsl_sf_lambert_Wm1, 1);
+  rb_define_module_function(module, "lambert_Wm1_e",  rb_gsl_sf_lambert_Wm1_e, 1);
+
+  mgsl_sf_lambert = rb_define_module_under(module, "Lambert");
+  rb_define_module_function(mgsl_sf_lambert, "W0",  rb_gsl_sf_lambert_W0, 1);
+  rb_define_module_function(mgsl_sf_lambert, "W0_e",  rb_gsl_sf_lambert_W0_e, 1);
+  rb_define_module_function(mgsl_sf_lambert, "Wm1",  rb_gsl_sf_lambert_Wm1, 1);
+  rb_define_module_function(mgsl_sf_lambert, "Wm1_e",  rb_gsl_sf_lambert_Wm1_e, 1);
+}