rb-gsl 1.16.0.5 → 1.16.0.6

data/ext/gsl_native/error.c

@@ -1,193 +0,0 @@
-/*
-  error.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.h"
-#include <gsl/gsl_errno.h>
-#include "include/rb_gsl_array.h"
-#include "include/rb_gsl_common.h"
-#include "include/rb_gsl_function.h"
-
-static VALUE eHandler;
-static VALUE cgsl_error[35];
-static VALUE *pgsl_error;
-
-static void Init_rb_gsl_define_GSL_CONST(VALUE module);
-void rb_gsl_error_handler(const char *reason, const char *file,
-                          int line, int gsl_errno);
-static void rb_gsl_my_error_handler(const char *reason, const char *file,
-                                    int line, int gsl_errno);
-
-void rb_gsl_error_handler(const char *reason, const char *file,
-                          int line, int gsl_errno)
-{
-  const char *emessage = gsl_strerror(gsl_errno);
-  rb_raise(pgsl_error[gsl_errno],
-           "Ruby/GSL error code %d, %s (file %s, line %d), %s",
-           gsl_errno, reason, file, line, emessage);
-}
-
-static void rb_gsl_my_error_handler(const char *reason, const char *file,
-                                    int line, int gsl_errno)
-{
-  VALUE vreason, vfile;
-  VALUE vline, verrno;
-  vreason = rb_str_new2(reason);
-  vfile = rb_str_new2(file);
-  vline = INT2FIX(line);
-  verrno = INT2FIX(gsl_errno);
-  rb_funcall(eHandler, RBGSL_ID_call, 4, vreason, vfile, vline, verrno);
-}
-
-static VALUE rb_gsl_set_error_handler(int argc, VALUE *argv, VALUE module)
-{
-  if (rb_block_given_p()) {
-    eHandler = rb_block_proc();
-    gsl_set_error_handler(&rb_gsl_my_error_handler);
-    return Qtrue;
-  }
-  switch (argc) {
-  case 0:
-    gsl_set_error_handler(&rb_gsl_error_handler);
-    return Qtrue;
-    break;
-  case 1:
-    CHECK_PROC(argv[0]);
-    eHandler = argv[0];
-    gsl_set_error_handler(&rb_gsl_my_error_handler);
-    return Qtrue;
-    break;
-  default:
-    rb_raise(rb_eArgError, "too many arguments (%d for 0 or 1 Proc)", argc);
-    break;
-  }
-}
-
-static VALUE rb_gsl_set_default_error_handler(VALUE module)
-{
-  gsl_set_error_handler(&rb_gsl_error_handler);
-  return Qtrue;
-}
-
-static void rb_gsl_define_exceptions(VALUE module)
-{
-  VALUE mgsl_error;
-  mgsl_error = rb_define_module_under(module, "ERROR");
-  pgsl_error = &cgsl_error[2];
-  pgsl_error[-2] = rb_define_class_under(mgsl_error, "CONTINUE", rb_cFixnum);
-  pgsl_error[-1] = rb_define_class_under(mgsl_error, "FAILURE", rb_eRuntimeError);
-  pgsl_error[0] = rb_define_class_under(mgsl_error, "SUCCESS", rb_cFixnum);
-  pgsl_error[1] = rb_define_class_under(mgsl_error, "EDOM", rb_eRangeError);
-  pgsl_error[2] = rb_define_class_under(mgsl_error, "ERANGE", rb_eRangeError);
-  pgsl_error[3] = rb_define_class_under(mgsl_error, "EFAULT", rb_eRuntimeError);
-  pgsl_error[4] = rb_define_class_under(mgsl_error, "EINVAL", rb_eIndexError);
-  pgsl_error[5] = rb_define_class_under(mgsl_error, "EFAILED", rb_eRuntimeError);
-  pgsl_error[6] = rb_define_class_under(mgsl_error, "EFACTOR", rb_eRuntimeError);
-  pgsl_error[7] = rb_define_class_under(mgsl_error, "ESANITY", rb_eRuntimeError);
-  pgsl_error[8] = rb_define_class_under(mgsl_error, "ENOMEM", rb_eNoMemError);
-  pgsl_error[9] = rb_define_class_under(mgsl_error, "EBADFUNC", rb_eRuntimeError);
-  pgsl_error[10] = rb_define_class_under(mgsl_error, "ERUNAWAY", rb_eRuntimeError);
-  pgsl_error[11] = rb_define_class_under(mgsl_error, "EMAXITER", rb_eRuntimeError);
-  pgsl_error[12] = rb_define_class_under(mgsl_error, "EZERODIV", rb_eZeroDivError);
-  pgsl_error[13] = rb_define_class_under(mgsl_error, "EBADTOL", rb_eRuntimeError);
-  pgsl_error[14] = rb_define_class_under(mgsl_error, "ETOL", rb_eRuntimeError);
-  pgsl_error[15] = rb_define_class_under(mgsl_error, "EUNDRFLW", rb_eRangeError);
-  pgsl_error[16] = rb_define_class_under(mgsl_error, "EOVRFLW", rb_eRangeError);
-  pgsl_error[17] = rb_define_class_under(mgsl_error, "ELOSS", rb_eRuntimeError);
-  pgsl_error[18] = rb_define_class_under(mgsl_error, "EROUND", rb_eRuntimeError);
-  pgsl_error[19] = rb_define_class_under(mgsl_error, "EBADLEN", rb_eIndexError);
-  pgsl_error[20] = rb_define_class_under(mgsl_error, "ENOTSQR", rb_eRuntimeError);
-  pgsl_error[21] = rb_define_class_under(mgsl_error, "ESING", rb_eRuntimeError);
-  pgsl_error[22] = rb_define_class_under(mgsl_error, "EDIVERGE", rb_eRuntimeError);
-  pgsl_error[23] = rb_define_class_under(mgsl_error, "EUNSUP", rb_eRuntimeError);
-  pgsl_error[24] = rb_define_class_under(mgsl_error, "EUNIMPL", rb_eNotImpError);
-  pgsl_error[25] = rb_define_class_under(mgsl_error, "ECACHE", rb_eRuntimeError);
-  pgsl_error[26] = rb_define_class_under(mgsl_error, "ETABLE", rb_eRuntimeError);
-  pgsl_error[27] = rb_define_class_under(mgsl_error, "ENOPROG", rb_eRuntimeError);
-  pgsl_error[28] = rb_define_class_under(mgsl_error, "ENOPROGJ", rb_eRuntimeError);
-  pgsl_error[29] = rb_define_class_under(mgsl_error, "ETOLF", rb_eRuntimeError);
-  pgsl_error[30] = rb_define_class_under(mgsl_error, "ETOLX", rb_eRuntimeError);
-  pgsl_error[31] = rb_define_class_under(mgsl_error, "ETOLG", rb_eRuntimeError);
-  pgsl_error[32] = rb_define_class_under(mgsl_error, "EOF", rb_eEOFError);
-}
-
-static void Init_rb_gsl_define_GSL_CONST(VALUE module)
-{
-  rb_define_const(module, "SUCCESS", INT2FIX(GSL_SUCCESS));
-  rb_define_const(module, "FAILURE", INT2FIX(GSL_FAILURE));
-  rb_define_const(module, "CONTINUE", INT2FIX(GSL_CONTINUE));
-  rb_define_const(module, "EDOM", INT2FIX(GSL_EDOM));
-  rb_define_const(module, "ERANGE", INT2FIX(GSL_ERANGE));
-  rb_define_const(module, "EFAULT", INT2FIX(GSL_EFAULT));
-  rb_define_const(module, "EINVAL", INT2FIX(GSL_EINVAL));
-  rb_define_const(module, "EFAILED", INT2FIX(GSL_EFAILED));
-  rb_define_const(module, "EFACTOR", INT2FIX(GSL_EFACTOR));
-  rb_define_const(module, "ESANITY", INT2FIX(GSL_ESANITY));
-  rb_define_const(module, "ENOMEM", INT2FIX(GSL_ENOMEM));
-  rb_define_const(module, "EBADFUNC", INT2FIX(GSL_EBADFUNC));
-  rb_define_const(module, "ERUNAWAY", INT2FIX(GSL_ERUNAWAY));
-  rb_define_const(module, "EMAXITER", INT2FIX(GSL_EMAXITER));
-  rb_define_const(module, "EZERODIV", INT2FIX(GSL_EZERODIV));
-  rb_define_const(module, "EBADTOL", INT2FIX(GSL_EBADTOL));
-  rb_define_const(module, "ETOL", INT2FIX(GSL_ETOL));
-  rb_define_const(module, "EUNDRFLW", INT2FIX(GSL_EUNDRFLW));
-  rb_define_const(module, "EOVRFLW", INT2FIX(GSL_EOVRFLW));
-  rb_define_const(module, "ELOSS", INT2FIX(GSL_ELOSS));
-  rb_define_const(module, "EROUND", INT2FIX(GSL_EROUND));
-  rb_define_const(module, "EBADLEN", INT2FIX(GSL_EBADLEN));
-  rb_define_const(module, "ENOTSQR", INT2FIX(GSL_ENOTSQR));
-  rb_define_const(module, "ESING", INT2FIX(GSL_ESING));
-  rb_define_const(module, "EDIVERGE", INT2FIX(GSL_EDIVERGE));
-  rb_define_const(module, "EUNSUP", INT2FIX(GSL_EUNSUP));
-  rb_define_const(module, "EUNIMPL", INT2FIX(GSL_EUNIMPL));
-  rb_define_const(module, "ECACHE", INT2FIX(GSL_ECACHE));
-  rb_define_const(module, "ETABLE", INT2FIX(GSL_ETABLE));
-  rb_define_const(module, "ENOPROG", INT2FIX(GSL_ENOPROG));
-  rb_define_const(module, "ENOPROGJ", INT2FIX(GSL_ENOPROGJ));
-  rb_define_const(module, "ETOLF", INT2FIX(GSL_ETOLF));
-  rb_define_const(module, "ETOLX", INT2FIX(GSL_ETOLX));
-  rb_define_const(module, "ETOLG", INT2FIX(GSL_ETOLG));
-  rb_define_const(module, "EOF", INT2FIX(GSL_EOF));
-}
-
-static VALUE rb_gsl_set_error_handler_off(VALUE module)
-{
-  gsl_set_error_handler_off();
-  return Qtrue;
-}
-
-static void define_module_functions(VALUE module);
-static VALUE rb_gsl_strerror(VALUE obj, VALUE errn);
-static void define_module_functions(VALUE module)
-{
-  rb_define_module_function(module, "set_error_handler_off",
-                            rb_gsl_set_error_handler_off, 0);
-  rb_define_module_function(module, "strerror",
-                            rb_gsl_strerror, 1);
-  rb_define_module_function(module, "set_error_handler",
-                            rb_gsl_set_error_handler, -1);
-  rb_define_module_function(module, "set_default_error_handler",
-                            rb_gsl_set_default_error_handler, 0);
-}
-
-static VALUE rb_gsl_strerror(VALUE obj, VALUE errn)
-{
-  int gsl_errno = FIX2INT(errn);
-  return rb_str_new2(gsl_strerror(gsl_errno));
-}
-
-void Init_gsl_error(VALUE module)
-{
-  Init_rb_gsl_define_GSL_CONST(module);
-
-  gsl_set_error_handler(&rb_gsl_error_handler);
-
-  define_module_functions(module);
-  rb_gsl_define_exceptions(module);
-}

data/ext/gsl_native/extconf.rb

@@ -1,118 +0,0 @@
-require 'mkmf'
-
-def gsl_config_arg(arg)
-  yield arg_config("--with-gsl-#{arg}") {
-    sh = 'sh ' if RUBY_PLATFORM =~ /mingw/
-    IO.popen("#{sh}gsl-config --#{arg}") { |f| f.gets.chomp }
-  }, lambda { |val| puts "checking gsl #{arg}... #{val}"; val }
-rescue => err
-  abort "*** ERROR: missing required library to compile this module: #{err}"
-end
-
-def gsl_def(const, value = nil)
-  value = "=#{value}" if value
-  $defs << "-D#{const}#{value}"
-end
-
-def gsl_have_header(library, header)
-  have_library(library) if have_header(header)
-end
-
-def gsl_have_library(func)
-  have_func(func) if have_library('gsl', func)
-end
-
-def gsl_dir_config(target, idir = nil, ldir = idir)
-  dir_config(target, idir || $sitearchdir, ldir || $sitearchdir)
-end
-
-def gsl_gem_config(target, dir = 'ext')
-  path = begin
-    require 'rubygems'
-
-    spec = Gem::Specification.find_by_path("#{target}.h")
-    File.join(spec.full_gem_path, dir) if spec
-  rescue LoadError
-  end
-
-  gsl_dir_config(target, path)
-
-  $LOCAL_LIBS += " -l:#{target}.so" if arg_config("--force-link-#{target}")   ||
-                                       $CFLAGS.include?('-Wl,--no-undefined') ||
-                                       $LDFLAGS.include?('-Wl,--no-undefined')
-end
-
-$CFLAGS += ' -Wall -Iinclude'
-
-gsl_config_arg(:version) { |version, check|
-  gsl_def(:GSL_VERSION, check[version])
-
-  ver = version.split('.').map { |x| x.to_i }
-
-  later = lambda { |other|
-    ary = other.split('.').map { |x| x.to_i }
-
-    gte = ver[0] > ary[0] ? true  :
-          ver[0] < ary[0] ? false :
-          ver[1] > ary[1] ? true  :
-          ver[1] < ary[1] ? false :
-          ver.size < ary.size ? false :
-          ver.size == 3 && ary.size == 3 ? ver[2] >= ary[2] : true
-
-    ary.pop && ary[-1] += 1 if ary.last == 90
-
-    gte && gsl_def("GSL_#{ary.join('_')}_LATER")
-  }
-
-  raise 'Ruby/GSL requires gsl-1.15 or later.' unless later['1.15']
-
-  %w[1.15 1.16].each { |v| later[v] }
-}
-
-gsl_config_arg(:cflags) { |cflags, check|
-  $CFLAGS += ' ' + check[cflags]
-}
-
-gsl_config_arg(:libs) { |libs, check|
-  libs.tr!(File::PATH_SEPARATOR, ' ')
-
-  if enable_config('atlas')
-    dir_config('cblas')
-    dir_config('atlas')
-
-    if have_library('cblas') && have_library('atlas')
-      libs.gsub!('-lgslcblas', '-lcblas -latlas')
-    end
-  end
-
-  $LOCAL_LIBS += ' ' + check[libs]
-}
-
-have_func('round')
-
-%w[alf qrngextra rngextra tensor].each { |library|
-  gsl_have_header(library, "#{library}/#{library}.h")
-}
-
-gsl_have_header('bundle_method', 'gsl/gsl_multimin_fsdf.h')
-gsl_have_header('cqp',           'gsl/gsl_cqp.h')
-gsl_have_header('jacobi',        'jacobi.h')
-gsl_have_header('ndlinear',      'ndlinear/gsl_multifit_ndlinear.h')
-gsl_have_header('ool',           'ool/ool_version.h')
-
-gsl_have_library('gsl_eigen_francis')
-gsl_have_library('gsl_poly_solve_quartic')
-
-gsl_def(:HAVE_GNU_GRAPH) if find_executable('graph')
-
-gsl_gem_config('narray', 'src')
-
-have_header('narray.h')
-have_library('narray') if RUBY_PLATFORM =~ /cygwin|mingw/
-
-unless arg_config('--disable-tamu-anova')
-  gsl_dir_config('tamu_anova')
-  gsl_have_header('tamuanova', 'tamu_anova/tamu_anova.h')
-end
-
-create_makefile('gsl_native')

data/ext/gsl_native/fft.c

@@ -1,1095 +0,0 @@
-/*
-  fft.c
-  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
-    (C) Copyright 2004 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_fft.h"
-
-VALUE mgsl_fft;
-VALUE cgsl_fft_wavetable;
-VALUE cgsl_fft_complex_wavetable, cgsl_fft_complex_workspace;
-VALUE cgsl_fft_real_wavetable, cgsl_fft_halfcomplex_wavetable;
-VALUE cgsl_fft_real_workspace;
-
-extern VALUE cgsl_vector_int;
-extern VALUE cgsl_vector_complex;
-
-static void GSL_FFT_Workspace_free(GSL_FFT_Workspace *space);
-
-static VALUE rb_gsl_fft_complex_wavetable_new(VALUE klass, VALUE n)
-{
-  CHECK_FIXNUM(n);
-  return Data_Wrap_Struct(cgsl_fft_complex_wavetable, 0,
-                          gsl_fft_complex_wavetable_free,
-                          gsl_fft_complex_wavetable_alloc(FIX2INT(n)));
-}
-
-static VALUE rb_gsl_fft_real_wavetable_new(VALUE klass, VALUE n)
-{
-  CHECK_FIXNUM(n);
-  return Data_Wrap_Struct(klass, 0, gsl_fft_real_wavetable_free,
-                          gsl_fft_real_wavetable_alloc(FIX2INT(n)));
-}
-
-static VALUE rb_gsl_fft_halfcomplex_wavetable_new(VALUE klass, VALUE n)
-{
-  CHECK_FIXNUM(n);
-  return Data_Wrap_Struct(klass, 0, gsl_fft_halfcomplex_wavetable_free,
-                          gsl_fft_halfcomplex_wavetable_alloc(FIX2INT(n)));
-
-}
-
-static void GSL_FFT_Wavetable_free(GSL_FFT_Wavetable *table)
-{
-  gsl_fft_complex_wavetable_free((gsl_fft_complex_wavetable *) table);
-}
-
-static VALUE rb_gsl_fft_complex_workspace_new(VALUE klass, VALUE n)
-{
-  CHECK_FIXNUM(n);
-  return Data_Wrap_Struct(klass, 0, gsl_fft_complex_workspace_free,
-                          gsl_fft_complex_workspace_alloc(FIX2INT(n)));
-}
-
-static VALUE rb_gsl_fft_real_workspace_new(VALUE klass, VALUE n)
-{
-  CHECK_FIXNUM(n);
-  return Data_Wrap_Struct(klass, 0, gsl_fft_real_workspace_free,
-                          gsl_fft_real_workspace_alloc(FIX2INT(n)));
-}
-
-static void GSL_FFT_Workspace_free(GSL_FFT_Workspace *space)
-{
-  gsl_fft_complex_workspace_free((gsl_fft_complex_workspace *) space);
-}
-
-// The FFT methods used to allow passing stride and n values as optional
-// parameters to control which elements get transformed.  This created problems
-// for Views which can have their own stride, so support for stride and n
-// parameters to the transform methods is being dropped.  This method used to
-// be called to determine the stride and n values to use based on the
-// parameters and/or the vector itself (depending on how many parameters were
-// passed). Now this function is somewhat unneceesary, but to simplify the code
-// refactoring, it has been left in place for the time being.  Eventually it
-// can be refactored away completely.
-static VALUE get_complex_stride_n(VALUE obj,
-                                  gsl_vector_complex **vin,
-                                  gsl_complex_packed_array *data, size_t *stride, size_t *n)
-{
-  gsl_vector_complex *v = NULL;
-
-  // obj must be a GSL::Vector::Complex
-  CHECK_VECTOR_COMPLEX(obj);
-  Data_Get_Struct(obj, gsl_vector_complex, v);
-
-  if(vin) *vin = v;
-  *data = (gsl_complex_packed_array) v->data;
-  *stride = v->stride;
-  *n = v->size;
-  return obj;
-}
-
-static VALUE rb_fft_complex_radix2(VALUE obj,
-                                   int (*trans)(gsl_complex_packed_array,
-                                                size_t, size_t), int flag)
-{
-  size_t stride, n;
-  gsl_complex_packed_array data;
-  gsl_vector_complex *vin, *vout;
-  VALUE ary;
-  ary = get_complex_stride_n(obj, &vin, &data, &stride, &n);
-  if (flag == RB_GSL_FFT_COPY) {
-    vout = gsl_vector_complex_alloc(n);
-    gsl_vector_complex_memcpy(vout, vin);
-    (*trans)(vout->data, vout->stride /*1*/, vout->size /*n*/);
-    return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vout);
-  } else { /* in-place */
-    (*trans)(data, stride, n);
-    return ary;
-  }
-}
-
-static VALUE rb_gsl_fft_complex_radix2_forward(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj, gsl_fft_complex_radix2_forward,
-                               RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_forward2(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj, gsl_fft_complex_radix2_forward,
-                               RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_transform(VALUE obj, VALUE val_sign)
-{
-  size_t stride, n;
-  gsl_complex_packed_array data;
-  gsl_fft_direction sign;
-  gsl_vector_complex *vin, *vout;
-  sign = NUM2INT(val_sign);
-  get_complex_stride_n(obj, &vin, &data, &stride, &n);
-  vout = gsl_vector_complex_alloc(n);
-  gsl_vector_complex_memcpy(vout, vin);
-  gsl_fft_complex_radix2_transform(vout->data, vout->stride /*1*/, vout->size /*n*/, sign);
-  return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vout);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_transform2(VALUE obj, VALUE val_sign)
-{
-  size_t stride, n;
-  gsl_complex_packed_array data;
-  gsl_fft_direction sign;
-  VALUE ary;
-  sign = NUM2INT(val_sign);
-  ary = get_complex_stride_n(obj, NULL, &data, &stride, &n);
-  gsl_fft_complex_radix2_transform(data, stride, n, sign);
-  return ary;
-}
-
-static VALUE rb_gsl_fft_complex_radix2_backward(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj, gsl_fft_complex_radix2_backward,
-                               RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_inverse(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj, gsl_fft_complex_radix2_inverse,
-                               RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_dif_forward(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj,
-                               gsl_fft_complex_radix2_dif_forward,
-                               RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_backward2(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj, gsl_fft_complex_radix2_backward,
-                               RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_inverse2(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj, gsl_fft_complex_radix2_inverse,
-                               RB_GSL_FFT_INPLACE);
-}
-
-
-static VALUE rb_gsl_fft_complex_radix2_dif_forward2(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj,
-                               gsl_fft_complex_radix2_dif_forward,
-                               RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_dif_transform(VALUE obj, VALUE val_sign)
-{
-  size_t stride, n;
-  gsl_complex_packed_array data;
-  gsl_vector_complex *vin, *vout;
-  gsl_fft_direction sign;
-  sign = NUM2INT(val_sign);
-  get_complex_stride_n(obj, &vin, &data, &stride, &n);
-  vout = gsl_vector_complex_alloc(n);
-  gsl_vector_complex_memcpy(vout, vin);
-  gsl_fft_complex_radix2_dif_transform(vout->data, vout->stride /*1*/, vout->size /*n*/, sign);
-  return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vout);
-}
-
-/* in-place */
-static VALUE rb_gsl_fft_complex_radix2_dif_transform2(VALUE obj, VALUE val_sign)
-{
-  size_t stride, n;
-  gsl_complex_packed_array data;
-  gsl_fft_direction sign;
-  VALUE ary;
-  sign = NUM2INT(val_sign);
-  ary = get_complex_stride_n(obj, NULL, &data, &stride, &n);
-  gsl_fft_complex_radix2_dif_transform(data, stride, n, sign);
-  return ary;
-}
-
-static VALUE rb_gsl_fft_complex_radix2_dif_backward(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj,
-                               gsl_fft_complex_radix2_dif_backward,
-                               RB_GSL_FFT_COPY);
-}
-static VALUE rb_gsl_fft_complex_radix2_dif_inverse(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj,
-                               gsl_fft_complex_radix2_dif_inverse,
-                               RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_complex_radix2_dif_backward2(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj,
-                               gsl_fft_complex_radix2_dif_backward,
-                               RB_GSL_FFT_INPLACE);
-}
-static VALUE rb_gsl_fft_complex_radix2_dif_inverse2(VALUE obj)
-{
-  return rb_fft_complex_radix2(obj,
-                               gsl_fft_complex_radix2_dif_inverse,
-                               RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_GSL_FFT_Wavetable_n(VALUE obj)
-{
-  GSL_FFT_Wavetable *table;
-  Data_Get_Struct(obj, GSL_FFT_Wavetable, table);
-  return INT2FIX(table->n);
-}
-
-static VALUE rb_GSL_FFT_Wavetable_nf(VALUE obj)
-{
-  GSL_FFT_Wavetable *table;
-  Data_Get_Struct(obj, GSL_FFT_Wavetable, table);
-  return INT2FIX(table->nf);
-}
-
-static VALUE rb_GSL_FFT_Wavetable_factor(VALUE obj)
-{
-  GSL_FFT_Wavetable *table;
-  gsl_vector_int *v;
-  size_t i;
-  Data_Get_Struct(obj, GSL_FFT_Wavetable, table);
-  v = gsl_vector_int_alloc(table->nf);
-  for (i = 0; i < table->nf; i++) gsl_vector_int_set(v, i, table->factor[i]);
-  return Data_Wrap_Struct(cgsl_vector_int, 0, gsl_vector_int_free, v);
-}
-
-enum {
-  NONE_OF_TWO = 0,
-  ALLOC_SPACE = 1,
-  ALLOC_TABLE = 2,
-  BOTH_OF_TWO = 3,
-} FFTComplexStructAllocFlag;
-
-static void gsl_fft_free(int flag, GSL_FFT_Wavetable *table,
-                         GSL_FFT_Workspace *space);
-
-// Parse argc, argv.  obj must be GSL::Vector::Complex.
-// This can be simplified at some point.
-// See comments preceding get_complex_stride_n()
-static int gsl_fft_get_argv_complex(int argc, VALUE *argv, VALUE obj,
-                                    gsl_vector_complex ** vin,
-                                    gsl_complex_packed_array *data, size_t *stride,
-                                    size_t *n, gsl_fft_complex_wavetable **table,
-                                    gsl_fft_complex_workspace **space)
-{
-  int flag = NONE_OF_TWO, flagtmp, i, itmp = argc, itmp2 = 0, ccc;
-  int flagw = 0;
-
-  CHECK_VECTOR_COMPLEX(obj);
-
-  ccc = argc;
-  flagtmp = 0;
-  flagw = 0;
-  for (i = argc-1; i >= itmp2; i--) {
-    if (rb_obj_is_kind_of(argv[i], cgsl_fft_complex_workspace)) {
-      Data_Get_Struct(argv[i], gsl_fft_complex_workspace, *space);
-      flagtmp = 1;
-      flagw = 1;
-      itmp = i;
-      ccc--;
-      break;
-    }
-  }
-  flagtmp = 0;
-  for (i = itmp-1; i >= itmp2; i--) {
-    if (rb_obj_is_kind_of(argv[i], cgsl_fft_complex_wavetable)) {
-      Data_Get_Struct(argv[i], gsl_fft_complex_wavetable, *table);
-      flagtmp = 1;
-      ccc--;
-      break;
-    }
-  }
-  get_complex_stride_n(obj, vin, data, stride, n);
-  if (flagw == 0) {
-    *space = gsl_fft_complex_workspace_alloc(*n);
-    flag += ALLOC_SPACE;
-  }
-  if (flagtmp == 0) {
-    *table = gsl_fft_complex_wavetable_alloc(*n);
-    flag += ALLOC_TABLE;
-  }
-  if (*table == NULL) {
-    rb_raise(rb_eRuntimeError, "something wrong with wavetable");
-  }
-  if (*space == NULL) {
-    rb_raise(rb_eRuntimeError, "something wrong with workspace");
-  }
-  return flag;
-}
-
-// Parse argc, argv.  obj must be GSL::Vector of real data
-static int gsl_fft_get_argv_real(int argc, VALUE *argv, VALUE obj,
-                                 double **ptr, size_t *stride,
-                                 size_t *n, gsl_fft_real_wavetable **table,
-                                 gsl_fft_real_workspace **space, int *naflag)
-{
-  int flag = NONE_OF_TWO, flagtmp, i, itmp = argc, itmp2 = 0, ccc;
-  int flagw = 0;
-  *naflag = 0;
-
-  *ptr = get_ptr_double3(obj, n, stride, naflag);
-
-  ccc = argc;
-  flagtmp = 0;
-  flagw = 0;
-  for (i = argc-1; i >= itmp2; i--) {
-    if (rb_obj_is_kind_of(argv[i], cgsl_fft_real_workspace)) {
-      Data_Get_Struct(argv[i], gsl_fft_real_workspace, *space);
-      flagtmp = 1;
-      flagw = 1;
-      itmp = i;
-      ccc--;
-      break;
-    }
-  }
-  flagtmp = 0;
-  for (i = itmp-1; i >= itmp2; i--) {
-    if (rb_obj_is_kind_of(argv[i], cgsl_fft_real_wavetable)) {
-      Data_Get_Struct(argv[i], gsl_fft_real_wavetable, *table);
-      flagtmp = 1;
-      ccc--;
-      break;
-    }
-  }
-  if (flagw == 0) {
-    *space = gsl_fft_real_workspace_alloc(*n);
-    flag += ALLOC_SPACE;
-  }
-  if (flagtmp == 0) {
-    *table = gsl_fft_real_wavetable_alloc(*n);
-    flag += ALLOC_TABLE;
-  }
-  if (*table == NULL) {
-    rb_raise(rb_eRuntimeError, "something wrong with wavetable");
-  }
-  if (*space == NULL) {
-    rb_raise(rb_eRuntimeError, "something wrong with workspace");
-  }
-  return flag;
-}
-
-// Parse argc, argv.  obj must be GSL::Vector of halfcomplex data
-static int gsl_fft_get_argv_halfcomplex(int argc, VALUE *argv, VALUE obj,
-                                        double **ptr, size_t *stride,
-                                        size_t *n, gsl_fft_halfcomplex_wavetable **table,
-                                        gsl_fft_real_workspace **space, int *naflag)
-{
-  int flag = NONE_OF_TWO, flagtmp, i, itmp = argc, itmp2 = 0, ccc;
-  int flagw = 0;
-
-  *ptr = get_ptr_double3(obj, n, stride, naflag);
-
-  ccc = argc;
-  flagtmp = 0;
-  flagw = 0;
-  for (i = argc-1; i >= itmp2; i--) {
-    if (rb_obj_is_kind_of(argv[i], cgsl_fft_real_workspace)) {
-      Data_Get_Struct(argv[i], gsl_fft_real_workspace, *space);
-      flagtmp = 1;
-      flagw = 1;
-      itmp = i;
-      ccc--;
-      break;
-    }
-  }
-  flagtmp = 0;
-  for (i = itmp-1; i >= itmp2; i--) {
-    if (rb_obj_is_kind_of(argv[i], cgsl_fft_halfcomplex_wavetable)) {
-      Data_Get_Struct(argv[i], gsl_fft_halfcomplex_wavetable, *table);
-      flagtmp = 1;
-      ccc--;
-      break;
-    }
-  }
-  if (flagw == 0) {
-    *space = gsl_fft_real_workspace_alloc(*n);
-    flag += ALLOC_SPACE;
-  }
-  if (flagtmp == 0) {
-    *table = gsl_fft_halfcomplex_wavetable_alloc(*n);
-    flag += ALLOC_TABLE;
-  }
-  if (*table == NULL) {
-    rb_raise(rb_eRuntimeError, "something wrong with wavetable");
-  }
-  if (*space == NULL) {
-    rb_raise(rb_eRuntimeError, "something wrong with workspace");
-  }
-  return flag;
-}
-
-static void gsl_fft_free(int flag, GSL_FFT_Wavetable *table,
-                         GSL_FFT_Workspace *space)
-{
-  switch (flag) {
-  case ALLOC_TABLE:
-    GSL_FFT_Wavetable_free(table);
-    break;
-  case ALLOC_SPACE:
-    GSL_FFT_Workspace_free(space);
-    break;
-  case BOTH_OF_TWO:
-    GSL_FFT_Wavetable_free(table);
-    GSL_FFT_Workspace_free(space);
-    break;
-  default:
-    /* never happens */
-    break;
-  }
-}
-
-static VALUE rb_fft_complex_trans(int argc, VALUE *argv, VALUE obj,
-                                  int (*transform)(gsl_complex_packed_array,
-                                                   size_t, size_t,
-                                                   const gsl_fft_complex_wavetable *,
-                                                   gsl_fft_complex_workspace *),
-                                  int sss)
-{
-  int flag = 0;
-  // local variable "status" was defined and set, but never used
-  //int status;
-  size_t stride, n;
-  gsl_complex_packed_array data;
-  gsl_vector_complex *vin, *vout;
-  gsl_fft_complex_wavetable *table = NULL;
-  gsl_fft_complex_workspace *space = NULL;
-  flag = gsl_fft_get_argv_complex(argc, argv, obj, &vin, &data, &stride, &n, &table, &space);
-  if (sss == RB_GSL_FFT_COPY) {
-    vout = gsl_vector_complex_alloc(n);
-    gsl_vector_complex_memcpy(vout, vin);
-    /*status =*/ (*transform)(vout->data, vout->stride /*1*/, vout->size /*n*/, table, space);
-    gsl_fft_free(flag, (GSL_FFT_Wavetable *) table, (GSL_FFT_Workspace *) space);
-    return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vout);
-  } else {    /* in-place */
-    /*status =*/ (*transform)(data, stride, n, table, space);
-    gsl_fft_free(flag, (GSL_FFT_Wavetable *) table, (GSL_FFT_Workspace *) space);
-    return obj;
-  }
-}
-
-static VALUE rb_gsl_fft_complex_forward(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_complex_trans(argc, argv, obj, gsl_fft_complex_forward,
-                              RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_complex_forward2(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_complex_trans(argc, argv, obj, gsl_fft_complex_forward,
-                              RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_complex_transform(int argc, VALUE *argv, VALUE obj)
-{
-  int flag = 0;
-  // local variable "status" was defined and set, but never used
-  //int status;
-  size_t stride, n;
-  gsl_vector_complex *vin, *vout;
-  gsl_fft_direction sign;
-  gsl_complex_packed_array data;
-  gsl_fft_complex_wavetable *table = NULL;
-  gsl_fft_complex_workspace *space = NULL;
-  CHECK_FIXNUM(argv[argc-1]);
-  sign = FIX2INT(argv[argc-1]);
-  flag = gsl_fft_get_argv_complex(argc-1, argv, obj, &vin, &data, &stride, &n, &table, &space);
-  vout = gsl_vector_complex_alloc(n);
-  gsl_vector_complex_memcpy(vout, vin);
-  /*status =*/ gsl_fft_complex_transform(vout->data, stride, n, table, space, sign);
-  gsl_fft_free(flag, (GSL_FFT_Wavetable *) table, (GSL_FFT_Workspace *) space);
-  return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vout);
-}
-
-/* in-place */
-static VALUE rb_gsl_fft_complex_transform2(int argc, VALUE *argv, VALUE obj)
-{
-  int flag = 0;
-  // local variable "status" was defined and set, but never used
-  //int status;
-  size_t stride, n;
-  gsl_fft_direction sign;
-  gsl_complex_packed_array data;
-  gsl_fft_complex_wavetable *table = NULL;
-  gsl_fft_complex_workspace *space = NULL;
-  CHECK_FIXNUM(argv[argc-1]);
-  sign = FIX2INT(argv[argc-1]);
-  flag = gsl_fft_get_argv_complex(argc-1, argv, obj, NULL, &data, &stride, &n, &table, &space);
-  /*status =*/ gsl_fft_complex_transform(data, stride, n, table, space, sign);
-  gsl_fft_free(flag, (GSL_FFT_Wavetable *) table, (GSL_FFT_Workspace *) space);
-  return obj;
-}
-
-static VALUE rb_gsl_fft_complex_backward(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_complex_trans(argc, argv, obj, gsl_fft_complex_backward,
-                              RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_complex_backward2(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_complex_trans(argc, argv, obj, gsl_fft_complex_backward,
-                              RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_complex_inverse(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_complex_trans(argc, argv, obj, gsl_fft_complex_inverse,
-                              RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_complex_inverse2(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_complex_trans(argc, argv, obj, gsl_fft_complex_inverse,
-                              RB_GSL_FFT_INPLACE);
-}
-
-// The FFT methods used to allow passing stride and n values as optional
-// parameters to control which elements get transformed.  This created problems
-// for Views which can have their own stride, so support for stride and n
-// parameters to the transform methods is being dropped.  This method used to
-// be called to determine the stride and n values to use based on the
-// parameters and/or the vector itself (depending on how many parameters were
-// passed). Now this function is somewhat unneceesary, but to simplify the code
-// refactoring, it has been left in place for the time being.  Eventually it
-// can be refactored away completely.
-//
-// obj must be GSL::Vector of real or halfcomplex data
-static VALUE get_ptr_stride_n(VALUE obj,
-                              double **ptr, size_t *stride, size_t *n, int *flag)
-{
-  *flag = 0;
-  *ptr = get_ptr_double3(obj, n, stride, flag);
-  return obj;
-}
-
-static VALUE rb_fft_radix2(VALUE obj,
-                           int (*trans)(double [], size_t, size_t),
-                           int sss)
-{
-  size_t stride, n;
-  gsl_vector *vnew;
-  gsl_vector_view vv;
-  double *ptr1, *ptr2;
-  int flag;
-  VALUE ary;
-  get_ptr_stride_n(obj, &ptr1, &stride, &n, &flag);
-  if (flag == 0) {
-    if (sss == RB_GSL_FFT_COPY) {
-      vnew = gsl_vector_alloc(n);
-      vv.vector.data = ptr1;
-      vv.vector.stride = stride;
-      vv.vector.size = n;
-      gsl_vector_memcpy(vnew, &vv.vector);
-      ptr2 = vnew->data;
-      stride = 1;
-      ary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vnew);
-    } else {
-      ary = obj;
-      ptr2 = ptr1;
-    }
-#ifdef HAVE_NARRAY_H
-  } else if (flag == 1) {
-    if (sss == RB_GSL_FFT_COPY) {
-      int shape[1];
-      shape[0] = n;
-      ary = na_make_object(NA_DFLOAT, 1, shape, cNArray);
-      ptr2 = NA_PTR_TYPE(ary, double*);
-      memcpy(ptr2, ptr1, sizeof(double)*n);
-      stride = 1;
-    } else {
-      ary = obj;
-      ptr2 = NA_PTR_TYPE(ary, double*);
-    }
-#endif
-  } else {
-    rb_raise(rb_eRuntimeError, "something wrong");
-  }
-  (*trans)(ptr2, stride, n);
-  return ary;
-}
-
-static VALUE rb_gsl_fft_real_radix2_transform(VALUE obj)
-{
-  return rb_fft_radix2(obj, gsl_fft_real_radix2_transform,
-                       RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_real_radix2_transform2(VALUE obj)
-{
-  return rb_fft_radix2(obj, gsl_fft_real_radix2_transform,
-                       RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_radix2_inverse(VALUE obj)
-{
-  return rb_fft_radix2(obj, gsl_fft_halfcomplex_radix2_inverse,
-                       RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_radix2_inverse2(VALUE obj)
-{
-  return rb_fft_radix2(obj, gsl_fft_halfcomplex_radix2_inverse,
-                       RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_radix2_backward(VALUE obj)
-{
-  return rb_fft_radix2(obj, gsl_fft_halfcomplex_radix2_backward,
-                       RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_radix2_backward2(VALUE obj)
-{
-  return rb_fft_radix2(obj, gsl_fft_halfcomplex_radix2_backward,
-                       RB_GSL_FFT_INPLACE);
-}
-
-/*****/
-
-static VALUE rb_fft_real_trans(int argc, VALUE *argv, VALUE obj,
-                               int (*trans)(double [], size_t, size_t,
-                                            const gsl_fft_real_wavetable *,
-                                            gsl_fft_real_workspace *),
-                               int sss)
-{
-  int flag = 0, naflag = 0;
-  // local variable "status" was defined and set, but never used
-  //int status;
-  size_t stride, n;
-  gsl_vector *vnew;
-  gsl_vector_view vv;
-  double *ptr1, *ptr2;
-  gsl_fft_real_wavetable *table = NULL;
-  gsl_fft_real_workspace *space = NULL;
-  VALUE ary;
-  flag = gsl_fft_get_argv_real(argc, argv, obj, &ptr1, &stride, &n, &table, &space, &naflag);
-  if (naflag == 0) {
-    if (sss == RB_GSL_FFT_COPY) {
-      vnew = gsl_vector_alloc(n);
-      vv.vector.data = ptr1;
-      vv.vector.stride = stride;
-      vv.vector.size = n;
-      gsl_vector_memcpy(vnew, &vv.vector);
-      ptr2 = vnew->data;
-      stride = 1;
-      ary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vnew);
-    } else {
-      ptr2 = ptr1;
-      ary = obj;
-    }
-#ifdef HAVE_NARRAY_H
-  } else if (naflag == 1) {
-    if (sss == RB_GSL_FFT_COPY) {
-      int shape[1];
-      shape[0] = n;
-      ary = na_make_object(NA_DFLOAT, 1, shape, cNArray);
-      ptr2 = NA_PTR_TYPE(ary, double*);
-      memcpy(ptr2, ptr1, sizeof(double)*n);
-      stride = 1;
-    } else {
-      ptr2 = ptr1;
-      ary = obj;
-    }
-#endif
-  } else {
-    rb_raise(rb_eRuntimeError, "something wrong");
-  }
-  /*status =*/ (*trans)(ptr2, stride, n, table, space);
-  gsl_fft_free(flag, (GSL_FFT_Wavetable *) table, (GSL_FFT_Workspace *) space);
-  return ary;
-}
-
-static VALUE rb_gsl_fft_real_transform(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_real_trans(argc, argv, obj, gsl_fft_real_transform,
-                           RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_real_transform2(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_real_trans(argc, argv, obj, gsl_fft_real_transform,
-                           RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_fft_halfcomplex_trans(int argc, VALUE *argv, VALUE obj,
-                                      int (*trans)(double [], size_t, size_t,
-                                                   const gsl_fft_halfcomplex_wavetable *, gsl_fft_real_workspace *),
-                                      int sss)
-{
-  int flag = 0, naflag = 0;
-  // local variable "status" was defined and set, but never used
-  //int status;
-  size_t stride, n;
-  gsl_vector *vnew;
-  gsl_vector_view vv;
-  double *ptr1, *ptr2;
-  gsl_fft_halfcomplex_wavetable *table = NULL;
-  gsl_fft_real_workspace *space = NULL;
-  VALUE ary;
-  flag = gsl_fft_get_argv_halfcomplex(argc, argv, obj, &ptr1, &stride, &n,
-                                      &table, &space, &naflag);
-  if (naflag == 0) {
-    if (sss == RB_GSL_FFT_COPY) {
-      vnew = gsl_vector_alloc(n);
-      vv.vector.data = ptr1;
-      vv.vector.stride = stride;
-      vv.vector.size = n;
-      gsl_vector_memcpy(vnew, &vv.vector);
-      ptr2 = vnew->data;
-      stride = 1;
-      ary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vnew);
-    } else {
-      ptr2 = ptr1;
-      ary = obj;
-    }
-#ifdef HAVE_NARRAY_H
-  } else if (naflag == 1) {
-    if (sss == RB_GSL_FFT_COPY) {
-      int shape[1];
-      shape[0] = n;
-      ary = na_make_object(NA_DFLOAT, 1, shape, cNArray);
-      ptr2 = NA_PTR_TYPE(ary, double*);
-      memcpy(ptr2, ptr1, sizeof(double)*n);
-      stride = 1;
-    } else {
-      ptr2 = ptr1;
-      ary = obj;
-    }
-#endif
-  } else {
-    rb_raise(rb_eRuntimeError, "something wrong");
-  }
-  /*status =*/ (*trans)(ptr2, stride, n, table, space);
-  gsl_fft_free(flag, (GSL_FFT_Wavetable *) table, (GSL_FFT_Workspace *) space);
-  return ary;
-}
-
-static VALUE rb_gsl_fft_halfcomplex_transform(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_halfcomplex_trans(argc, argv, obj,
-                                  gsl_fft_halfcomplex_transform,
-                                  RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_transform2(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_halfcomplex_trans(argc, argv, obj,
-                                  gsl_fft_halfcomplex_transform,
-                                  RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_backward(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_halfcomplex_trans(argc, argv, obj,
-                                  gsl_fft_halfcomplex_backward,
-                                  RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_backward2(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_halfcomplex_trans(argc, argv, obj,
-                                  gsl_fft_halfcomplex_backward,
-                                  RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_inverse(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_halfcomplex_trans(argc, argv, obj,
-                                  gsl_fft_halfcomplex_inverse,
-                                  RB_GSL_FFT_COPY);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_inverse2(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_fft_halfcomplex_trans(argc, argv, obj,
-                                  gsl_fft_halfcomplex_inverse,
-                                  RB_GSL_FFT_INPLACE);
-}
-
-static VALUE rb_gsl_fft_real_unpack(VALUE obj)
-{
-  gsl_vector *v;
-  gsl_vector_complex *vout;
-
-  CHECK_VECTOR(obj);
-  Data_Get_Struct(obj, gsl_vector, v);
-
-  vout = gsl_vector_complex_alloc(v->size);
-  gsl_fft_real_unpack(v->data, (gsl_complex_packed_array) vout->data, v->stride, v->size);
-  return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vout);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_unpack(VALUE obj)
-{
-  gsl_vector *v;
-  gsl_vector_complex *vout;
-
-  CHECK_VECTOR(obj);
-  Data_Get_Struct(obj, gsl_vector, v);
-
-  vout = gsl_vector_complex_alloc(v->size);
-  gsl_fft_halfcomplex_unpack(v->data, (gsl_complex_packed_array) vout->data, v->stride, v->size);
-  return Data_Wrap_Struct(cgsl_vector_complex, 0, gsl_vector_complex_free, vout);
-}
-
-/* Convert a halfcomplex data to Numerical Recipes style */
-static VALUE rb_gsl_fft_halfcomplex_to_nrc(VALUE obj)
-{
-  gsl_vector *v, *vnew;
-  size_t i, k;
-
-  CHECK_VECTOR(obj);
-  Data_Get_Struct(obj, gsl_vector, v);
-
-  vnew = gsl_vector_alloc(v->size);
-  gsl_vector_set(vnew, 0, gsl_vector_get(v, 0));  /* DC */
-  gsl_vector_set(vnew, 1, gsl_vector_get(v, v->size/2));  /* Nyquist freq */
-  for (i = 2, k = 1; i < vnew->size; i += 2, k++) {
-    gsl_vector_set(vnew, i, gsl_vector_get(v, k));
-    gsl_vector_set(vnew, i+1, -gsl_vector_get(v, v->size-k));
-  }
-  return Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, vnew);
-}
-
-static VALUE rb_gsl_fft_halfcomplex_amp_phase(VALUE obj)
-{
-  gsl_vector *v;
-  gsl_vector *amp, *phase;
-  double re, im;
-  VALUE vamp, vphase;
-  size_t i;
-  CHECK_VECTOR(obj);
-  Data_Get_Struct(obj, gsl_vector, v);
-  amp = gsl_vector_alloc(v->size/2);
-  phase = gsl_vector_alloc(v->size/2);
-  gsl_vector_set(amp, 0, gsl_vector_get(v, 0));
-  gsl_vector_set(phase, 0, 0);
-  gsl_vector_set(amp, amp->size-1, gsl_vector_get(v, v->size-1));
-  gsl_vector_set(phase, phase->size-1, 0);
-  for (i = 1; i < v->size-1; i += 2) {
-    re = gsl_vector_get(v, i);
-    im = gsl_vector_get(v, i+1);
-    gsl_vector_set(amp, i/2+1, sqrt(re*re + im*im));
-    gsl_vector_set(phase, i/2+1, atan2(im, re));
-  }
-  vamp = Data_Wrap_Struct(VECTOR_ROW_COL(obj), 0, gsl_vector_free, amp);
-  vphase = Data_Wrap_Struct(VECTOR_ROW_COL(obj), 0, gsl_vector_free, phase);
-  return rb_ary_new3(2, vamp, vphase);
-}
-
-void Init_gsl_fft(VALUE module)
-{
-  mgsl_fft = rb_define_module_under(module, "FFT");
-
-  /*****/
-
-  rb_define_const(mgsl_fft, "Forward", INT2FIX(gsl_fft_forward));
-  rb_define_const(mgsl_fft, "FORWARD", INT2FIX(gsl_fft_forward));
-  rb_define_const(mgsl_fft, "Backward", INT2FIX(gsl_fft_backward));
-  rb_define_const(mgsl_fft, "BACKWARD", INT2FIX(gsl_fft_backward));
-
-  /* Transforms for complex vectors */
-  rb_define_method(cgsl_vector_complex, "radix2_forward",
-                   rb_gsl_fft_complex_radix2_forward, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_transform",
-                   rb_gsl_fft_complex_radix2_transform, 1);
-  rb_define_method(cgsl_vector_complex, "radix2_backward",
-                   rb_gsl_fft_complex_radix2_backward, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_inverse",
-                   rb_gsl_fft_complex_radix2_inverse, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_dif_forward",
-                   rb_gsl_fft_complex_radix2_dif_forward, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_dif_transform",
-                   rb_gsl_fft_complex_radix2_dif_transform, 1);
-  rb_define_method(cgsl_vector_complex, "radix2_dif_backward",
-                   rb_gsl_fft_complex_radix2_dif_backward, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_dif_inverse",
-                   rb_gsl_fft_complex_radix2_dif_inverse, 0);
-
-  /* In-place radix-2 transforms for complex vectors */
-  rb_define_method(cgsl_vector_complex, "radix2_forward!",
-                   rb_gsl_fft_complex_radix2_forward2, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_transform!",
-                   rb_gsl_fft_complex_radix2_transform2, 1);
-  rb_define_method(cgsl_vector_complex, "radix2_backward!",
-                   rb_gsl_fft_complex_radix2_backward2, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_inverse!",
-                   rb_gsl_fft_complex_radix2_inverse2, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_dif_forward!",
-                   rb_gsl_fft_complex_radix2_dif_forward2, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_dif_transform!",
-                   rb_gsl_fft_complex_radix2_dif_transform2, 1);
-  rb_define_method(cgsl_vector_complex, "radix2_dif_backward!",
-                   rb_gsl_fft_complex_radix2_dif_backward2, 0);
-  rb_define_method(cgsl_vector_complex, "radix2_dif_inverse!",
-                   rb_gsl_fft_complex_radix2_dif_inverse2, 0);
-
-  // class GSL::FFT::Wavetable < GSL::Object
-  //
-  // Useful since some functionality is shared among
-  // GSL::FFT::Complex::Wavetable
-  // GSL::FFT::Real::Wavetable
-  // GSL::FFT::HalfComplex::Wavetable
-  cgsl_fft_wavetable = rb_define_class_under(mgsl_fft, "Wavetable", cGSL_Object);
-  // No alloc
-  // TODO Make GSL::FFT::Wavetable#initialize private?
-  rb_define_method(cgsl_fft_wavetable, "n",
-                   rb_GSL_FFT_Wavetable_n, 0);
-  rb_define_method(cgsl_fft_wavetable, "nf",
-                   rb_GSL_FFT_Wavetable_nf, 0);
-  rb_define_method(cgsl_fft_wavetable, "factor",
-                   rb_GSL_FFT_Wavetable_factor, 0);
-
-  // class GSL::FFT::ComplexWavetable < GSL::FFT::Wavetable
-  cgsl_fft_complex_wavetable = rb_define_class_under(mgsl_fft, "ComplexWavetable",
-                                                     cgsl_fft_wavetable);
-  rb_define_singleton_method(cgsl_fft_complex_wavetable, "alloc",
-                             rb_gsl_fft_complex_wavetable_new, 1);
-
-  // class GSL::FFT::ComplexWorkspace < GSL::Object
-  cgsl_fft_complex_workspace = rb_define_class_under(mgsl_fft, "ComplexWorkspace",
-                                                     cGSL_Object);
-  rb_define_singleton_method(cgsl_fft_complex_workspace, "alloc",
-                             rb_gsl_fft_complex_workspace_new, 1);
-
-  rb_define_method(cgsl_vector_complex, "forward", rb_gsl_fft_complex_forward, -1);
-  rb_define_method(cgsl_vector_complex, "transform", rb_gsl_fft_complex_transform, -1);
-  rb_define_method(cgsl_vector_complex, "backward", rb_gsl_fft_complex_backward, -1);
-  rb_define_method(cgsl_vector_complex, "inverse", rb_gsl_fft_complex_inverse, -1);
-
-  rb_define_method(cgsl_vector_complex, "forward!", rb_gsl_fft_complex_forward2, -1);
-  rb_define_method(cgsl_vector_complex, "transform!", rb_gsl_fft_complex_transform2, -1);
-  rb_define_method(cgsl_vector_complex, "backward!", rb_gsl_fft_complex_backward2, -1);
-  rb_define_method(cgsl_vector_complex, "inverse!", rb_gsl_fft_complex_inverse2, -1);
-
-  /*****/
-
-  // TODO Do these method names need the "real_" and "halfcomplex_" prefixes?
-  rb_define_method(cgsl_vector, "real_radix2_transform",
-                   rb_gsl_fft_real_radix2_transform, 0);
-  rb_define_alias(cgsl_vector, "radix2_transform", "real_radix2_transform");
-  rb_define_alias(cgsl_vector, "radix2_forward", "real_radix2_transform");
-  rb_define_method(cgsl_vector, "real_radix2_inverse",
-                   rb_gsl_fft_halfcomplex_radix2_inverse, 0);
-  rb_define_alias(cgsl_vector, "radix2_inverse", "real_radix2_inverse");
-  rb_define_alias(cgsl_vector, "halfcomplex_radix2_inverse",
-                  "real_radix2_inverse");
-  rb_define_method(cgsl_vector, "real_radix2_backward",
-                   rb_gsl_fft_halfcomplex_radix2_backward, 0);
-  rb_define_alias(cgsl_vector, "radix2_backward", "real_radix2_backward");
-  rb_define_alias(cgsl_vector, "halfcomplex_radix2_backward",
-                  "real_radix2_backward");
-
-  // TODO Do these method names need the "real_" and "halfcomplex_" prefixes?
-  rb_define_method(cgsl_vector, "real_radix2_transform!",
-                   rb_gsl_fft_real_radix2_transform2, 0);
-  rb_define_alias(cgsl_vector, "radix2_transform!", "real_radix2_transform!");
-  rb_define_alias(cgsl_vector, "radix2_forward!", "real_radix2_transform!");
-  rb_define_method(cgsl_vector, "real_radix2_inverse!",
-                   rb_gsl_fft_halfcomplex_radix2_inverse2, 0);
-  rb_define_alias(cgsl_vector, "radix2_inverse!", "real_radix2_inverse!");
-  rb_define_alias(cgsl_vector, "halfcomplex_radix2_inverse!",
-                  "real_radix2_inverse!");
-  rb_define_method(cgsl_vector, "real_radix2_backward!",
-                   rb_gsl_fft_halfcomplex_radix2_backward2, 0);
-  rb_define_alias(cgsl_vector, "radix2_backward!", "real_radix2_backward!");
-  rb_define_alias(cgsl_vector, "halfcomplex_radix2_backward!",
-                  "real_radix2_backward!");
-
-  /*****/
-
-  // class GSL::FFT::RealWavetable < GSL::FFT::Wavetable
-  cgsl_fft_real_wavetable = rb_define_class_under(mgsl_fft, "RealWavetable",
-                                                  cgsl_fft_wavetable);
-  rb_define_singleton_method(cgsl_fft_real_wavetable, "alloc",
-                             rb_gsl_fft_real_wavetable_new, 1);
-
-  // class GSL::FFT::HalfComplexWavetable < GSL::FFT::Wavetable
-  cgsl_fft_halfcomplex_wavetable = rb_define_class_under(mgsl_fft,
-                                                         "HalfComplexWavetable", cgsl_fft_wavetable);
-  rb_define_singleton_method(cgsl_fft_halfcomplex_wavetable, "alloc",
-                             rb_gsl_fft_halfcomplex_wavetable_new, 1);
-
-  /*****/
-
-  // class GSL::FFT::RealWorkspace < GSL::Object
-  cgsl_fft_real_workspace = rb_define_class_under(mgsl_fft, "RealWorkspace",
-                                                  cGSL_Object);
-  rb_define_singleton_method(cgsl_fft_real_workspace, "alloc",
-                             rb_gsl_fft_real_workspace_new, 1);
-
-  /*****/
-
-  // TODO Do these method names need the "real_" and "halfcomplex_" prefixes?
-  rb_define_method(cgsl_vector, "real_transform", rb_gsl_fft_real_transform, -1);
-  rb_define_alias(cgsl_vector, "transform", "real_transform");
-  rb_define_alias(cgsl_vector, "forward", "real_transform");
-  rb_define_alias(cgsl_vector, "fft_forward", "real_transform");
-  rb_define_alias(cgsl_vector, "fft", "real_transform");
-  rb_define_method(cgsl_vector, "halfcomplex_transform",
-                   rb_gsl_fft_halfcomplex_transform, -1);
-  rb_define_method(cgsl_vector, "halfcomplex_backward",
-                   rb_gsl_fft_halfcomplex_backward, -1);
-  rb_define_alias(cgsl_vector, "backward", "halfcomplex_backward");
-  rb_define_alias(cgsl_vector, "fft_backward", "halfcomplex_backward");
-  rb_define_method(cgsl_vector, "halfcomplex_inverse",
-                   rb_gsl_fft_halfcomplex_inverse, -1);
-  rb_define_alias(cgsl_vector, "fft_inverse", "halfcomplex_inverse");
-  rb_define_alias(cgsl_vector, "ifft", "halfcomplex_inverse");
-  rb_define_alias(cgsl_vector, "inverse", "halfcomplex_inverse");
-
-  rb_define_method(cgsl_vector, "real_transform!", rb_gsl_fft_real_transform2, -1);
-  rb_define_alias(cgsl_vector, "transform!", "real_transform!");
-  rb_define_alias(cgsl_vector, "forward!", "real_transform!");
-  rb_define_alias(cgsl_vector, "fft_forward!", "real_transform!");
-  rb_define_alias(cgsl_vector, "fft!", "real_transform!");
-  rb_define_method(cgsl_vector, "halfcomplex_transform!",
-                   rb_gsl_fft_halfcomplex_transform2, -1);
-  rb_define_method(cgsl_vector, "halfcomplex_backward!",
-                   rb_gsl_fft_halfcomplex_backward2, -1);
-  rb_define_alias(cgsl_vector, "backward!", "halfcomplex_backward!");
-  rb_define_alias(cgsl_vector, "fft_backward!", "halfcomplex_backward!");
-  rb_define_method(cgsl_vector, "halfcomplex_inverse!",
-                   rb_gsl_fft_halfcomplex_inverse2, -1);
-  rb_define_alias(cgsl_vector, "fft_inverse!", "halfcomplex_inverse!");
-  rb_define_alias(cgsl_vector, "ifft!", "halfcomplex_inverse!");
-  rb_define_alias(cgsl_vector, "inverse!", "halfcomplex_inverse!");
-
-  /***/
-  rb_define_method(cgsl_vector, "fft_real_unpack", rb_gsl_fft_real_unpack, 0);
-  rb_define_alias(cgsl_vector, "real_unpack", "fft_real_unpack");
-  rb_define_alias(cgsl_vector, "real_to_complex", "fft_real_unpack");
-  rb_define_alias(cgsl_vector, "r_to_c", "fft_real_unpack");
-
-  rb_define_method(cgsl_vector, "fft_halfcomplex_unpack",
-                   rb_gsl_fft_halfcomplex_unpack, 0);
-  rb_define_alias(cgsl_vector, "halfcomplex_unpack", "fft_halfcomplex_unpack");
-  rb_define_alias(cgsl_vector, "halfcomplex_to_complex", "fft_halfcomplex_unpack");
-  rb_define_alias(cgsl_vector, "hc_to_c", "fft_halfcomplex_unpack");
-
-  /*****/
-
-  rb_define_method(cgsl_vector, "to_nrc_order",
-                   rb_gsl_fft_halfcomplex_to_nrc, 0);
-
-  rb_define_method(cgsl_vector, "halfcomplex_amp_phase",
-                   rb_gsl_fft_halfcomplex_amp_phase, 0);
-  rb_define_alias(cgsl_vector, "hc_amp_phase", "halfcomplex_amp_phase");
-}