gsl 1.12.108

data/ext/sf_dilog.c

@@ -0,0 +1,42 @@
+/*
+  sf_dilog.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_dilog(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_dilog, x);
+}
+
+static VALUE rb_gsl_sf_dilog_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_dilog_e, x);
+}
+
+static VALUE rb_gsl_sf_complex_dilog_e(VALUE obj, VALUE r, VALUE theta)
+{
+  gsl_sf_result *re, *im;
+  VALUE vre, vim;
+  int status;
+  Need_Float(r); Need_Float(theta);
+  vre = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, re);
+  vim = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, im);
+  status = gsl_sf_complex_dilog_e(NUM2DBL(r), NUM2DBL(theta), re, im);
+  return rb_ary_new3(2, vre, vim);
+}
+
+void Init_gsl_sf_dilog(VALUE module)
+{
+  rb_define_module_function(module, "dilog",  rb_gsl_sf_dilog, 1);
+  rb_define_module_function(module, "dilog_e",  rb_gsl_sf_dilog_e, 1);
+  rb_define_module_function(module, "complex_dilog_e",  rb_gsl_sf_complex_dilog_e, 2);
+}

data/ext/sf_elementary.c

@@ -0,0 +1,44 @@
+/*
+  sf_elementary.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_multiply_e(VALUE obj, VALUE x, VALUE y)
+{
+  gsl_sf_result *r;
+  VALUE v;
+  int status;
+  Need_Float(x); Need_Float(y);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r);
+  status = gsl_sf_multiply_e(NUM2DBL(x), NUM2DBL(y), r);
+  return v;
+}
+
+static VALUE rb_gsl_sf_multiply_err_e(VALUE obj, VALUE x, VALUE dx,
+				      VALUE y, VALUE dy)
+{
+  gsl_sf_result *r;
+  VALUE v;
+  int status;
+  Need_Float(x); Need_Float(y);
+  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);
+  return v;
+}
+
+void Init_gsl_sf_elementary(VALUE module)
+{
+  rb_define_module_function(module, "multiply_e",  rb_gsl_sf_multiply_e, 2);
+  rb_define_module_function(module, "multiply_err_e",  rb_gsl_sf_multiply_err_e, 4);
+}

data/ext/sf_ellint.c

@@ -0,0 +1,206 @@
+/*
+  sf_ellint.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_ellint_Kcomp(int argc, VALUE *argv, VALUE obj)
+{
+  if (argc == 1) return eval_sf(gsl_sf_ellint_Kcomp, argv[0]);
+  else return rb_gsl_sf_eval_double_m(gsl_sf_ellint_Kcomp, argv[0], argv[1]);
+}
+
+static VALUE rb_gsl_sf_ellint_Kcomp_e(VALUE obj, VALUE k, VALUE m)
+{
+  return rb_gsl_sf_eval_e_m(gsl_sf_ellint_Kcomp_e, k, m);
+}
+
+static VALUE rb_gsl_sf_ellint_Ecomp(int argc, VALUE *argv, VALUE obj)
+{
+  if (argc == 1) return eval_sf(gsl_sf_ellint_Ecomp, argv[0]);
+  else return rb_gsl_sf_eval_double_m(gsl_sf_ellint_Ecomp, argv[0], argv[1]);
+}
+
+static VALUE rb_gsl_sf_ellint_Ecomp_e(VALUE obj, VALUE k, VALUE m)
+{
+  return rb_gsl_sf_eval_e_m(gsl_sf_ellint_Ecomp_e, k, 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 
+    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_F, argv[0], argv[1], argv[2]);
+}
+
+static VALUE rb_gsl_sf_ellint_F_e(VALUE obj, VALUE phi, VALUE k, VALUE m)
+{
+  return rb_gsl_sf_eval_e_double2_m(gsl_sf_ellint_F_e, phi, k, 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 
+    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_E, argv[0], argv[1], argv[2]);
+}
+
+static VALUE rb_gsl_sf_ellint_E_e(VALUE obj, VALUE phi, VALUE k, VALUE m)
+{
+  return rb_gsl_sf_eval_e_double2_m(gsl_sf_ellint_E_e, phi, k, m);
+}
+
+static VALUE rb_gsl_sf_ellint_P(int argc, VALUE *argv, VALUE obj)
+
+{
+  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 
+    return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_P, argv[0], argv[1], argv[2],
+				    argv[3]);
+}
+
+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) 
+    return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_D, argv[0], argv[1], argv[2],
+				    INT2FIX(GSL_PREC_DOUBLE));
+  else 
+    return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_D, argv[0], argv[1], argv[2],
+				    argv[3]);
+}
+
+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 
+    return rb_gsl_sf_eval_double2_m(gsl_sf_ellint_RC, argv[0], argv[1], argv[2]);
+}
+
+static VALUE rb_gsl_sf_ellint_RC_e(VALUE obj, VALUE x, VALUE y, VALUE m)
+{
+  return rb_gsl_sf_eval_e_double2_m(gsl_sf_ellint_RC_e, x, y, m);
+}
+
+static VALUE rb_gsl_sf_ellint_RD(int argc, VALUE *argv, VALUE obj)
+{
+  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 
+    return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_RD, argv[0], argv[1], argv[2],
+				    argv[3]);
+}
+
+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) 
+    return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_RF, argv[0], argv[1], argv[2],
+				    INT2FIX(GSL_PREC_DOUBLE));
+  else 
+    return rb_gsl_sf_eval_double3_m(gsl_sf_ellint_RF, argv[0], argv[1], argv[2],
+				    argv[3]);
+}
+
+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) 
+    return rb_gsl_sf_eval_double4_m(gsl_sf_ellint_RJ, argv[0], argv[1], argv[2],
+				    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]);
+}
+
+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);
+}
+
+void Init_gsl_sf_ellint(VALUE module)
+{
+  VALUE mgsl_sf_ellint;
+
+  rb_define_module_function(module, "ellint_Kcomp",  rb_gsl_sf_ellint_Kcomp, -1);
+  rb_define_module_function(module, "ellint_Kcomp_e",  rb_gsl_sf_ellint_Kcomp_e, 2);
+  rb_define_module_function(module, "ellint_Ecomp",  rb_gsl_sf_ellint_Ecomp, -1);
+  rb_define_module_function(module, "ellint_Ecomp_e",  rb_gsl_sf_ellint_Ecomp_e, 2);
+  rb_define_module_function(module, "ellint_F",  rb_gsl_sf_ellint_F, -1);
+  rb_define_module_function(module, "ellint_F_e",  rb_gsl_sf_ellint_F_e, 3);
+  rb_define_module_function(module, "ellint_E",  rb_gsl_sf_ellint_E, -1);
+  rb_define_module_function(module, "ellint_E_e",  rb_gsl_sf_ellint_E_e, 3);
+  rb_define_module_function(module, "ellint_P",  rb_gsl_sf_ellint_P, -1);
+  rb_define_module_function(module, "ellint_P_e",  rb_gsl_sf_ellint_P_e, 4);
+  rb_define_module_function(module, "ellint_D",  rb_gsl_sf_ellint_D, -1);
+  rb_define_module_function(module, "ellint_D_e",  rb_gsl_sf_ellint_D_e, 4);
+  rb_define_module_function(module, "ellint_RC",  rb_gsl_sf_ellint_RC, -1);
+  rb_define_module_function(module, "ellint_RC_e",  rb_gsl_sf_ellint_RC_e, 3);
+  rb_define_module_function(module, "ellint_RD",  rb_gsl_sf_ellint_RD, -1);
+  rb_define_module_function(module, "ellint_RD_e",  rb_gsl_sf_ellint_RD_e, 4);
+  rb_define_module_function(module, "ellint_RF",  rb_gsl_sf_ellint_RF, -1);
+  rb_define_module_function(module, "ellint_RF_e",  rb_gsl_sf_ellint_RF_e, 4);
+  rb_define_module_function(module, "ellint_RJ",  rb_gsl_sf_ellint_RJ, -1);
+  rb_define_module_function(module, "ellint_RJ_e",  rb_gsl_sf_ellint_RJ_e, 5);
+
+  mgsl_sf_ellint = rb_define_module_under(module, "Ellint");
+  rb_define_module_function(mgsl_sf_ellint, "Kcomp",  rb_gsl_sf_ellint_Kcomp, -1);
+  rb_define_module_function(mgsl_sf_ellint, "Kcomp_e",  rb_gsl_sf_ellint_Kcomp_e, 2);
+  rb_define_module_function(mgsl_sf_ellint, "Ecomp",  rb_gsl_sf_ellint_Ecomp, -1);
+  rb_define_module_function(mgsl_sf_ellint, "Ecomp_e",  rb_gsl_sf_ellint_Ecomp_e, 2);
+  rb_define_module_function(mgsl_sf_ellint, "F",  rb_gsl_sf_ellint_F, -1);
+  rb_define_module_function(mgsl_sf_ellint, "F_e",  rb_gsl_sf_ellint_F_e, 3);
+  rb_define_module_function(mgsl_sf_ellint, "E",  rb_gsl_sf_ellint_E, -1);
+  rb_define_module_function(mgsl_sf_ellint, "E_e",  rb_gsl_sf_ellint_E_e, 3);
+  rb_define_module_function(mgsl_sf_ellint, "P",  rb_gsl_sf_ellint_P, -1);
+  rb_define_module_function(mgsl_sf_ellint, "P_e",  rb_gsl_sf_ellint_P_e, 4);
+  rb_define_module_function(mgsl_sf_ellint, "D",  rb_gsl_sf_ellint_D, -1);
+  rb_define_module_function(mgsl_sf_ellint, "D_e",  rb_gsl_sf_ellint_D_e, 4);
+  rb_define_module_function(mgsl_sf_ellint, "RC",  rb_gsl_sf_ellint_RC, -1);
+  rb_define_module_function(mgsl_sf_ellint, "RC_e",  rb_gsl_sf_ellint_RC_e, 3);
+  rb_define_module_function(mgsl_sf_ellint, "RD",  rb_gsl_sf_ellint_RD, -1);
+  rb_define_module_function(mgsl_sf_ellint, "RD_e",  rb_gsl_sf_ellint_RD_e, 4);
+  rb_define_module_function(mgsl_sf_ellint, "RF",  rb_gsl_sf_ellint_RF, -1);
+  rb_define_module_function(mgsl_sf_ellint, "RF_e",  rb_gsl_sf_ellint_RF_e, 4);
+  rb_define_module_function(mgsl_sf_ellint, "RJ",  rb_gsl_sf_ellint_RJ, -1);
+  rb_define_module_function(mgsl_sf_ellint, "RJ_e",  rb_gsl_sf_ellint_RJ_e, 5);
+}

data/ext/sf_elljac.c

@@ -0,0 +1,29 @@
+/*
+  sf_elljac.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_elljac_e(VALUE obj, VALUE n, VALUE m)
+{
+  double sn, cn, dn;
+  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));
+}
+
+void Init_gsl_sf_elljac(VALUE module)
+{
+  rb_define_module_function(module, "elljac_e",  rb_gsl_sf_elljac_e, 2);
+  rb_define_module_function(module, "elljac",  rb_gsl_sf_elljac_e, 2);
+}

data/ext/sf_erfc.c

@@ -0,0 +1,93 @@
+/*
+  sf_erfc.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_erf(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_erf, x);
+}
+
+static VALUE rb_gsl_sf_erf_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_erf_e, x);
+}
+
+static VALUE rb_gsl_sf_erfc(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_erfc, x);
+}
+
+static VALUE rb_gsl_sf_erfc_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_erfc_e, x);
+}
+
+static VALUE rb_gsl_sf_log_erfc(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_log_erfc, x);
+}
+
+static VALUE rb_gsl_sf_log_erfc_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_log_erfc_e, x);
+}
+
+static VALUE rb_gsl_sf_erf_Z(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_erf_Z, x);
+}
+
+static VALUE rb_gsl_sf_erf_Z_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_erf_Z_e, x);
+}
+
+static VALUE rb_gsl_sf_erf_Q(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_erf_Q, x);
+}
+
+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);
+}
+
+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)
+{
+  rb_define_module_function(module, "erf",  rb_gsl_sf_erf, 1);
+  rb_define_module_function(module, "erf_e",  rb_gsl_sf_erf_e, 1);
+  rb_define_module_function(module, "erfc",  rb_gsl_sf_erfc, 1);
+  rb_define_module_function(module, "erfc_e",  rb_gsl_sf_erfc_e, 1);
+  rb_define_module_function(module, "log_erfc",  rb_gsl_sf_log_erfc, 1);
+  rb_define_module_function(module, "log_erfc_e",  rb_gsl_sf_log_erfc_e, 1);
+  rb_define_module_function(module, "erf_Z",  rb_gsl_sf_erf_Z, 1);
+  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

@@ -0,0 +1,164 @@
+/*
+  sf_exp.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_exp(VALUE obj, VALUE x)
+{
+  if (COMPLEX_P(x) || VECTOR_COMPLEX_P(x) || MATRIX_COMPLEX_P(x)) {
+    return rb_gsl_math_complex_eval(gsl_complex_exp, x);
+  } else {
+    return rb_gsl_sf_eval1(gsl_sf_exp, x);
+  }
+}
+
+static VALUE rb_gsl_sf_exp_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_exp_e, x);
+}
+
+static VALUE rb_gsl_sf_exp_e10_e(VALUE obj, VALUE x)
+{
+  gsl_sf_result_e10 *rslt = NULL;
+  VALUE v;
+  int status;
+  Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
+  status = gsl_sf_exp_e10_e(NUM2DBL(x), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_exp_mult(VALUE obj, VALUE x, VALUE y)
+{
+  Need_Float(x);   Need_Float(y);
+  return rb_float_new(gsl_sf_exp_mult(NUM2DBL(x), NUM2DBL(y)));
+}
+
+static VALUE rb_gsl_sf_exp_mult_e(VALUE obj, VALUE x, VALUE y)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_exp_mult_e, x, y);
+}
+
+static VALUE rb_gsl_sf_exp_mult_e10_e(VALUE obj, VALUE x, VALUE y)
+{
+  gsl_sf_result_e10 *rslt = NULL;
+  VALUE v;
+  int status;
+  Need_Float(x);   Need_Float(y);
+  v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
+  status = gsl_sf_exp_mult_e10_e(NUM2DBL(x), NUM2DBL(y), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_expm1(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_expm1, x);
+}
+
+static VALUE rb_gsl_sf_expm1_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_expm1_e, x);
+}
+
+static VALUE rb_gsl_sf_exprel(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_exprel, x);
+}
+
+static VALUE rb_gsl_sf_exprel_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_exprel_e, x);
+}
+
+static VALUE rb_gsl_sf_exprel_2(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_exprel_2, x);
+}
+
+static VALUE rb_gsl_sf_exprel_2_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_exprel_2_e, x);
+}
+
+static VALUE rb_gsl_sf_exprel_n(VALUE obj, VALUE n, VALUE x)
+{
+  return rb_gsl_sf_eval_int_double(gsl_sf_exprel_n, n, x);
+}
+
+static VALUE rb_gsl_sf_exprel_n_e(VALUE obj, VALUE n, VALUE x)
+{
+  return rb_gsl_sf_eval_e_int_double(gsl_sf_exprel_n_e, n, x);
+}
+
+static VALUE rb_gsl_sf_exp_err_e(VALUE obj, VALUE x, VALUE dx)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_exp_err_e, x, dx);
+}
+
+static VALUE rb_gsl_sf_exp_err_e10_e(VALUE obj, VALUE x, VALUE dx)
+{
+  gsl_sf_result_e10 *rslt = NULL;
+  VALUE v;
+  int status;
+  Need_Float(x);   Need_Float(dx);
+  v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
+  status = gsl_sf_exp_err_e10_e(NUM2DBL(x), NUM2DBL(dx), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_exp_mult_err_e(VALUE obj, VALUE x, VALUE dx,
+				      VALUE y, VALUE dy)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  int status;
+  Need_Float(x);   Need_Float(y);
+  Need_Float(dx);   Need_Float(dy);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  status = gsl_sf_exp_mult_err_e(NUM2DBL(x), NUM2DBL(dx), NUM2DBL(y), NUM2DBL(dy), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_exp_mult_err_e10_e(VALUE obj, VALUE x, VALUE dx,
+					  VALUE y, VALUE dy)
+{
+  gsl_sf_result_e10 *rslt = NULL;
+  VALUE v;
+  int status;
+  Need_Float(x);   Need_Float(y);
+  Need_Float(dx);   Need_Float(dy);
+  v = Data_Make_Struct(cgsl_sf_result_e10, gsl_sf_result_e10, 0, free, rslt);
+  status = gsl_sf_exp_mult_err_e10_e(NUM2DBL(x), NUM2DBL(dx), NUM2DBL(y), NUM2DBL(dy), rslt);
+  return v;
+}
+
+void Init_gsl_sf_exp(VALUE module)
+{
+  rb_define_module_function(module, "exp",  rb_gsl_sf_exp, 1);
+  rb_define_module_function(module, "exp_e",  rb_gsl_sf_exp_e, 1);
+  rb_define_module_function(module, "exp_e10_e",  rb_gsl_sf_exp_e10_e, 1);
+  rb_define_module_function(module, "exp_mult",  rb_gsl_sf_exp_mult, 2);
+  rb_define_module_function(module, "exp_mult_e",  rb_gsl_sf_exp_mult_e, 2);
+  rb_define_module_function(module, "exp_mult_e10_e",  rb_gsl_sf_exp_mult_e10_e, 2);
+  rb_define_module_function(module, "expm1",  rb_gsl_sf_expm1, 1);
+  rb_define_module_function(module, "expm1_e",  rb_gsl_sf_expm1_e, 1);
+  rb_define_module_function(module, "exprel",  rb_gsl_sf_exprel, 1);
+  rb_define_module_function(module, "exprel_e",  rb_gsl_sf_exprel_e, 1);
+  rb_define_module_function(module, "exprel_2",  rb_gsl_sf_exprel_2, 1);
+  rb_define_module_function(module, "exprel_2_e",  rb_gsl_sf_exprel_2_e, 1);
+  rb_define_module_function(module, "exprel_n",  rb_gsl_sf_exprel_n, 2);
+  rb_define_module_function(module, "exprel_n_e",  rb_gsl_sf_exprel_n_e, 2);
+  rb_define_module_function(module, "exp_err_e",  rb_gsl_sf_exp_err_e, 2);
+  rb_define_module_function(module, "exp_err_e10_e",  rb_gsl_sf_exp_err_e10_e, 2);
+  rb_define_module_function(module, "exp_mult_err_e",  rb_gsl_sf_exp_mult_err_e, 4);
+  rb_define_module_function(module, "exp_mult_err_e10_e",  rb_gsl_sf_exp_mult_err_e10_e, 4);
+}