rb-gsl 1.15.3.1

data/ext/sf_legendre.c

@@ -0,0 +1,368 @@
+/*
+  legendre.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"
+EXTERN VALUE cgsl_vector;
+
+static VALUE rb_gsl_sf_legendre_P1(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_legendre_P1, x);
+}
+
+static VALUE rb_gsl_sf_legendre_P1_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_legendre_P1_e, x);
+}
+
+static VALUE rb_gsl_sf_legendre_P2(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_legendre_P2, x);
+}
+
+static VALUE rb_gsl_sf_legendre_P2_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_legendre_P2_e, x);
+}
+
+static VALUE rb_gsl_sf_legendre_P3(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_legendre_P3, x);
+}
+
+static VALUE rb_gsl_sf_legendre_P3_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_legendre_P3_e, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Pl(VALUE obj, VALUE l, VALUE x)
+{
+  return rb_gsl_sf_eval_int_double(gsl_sf_legendre_Pl, l, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Pl_e(VALUE obj, VALUE l, VALUE x)
+{
+  return rb_gsl_sf_eval_e_int_double(gsl_sf_legendre_Pl_e, l, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Pl_array(VALUE obj, VALUE lmax, VALUE x)
+{
+  gsl_vector *v = NULL;
+  CHECK_FIXNUM(lmax);
+  Need_Float(x);
+  v = gsl_vector_alloc(FIX2INT(lmax) + 1);
+  gsl_sf_legendre_Pl_array(FIX2INT(lmax), NUM2DBL(x), v->data);
+  return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
+}
+
+static VALUE rb_gsl_sf_legendre_Q0(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_legendre_Q0, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Q0_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_legendre_Q0_e, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Q1(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_legendre_Q1, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Q1_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_legendre_Q1_e, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Ql(VALUE obj, VALUE l, VALUE x)
+{
+  return rb_gsl_sf_eval_int_double(gsl_sf_legendre_Ql, l, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Ql_e(VALUE obj, VALUE l, VALUE x)
+{
+  return rb_gsl_sf_eval_e_int_double(gsl_sf_legendre_Ql_e, l, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Plm(VALUE obj, VALUE l, VALUE m, VALUE x)
+{
+  return rb_gsl_sf_eval_int_int_double(gsl_sf_legendre_Plm, l, m, x);
+}
+
+static VALUE rb_gsl_sf_legendre_Plm_e(VALUE obj, VALUE l, VALUE m, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  int status;
+  CHECK_FIXNUM(l);
+  CHECK_FIXNUM(m);
+  Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  status = gsl_sf_legendre_Plm_e(FIX2INT(l), FIX2INT(m), NUM2DBL(x), rslt);
+  return rb_ary_new3(2, v, INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_legendre_Plm_array(VALUE obj, VALUE lmax, VALUE m, VALUE x)
+{
+  gsl_vector *v = NULL;
+  int size;
+  int ll, mm;
+  CHECK_FIXNUM(lmax); CHECK_FIXNUM(m);
+  Need_Float(x);
+  ll = FIX2INT(lmax);
+  mm = FIX2INT(m);
+  size = gsl_sf_legendre_array_size(ll, mm);
+  v = gsl_vector_alloc(size);
+  gsl_sf_legendre_Plm_array(ll, mm, NUM2DBL(x), v->data);
+  return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
+}
+
+static VALUE rb_gsl_sf_legendre_sphPlm(VALUE obj, VALUE l, VALUE m, VALUE x)
+{
+  return rb_gsl_sf_eval_int_int_double(gsl_sf_legendre_sphPlm, l, m, x);
+}
+
+static VALUE rb_gsl_sf_legendre_sphPlm_e(VALUE obj, VALUE l, VALUE m, VALUE x)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  int status;
+  CHECK_FIXNUM(l); CHECK_FIXNUM(m);
+  Need_Float(x);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  status = gsl_sf_legendre_sphPlm_e(FIX2INT(l), FIX2INT(m), NUM2DBL(x), rslt);
+  return rb_ary_new3(2, v, INT2FIX(status));
+}
+
+static VALUE rb_gsl_sf_legendre_sphPlm_array(VALUE obj, VALUE lmax, VALUE m, VALUE x)
+{
+  gsl_vector *v = NULL;
+  int size;
+  int ll, mm;
+  CHECK_FIXNUM(lmax); CHECK_FIXNUM(m);
+  Need_Float(x);
+  ll = FIX2INT(lmax);
+  mm = FIX2INT(m);
+  size = gsl_sf_legendre_array_size(ll, mm);
+  v = gsl_vector_alloc(size);
+  gsl_sf_legendre_sphPlm_array(ll, mm, NUM2DBL(x), v->data);
+  return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
+}
+
+static VALUE rb_gsl_sf_legendre_array_size(VALUE obj, VALUE lmax, VALUE m)
+{
+  CHECK_FIXNUM(lmax);   CHECK_FIXNUM(m);
+  return INT2FIX(gsl_sf_legendre_array_size(FIX2INT(lmax), FIX2INT(m)));
+}
+
+static VALUE rb_gsl_sf_conicalP_half(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_conicalP_half, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_half_e(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_conicalP_half_e, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_mhalf(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_conicalP_mhalf, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_mhalf_e(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_conicalP_mhalf_e, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_0(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_conicalP_0, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_0_e(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_conicalP_0_e, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_1(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_double_double(gsl_sf_conicalP_1, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_1_e(VALUE obj, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_conicalP_1_e, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_sph_reg(VALUE obj, VALUE l, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_int_double_double(gsl_sf_conicalP_sph_reg, l, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_sph_reg_e(VALUE obj, VALUE l, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_int_double2(gsl_sf_conicalP_sph_reg_e, l, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_cyl_reg(VALUE obj, VALUE m, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_int_double_double(gsl_sf_conicalP_cyl_reg, m, lambda, x);
+}
+
+static VALUE rb_gsl_sf_conicalP_cyl_reg_e(VALUE obj, VALUE m, VALUE lambda, VALUE x)
+{
+  return rb_gsl_sf_eval_e_int_double2(gsl_sf_conicalP_cyl_reg_e, m, lambda, x);
+}
+
+static VALUE rb_gsl_sf_legendre_H3d_0(VALUE obj, VALUE lambda, VALUE eta)
+{
+  Need_Float(lambda); Need_Float(eta);
+  return rb_float_new(gsl_sf_legendre_H3d_0(NUM2DBL(lambda), NUM2DBL(eta)));
+}
+
+static VALUE rb_gsl_sf_legendre_H3d_0_e(VALUE obj, VALUE lambda, VALUE eta)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_legendre_H3d_0_e, lambda, eta);
+}
+
+static VALUE rb_gsl_sf_legendre_H3d_1(VALUE obj, VALUE lambda, VALUE eta)
+{
+  Need_Float(lambda); Need_Float(eta);
+  return rb_float_new(gsl_sf_legendre_H3d_1(NUM2DBL(lambda), NUM2DBL(eta)));
+}
+
+static VALUE rb_gsl_sf_legendre_H3d_1_e(VALUE obj, VALUE lambda, VALUE eta)
+{
+  return rb_gsl_sf_eval_e_double2(gsl_sf_legendre_H3d_1_e, lambda, eta);
+}
+
+
+static VALUE rb_gsl_sf_legendre_H3d(VALUE obj, VALUE l, VALUE lambda, VALUE eta)
+{
+  return rb_float_new(gsl_sf_legendre_H3d(FIX2INT(l), NUM2DBL(lambda), NUM2DBL(eta)));
+}
+
+static VALUE rb_gsl_sf_legendre_H3d_e(VALUE obj,VALUE l,  VALUE lambda, VALUE eta)
+{
+  gsl_sf_result *rslt = NULL;
+  VALUE v;
+  // local variable "status" declared and set, but never used
+  //int status;
+  CHECK_FIXNUM(l);
+  Need_Float(lambda);
+  Need_Float(eta);
+  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
+  /*status =*/ gsl_sf_legendre_H3d_e(FIX2INT(l), NUM2DBL(lambda), NUM2DBL(eta), rslt);
+  return v;
+}
+
+static VALUE rb_gsl_sf_legendre_H3d_array(VALUE obj, VALUE lmax, VALUE lambda, VALUE eta)
+{
+  gsl_vector *v = NULL;
+  CHECK_FIXNUM(lmax);
+  Need_Float(lambda);
+  Need_Float(eta);
+  v = gsl_vector_alloc(FIX2INT(lmax) + 1);
+  gsl_sf_legendre_H3d_array(FIX2INT(lmax), NUM2DBL(lambda), NUM2DBL(eta), v->data);
+  return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
+}
+
+void Init_gsl_sf_legendre(VALUE module)
+{
+  VALUE mgsl_sf_leg;
+
+  rb_define_module_function(module, "legendre_P1",  rb_gsl_sf_legendre_P1, 1);
+  rb_define_module_function(module, "legendre_P1_e",  rb_gsl_sf_legendre_P1_e, 1);
+  rb_define_module_function(module, "legendre_P2",  rb_gsl_sf_legendre_P2, 1);
+  rb_define_module_function(module, "legendre_P2_e",  rb_gsl_sf_legendre_P2_e, 1);
+  rb_define_module_function(module, "legendre_P3",  rb_gsl_sf_legendre_P3, 1);
+  rb_define_module_function(module, "legendre_P3_e",  rb_gsl_sf_legendre_P3_e, 1);
+  rb_define_module_function(module, "legendre_Pl",  rb_gsl_sf_legendre_Pl, 2);
+  rb_define_module_function(module, "legendre_Pl_e",  rb_gsl_sf_legendre_Pl_e, 2);
+  rb_define_module_function(module, "legendre_Pl_array",  rb_gsl_sf_legendre_Pl_array, 2);
+  rb_define_module_function(module, "legendre_Q0",  rb_gsl_sf_legendre_Q0, 1);
+  rb_define_module_function(module, "legendre_Q0_e",  rb_gsl_sf_legendre_Q0_e, 1);
+  rb_define_module_function(module, "legendre_Q1",  rb_gsl_sf_legendre_Q1, 1);
+  rb_define_module_function(module, "legendre_Q1_e",  rb_gsl_sf_legendre_Q1_e, 1);
+  rb_define_module_function(module, "legendre_Ql",  rb_gsl_sf_legendre_Ql, 2);
+  rb_define_module_function(module, "legendre_Ql_e",  rb_gsl_sf_legendre_Ql_e, 2);
+  rb_define_module_function(module, "legendre_Plm",  rb_gsl_sf_legendre_Plm, 3);
+  rb_define_module_function(module, "legendre_Plm_e",  rb_gsl_sf_legendre_Plm_e, 3);
+  rb_define_module_function(module, "legendre_Plm_array",  rb_gsl_sf_legendre_Plm_array, 3);
+  rb_define_module_function(module, "legendre_sphPlm",  rb_gsl_sf_legendre_sphPlm, 3);
+  rb_define_module_function(module, "legendre_sphPlm_e",  rb_gsl_sf_legendre_sphPlm_e, 3);
+  rb_define_module_function(module, "legendre_sphPlm_array",  rb_gsl_sf_legendre_sphPlm_array, 3);
+  rb_define_module_function(module, "legendre_array_size",  rb_gsl_sf_legendre_array_size, 2);
+  rb_define_module_function(module, "conicalP_half",  rb_gsl_sf_conicalP_half, 2);
+  rb_define_module_function(module, "conicalP_half_e",  rb_gsl_sf_conicalP_half_e, 2);
+  rb_define_module_function(module, "conicalP_mhalf",  rb_gsl_sf_conicalP_mhalf, 2);
+  rb_define_module_function(module, "conicalP_mhalf_e",  rb_gsl_sf_conicalP_mhalf_e, 2);
+  rb_define_module_function(module, "conicalP_0",  rb_gsl_sf_conicalP_0, 2);
+  rb_define_module_function(module, "conicalP_0_e",  rb_gsl_sf_conicalP_0_e, 2);
+  rb_define_module_function(module, "conicalP_1",  rb_gsl_sf_conicalP_1, 2);
+  rb_define_module_function(module, "conicalP_1_e",  rb_gsl_sf_conicalP_1_e, 2);
+  rb_define_module_function(module, "conicalP_sph_reg",  rb_gsl_sf_conicalP_sph_reg, 3);
+  rb_define_module_function(module, "conicalP_sph_reg_e",  rb_gsl_sf_conicalP_sph_reg_e, 3);
+  rb_define_module_function(module, "conicalP_cyl_reg",  rb_gsl_sf_conicalP_cyl_reg, 3);
+  rb_define_module_function(module, "conicalP_cyl_reg_e",  rb_gsl_sf_conicalP_cyl_reg_e, 3);
+  rb_define_module_function(module, "legendre_H3d_0",  rb_gsl_sf_legendre_H3d_0, 2);
+  rb_define_module_function(module, "legendre_H3d_0_e",  rb_gsl_sf_legendre_H3d_0_e, 2);
+  rb_define_module_function(module, "legendre_H3d_1",  rb_gsl_sf_legendre_H3d_1, 2);
+  rb_define_module_function(module, "legendre_H3d_1_e",  rb_gsl_sf_legendre_H3d_1_e, 2);
+  rb_define_module_function(module, "legendre_H3d",  rb_gsl_sf_legendre_H3d, 3);
+  rb_define_module_function(module, "legendre_H3d_e",  rb_gsl_sf_legendre_H3d_e, 3);
+  rb_define_module_function(module, "legendre_H3d_array",  rb_gsl_sf_legendre_H3d_array, 3);
+
+  /*****/
+
+  mgsl_sf_leg = rb_define_module_under(module, "Legendre");
+  rb_define_module_function(mgsl_sf_leg, "P1",  rb_gsl_sf_legendre_P1, 1);
+  rb_define_module_function(mgsl_sf_leg, "P1_e",  rb_gsl_sf_legendre_P1_e, 1);
+  rb_define_module_function(mgsl_sf_leg, "P2",  rb_gsl_sf_legendre_P2, 1);
+  rb_define_module_function(mgsl_sf_leg, "P2_e",  rb_gsl_sf_legendre_P2_e, 1);
+  rb_define_module_function(mgsl_sf_leg, "P3",  rb_gsl_sf_legendre_P3, 1);
+  rb_define_module_function(mgsl_sf_leg, "P3_e",  rb_gsl_sf_legendre_P3_e, 1);
+  rb_define_module_function(mgsl_sf_leg, "Pl",  rb_gsl_sf_legendre_Pl, 2);
+  rb_define_module_function(mgsl_sf_leg, "Pl_e",  rb_gsl_sf_legendre_Pl_e, 2);
+  rb_define_module_function(mgsl_sf_leg, "Pl_array",  rb_gsl_sf_legendre_Pl_array, 2);
+  rb_define_module_function(mgsl_sf_leg, "Q0",  rb_gsl_sf_legendre_Q0, 1);
+  rb_define_module_function(mgsl_sf_leg, "Q0_e",  rb_gsl_sf_legendre_Q0_e, 1);
+  rb_define_module_function(mgsl_sf_leg, "Q1",  rb_gsl_sf_legendre_Q1, 1);
+  rb_define_module_function(mgsl_sf_leg, "Q1_e",  rb_gsl_sf_legendre_Q1_e, 1);
+  rb_define_module_function(mgsl_sf_leg, "Plm",  rb_gsl_sf_legendre_Plm, 3);
+  rb_define_module_function(mgsl_sf_leg, "Plm_e",  rb_gsl_sf_legendre_Plm_e, 3);
+  rb_define_module_function(mgsl_sf_leg, "Plm_array",  rb_gsl_sf_legendre_Plm_array, 3);
+  rb_define_module_function(mgsl_sf_leg, "sphPlm",  rb_gsl_sf_legendre_sphPlm, 3);
+  rb_define_module_function(mgsl_sf_leg, "sphPlm_e",  rb_gsl_sf_legendre_sphPlm_e, 3);
+  rb_define_module_function(mgsl_sf_leg, "sphPlm_array",  rb_gsl_sf_legendre_sphPlm_array, 3);
+  rb_define_module_function(mgsl_sf_leg, "array_size",  rb_gsl_sf_legendre_array_size, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_half",  rb_gsl_sf_conicalP_half, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_half_e",  rb_gsl_sf_conicalP_half_e, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_mhalf",  rb_gsl_sf_conicalP_mhalf, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_mhalf_e",  rb_gsl_sf_conicalP_mhalf_e, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_0",  rb_gsl_sf_conicalP_0, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_0_e",  rb_gsl_sf_conicalP_0_e, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_1",  rb_gsl_sf_conicalP_1, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_1_e",  rb_gsl_sf_conicalP_1_e, 2);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_sph_reg",  rb_gsl_sf_conicalP_sph_reg, 3);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_sph_reg_e",  rb_gsl_sf_conicalP_sph_reg_e, 3);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_cyl_reg",  rb_gsl_sf_conicalP_cyl_reg, 3);
+  rb_define_module_function(mgsl_sf_leg, "conicalP_cyl_reg_e",  rb_gsl_sf_conicalP_cyl_reg_e, 3);
+  rb_define_module_function(mgsl_sf_leg, "H3d_0",  rb_gsl_sf_legendre_H3d_0, 2);
+  rb_define_module_function(mgsl_sf_leg, "H3d_0_e",  rb_gsl_sf_legendre_H3d_0_e, 2);
+  rb_define_module_function(mgsl_sf_leg, "H3d_1",  rb_gsl_sf_legendre_H3d_1, 2);
+  rb_define_module_function(mgsl_sf_leg, "H3d_1_e",  rb_gsl_sf_legendre_H3d_1_e, 2);
+  rb_define_module_function(mgsl_sf_leg, "H3d",  rb_gsl_sf_legendre_H3d, 3);
+  rb_define_module_function(mgsl_sf_leg, "H3d_e",  rb_gsl_sf_legendre_H3d_e, 3);
+  rb_define_module_function(mgsl_sf_leg, "H3d_array",  rb_gsl_sf_legendre_H3d_array, 3);
+
+}

data/ext/sf_log.c

@@ -0,0 +1,105 @@
+/*
+  sf_log.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_log(VALUE obj, VALUE x)
+{
+  if (COMPLEX_P(x) || VECTOR_COMPLEX_P(x) || MATRIX_COMPLEX_P(x)) 
+    return rb_gsl_math_complex_eval(gsl_complex_log, x);
+  return rb_gsl_sf_eval1(gsl_sf_log, x);
+}
+
+static VALUE rb_gsl_sf_log10(VALUE obj, VALUE x)
+{
+  if (COMPLEX_P(x) || VECTOR_COMPLEX_P(x) || MATRIX_COMPLEX_P(x)) 
+    return rb_gsl_math_complex_eval(gsl_complex_log10, x);
+  return rb_gsl_sf_eval1(log10, x);
+}
+
+
+static VALUE rb_gsl_sf_log_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_log_e, x);
+}
+
+static VALUE rb_gsl_sf_log_abs(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_log_abs, x);
+}
+
+static VALUE rb_gsl_sf_log_abs_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_log_abs_e, x);
+}
+
+static VALUE rb_gsl_sf_complex_log_e(int argc, VALUE *argv, VALUE obj)
+{
+  gsl_sf_result *rslt1 = NULL, *rslt2 = NULL;
+  gsl_complex *z = NULL;
+  VALUE vlnr, vtheta;
+  double re, im;
+  // local variable "status" was defined and set, but never used
+  //int status;
+  switch (argc) {
+  case 1:
+    CHECK_COMPLEX(argv[0]);
+    Data_Get_Struct(argv[0], gsl_complex, z);
+    re = GSL_REAL(*z);
+    im = GSL_IMAG(*z);
+    break;
+  case 2:
+    Need_Float(argv[0]);     Need_Float(argv[1]);
+    re = NUM2DBL(argv[0]);
+    im = NUM2DBL(argv[1]);
+  default:
+    rb_raise(rb_eArgError, "wrong number of arguments (%d for 1 or 2)", argc);
+    break;
+  }
+  vlnr = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt1);
+  vtheta = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt2);
+  /*status =*/ gsl_sf_complex_log_e(re, im, rslt1, rslt2);
+  return rb_ary_new3(2, vlnr, vtheta);
+}
+
+static VALUE rb_gsl_sf_log_1plusx(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_log_1plusx, x);
+}
+
+static VALUE rb_gsl_sf_log_1plusx_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_log_1plusx_e, x);
+}
+
+static VALUE rb_gsl_sf_log_1plusx_mx(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval1(gsl_sf_log_1plusx_mx, x);
+}
+
+static VALUE rb_gsl_sf_log_1plusx_mx_e(VALUE obj, VALUE x)
+{
+  return rb_gsl_sf_eval_e(gsl_sf_log_1plusx_mx_e, x);
+}
+
+void Init_gsl_sf_log(VALUE module)
+{
+  rb_define_module_function(module, "log",  rb_gsl_sf_log, 1);
+  rb_define_module_function(module, "log10",  rb_gsl_sf_log10, 1);  
+  rb_define_module_function(module, "log_e",  rb_gsl_sf_log_e, 1);
+  rb_define_module_function(module, "log_abs",  rb_gsl_sf_log_abs, 1);
+  rb_define_module_function(module, "log_abs_e",  rb_gsl_sf_log_abs_e, 1);
+  rb_define_module_function(module, "complex_log_e",  rb_gsl_sf_complex_log_e, -1);
+  rb_define_module_function(module, "log_1plusx",  rb_gsl_sf_log_1plusx, 1);
+  rb_define_module_function(module, "log_1plusx_e",  rb_gsl_sf_log_1plusx_e, 1);
+  rb_define_module_function(module, "log_1plusx_mx",  rb_gsl_sf_log_1plusx_mx, 1);
+  rb_define_module_function(module, "log_1plusx_mx_e",  rb_gsl_sf_log_1plusx_mx_e, 1);
+}