rb-gsl 1.16.0.5 → 1.16.0.6

data/ext/gsl_native/sf_trigonometric.c

@@ -1,210 +0,0 @@
-/*
-  sf_trigonometric.c
-  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
-    (C) Copyright 2001-2006 by Yoshiki Tsunesada
-
-  Ruby/GSL is free software: you can redistribute it and/or modify it
-  under the terms of the GNU General Public License.
-  This library is distributed in the hope that it will be useful, but
-  WITHOUT ANY WARRANTY.
-*/
-
-#include "include/rb_gsl_sf.h"
-
-static VALUE rb_gsl_sf_sin(VALUE obj, VALUE x)
-{
-  if (COMPLEX_P(x) || VECTOR_COMPLEX_P(x) || MATRIX_COMPLEX_P(x))
-    return rb_gsl_math_complex_eval(gsl_complex_sin, x);
-  return rb_gsl_sf_eval1(gsl_sf_sin, x);
-}
-
-static VALUE rb_gsl_sf_sin_e(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval_e(gsl_sf_sin_e, x);
-}
-
-static VALUE rb_gsl_sf_cos(VALUE obj, VALUE x)
-{
-  if (COMPLEX_P(x) || VECTOR_COMPLEX_P(x) || MATRIX_COMPLEX_P(x))
-    return rb_gsl_math_complex_eval(gsl_complex_cos, x);
-  return rb_gsl_sf_eval1(gsl_sf_cos, x);
-}
-
-static VALUE rb_gsl_sf_cos_e(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval_e(gsl_sf_cos_e, x);
-}
-
-static VALUE rb_gsl_sf_hypot(VALUE obj, VALUE x, VALUE y)
-{
-  Need_Float(x); Need_Float(y);
-  return rb_float_new(gsl_sf_hypot(NUM2DBL(x), NUM2DBL(y)));
-}
-
-static VALUE rb_gsl_sf_hypot_e(VALUE obj, VALUE x, VALUE y)
-{
-  return rb_gsl_sf_eval_e_double2(gsl_sf_hypot_e, x, y);
-}
-
-static VALUE rb_gsl_sf_sinc(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval1(gsl_sf_sinc, x);
-}
-
-static VALUE rb_gsl_sf_sinc_e(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval_e(gsl_sf_sinc_e, x);
-}
-
-static VALUE rb_gsl_sf_complex_XXX_e(int argc, VALUE *argv, VALUE obj,
-                                     int (*f)(double, double, gsl_sf_result*, gsl_sf_result*))
-{
-  gsl_sf_result *r1, *r2;
-  gsl_complex *z;
-  double re, im;
-  VALUE v1, v2;
-  // local variable "status" declared 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]);
-    break;
-  default:
-    rb_raise(rb_eArgError, "wrong number of arguments (%d for 1 or 2)", argc);
-    break;
-  }
-  v1 = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r1);
-  v2 = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, r2);
-  /*status =*/ (*f)(re, im, r1, r2);
-  return rb_ary_new3(2, v1, v2);
-}
-
-static VALUE rb_gsl_sf_complex_sin_e(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_sf_complex_XXX_e(argc, argv, obj, gsl_sf_complex_sin_e);
-}
-
-static VALUE rb_gsl_sf_complex_cos_e(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_sf_complex_XXX_e(argc, argv, obj, gsl_sf_complex_cos_e);
-}
-
-static VALUE rb_gsl_sf_complex_logsin_e(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_sf_complex_XXX_e(argc, argv, obj, gsl_sf_complex_logsin_e);
-}
-
-static VALUE rb_gsl_sf_lnsinh(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval1(gsl_sf_lnsinh, x);
-}
-
-static VALUE rb_gsl_sf_lnsinh_e(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval_e(gsl_sf_lnsinh_e, x);
-}
-
-static VALUE rb_gsl_sf_lncosh(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval1(gsl_sf_lncosh, x);
-}
-
-static VALUE rb_gsl_sf_lncosh_e(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval_e(gsl_sf_lncosh_e, x);
-}
-
-static VALUE rb_gsl_sf_polar_to_rect(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_sf_complex_XXX_e(argc, argv, obj, gsl_sf_polar_to_rect);
-}
-
-static VALUE rb_gsl_sf_rect_to_polar(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_sf_complex_XXX_e(argc, argv, obj, gsl_sf_rect_to_polar);
-}
-
-static VALUE rb_gsl_sf_angle_restrict_symm(VALUE obj, VALUE theta)
-{
-  return rb_gsl_sf_eval1(gsl_sf_angle_restrict_symm, theta);
-}
-
-static VALUE rb_gsl_sf_angle_restrict_pos(VALUE obj, VALUE theta)
-{
-  return rb_gsl_sf_eval1(gsl_sf_angle_restrict_pos, theta);
-}
-
-/*
-static VALUE rb_gsl_sf_sin_err(VALUE obj, VALUE x, VALUE dx)
-{
-  return rb_float_new(gsl_sf_sin_err(NUM2DBL(x), NUM2DBL(dx)));
-}
-
-static VALUE rb_gsl_sf_cos_err(VALUE obj, VALUE x, VALUE dx)
-{
-  return rb_float_new(gsl_sf_cos_err(NUM2DBL(x), NUM2DBL(dx)));
-}
-*/
-
-static VALUE rb_gsl_sf_sin_err_e(VALUE obj, VALUE x, VALUE dx)
-{
-  gsl_sf_result *rslt = NULL;
-  VALUE v;
-  // local variable "status" declared and set, but never used
-  //int status;
-  Need_Float(x);
-  Need_Float(dx);
-  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  /*status =*/ gsl_sf_sin_err_e(NUM2DBL(x), NUM2DBL(dx), rslt);
-  return v;
-}
-
-static VALUE rb_gsl_sf_cos_err_e(VALUE obj, VALUE x, VALUE dx)
-{
-  gsl_sf_result *rslt = NULL;
-  VALUE v;
-  // local variable "status" declared and set, but never used
-  //int status;
-  Need_Float(x);
-  Need_Float(dx);
-  v = Data_Make_Struct(cgsl_sf_result, gsl_sf_result, 0, free, rslt);
-  /*status =*/ gsl_sf_cos_err_e(NUM2DBL(x), NUM2DBL(dx), rslt);
-  return v;
-}
-
-void Init_gsl_sf_trigonometric(VALUE module)
-{
-  rb_define_module_function(module, "sin",  rb_gsl_sf_sin, 1);
-  rb_define_module_function(module, "sin_e",  rb_gsl_sf_sin_e, 1);
-  rb_define_module_function(module, "cos",  rb_gsl_sf_cos, 1);
-  rb_define_module_function(module, "cos_e",  rb_gsl_sf_cos_e, 1);
-  rb_define_module_function(module, "hypot",  rb_gsl_sf_hypot, 2);
-  rb_define_module_function(module, "hypot_e",  rb_gsl_sf_hypot_e, 2);
-  rb_define_module_function(module, "sinc",  rb_gsl_sf_sinc, 1);
-  rb_define_module_function(module, "sinc_e",  rb_gsl_sf_sinc_e, 1);
-  rb_define_module_function(module, "complex_sin_e",  rb_gsl_sf_complex_sin_e, -1);
-  rb_define_module_function(module, "complex_cos_e",  rb_gsl_sf_complex_cos_e, -1);
-  rb_define_module_function(module, "complex_logsin_e",  rb_gsl_sf_complex_logsin_e, -1);
-  rb_define_module_function(module, "lnsinh",  rb_gsl_sf_lnsinh, 1);
-  rb_define_module_function(module, "lnsinh_e",  rb_gsl_sf_lnsinh_e, 1);
-  rb_define_module_function(module, "lncosh",  rb_gsl_sf_lncosh, 1);
-  rb_define_module_function(module, "lncosh_e",  rb_gsl_sf_lncosh_e, 1);
-  rb_define_module_function(module, "polar_to_rect",  rb_gsl_sf_polar_to_rect, -1);
-  rb_define_module_function(module, "rect_to_polar",  rb_gsl_sf_rect_to_polar, -1);
-  rb_define_module_function(module, "angle_restrict_symm",  rb_gsl_sf_angle_restrict_symm, 1);
-  rb_define_module_function(module, "angle_restrict_pos",  rb_gsl_sf_angle_restrict_pos, 1);
-
-  /*  rb_define_module_function(module, "sin_err",  rb_gsl_sf_sin_err, 2);
-      rb_define_module_function(module, "cos_err",  rb_gsl_sf_cos_err, 2);*/
-  rb_define_module_function(module, "sin_err_e",  rb_gsl_sf_sin_err_e, 2);
-  rb_define_module_function(module, "cos_err_e",  rb_gsl_sf_cos_err_e, 2);
-
-}

data/ext/gsl_native/sf_zeta.c

@@ -1,115 +0,0 @@
-/*
-  sf_zeta.c
-  Ruby/GSL: Ruby extension library for GSL (GNU Scientific Library)
-    (C) Copyright 2001-2006 by Yoshiki Tsunesada
-
-  Ruby/GSL is free software: you can redistribute it and/or modify it
-  under the terms of the GNU General Public License.
-  This library is distributed in the hope that it will be useful, but
-  WITHOUT ANY WARRANTY.
-*/
-
-#include "include/rb_gsl_sf.h"
-
-static VALUE rb_gsl_sf_zeta_int(VALUE obj, VALUE n)
-{
-  VALUE nn;
-  if (TYPE(n) != T_FIXNUM) nn = INT2FIX(NUM2INT(n));
-  else nn = n;
-  return rb_gsl_sf_eval1_int(gsl_sf_zeta_int, nn);
-}
-
-static VALUE rb_gsl_sf_zeta_int_e(VALUE obj, VALUE n)
-{
-  VALUE nn;
-  if (TYPE(n) != T_FIXNUM) nn = INT2FIX(NUM2INT(n));
-  else nn = n;
-  return rb_gsl_sf_eval_e_int(gsl_sf_zeta_int_e, nn);
-}
-
-static VALUE rb_gsl_sf_zeta(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval1(gsl_sf_zeta, x);
-}
-
-static VALUE rb_gsl_sf_zeta_e(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval_e(gsl_sf_zeta_e, x);
-}
-
-static VALUE rb_gsl_sf_hzeta(VALUE obj, VALUE s, VALUE q)
-{
-  Need_Float(s); Need_Float(q);
-  return rb_float_new(gsl_sf_hzeta(NUM2DBL(s), NUM2DBL(q)));
-}
-
-static VALUE rb_gsl_sf_hzeta_e(VALUE obj, VALUE s, VALUE q)
-{
-  return rb_gsl_sf_eval_e_double2(gsl_sf_hzeta_e, s, q);
-}
-
-static VALUE rb_gsl_sf_eta_int(VALUE obj, VALUE n)
-{
-  return rb_gsl_sf_eval1_int(gsl_sf_eta_int, n);
-}
-
-static VALUE rb_gsl_sf_eta_int_e(VALUE obj, VALUE n)
-{
-  return rb_gsl_sf_eval_e_int(gsl_sf_eta_int_e, n);
-}
-
-static VALUE rb_gsl_sf_eta(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval1(gsl_sf_eta, x);
-}
-
-static VALUE rb_gsl_sf_eta_e(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval_e(gsl_sf_eta_e, x);
-}
-
-static VALUE rb_gsl_sf_zetam1_int(VALUE obj, VALUE n)
-{
-  VALUE nn;
-  if (TYPE(n) != T_FIXNUM) nn = INT2FIX(NUM2INT(n));
-  else nn = n;
-  return rb_gsl_sf_eval1_int(gsl_sf_zetam1_int, nn);
-}
-
-static VALUE rb_gsl_sf_zetam1_int_e(VALUE obj, VALUE n)
-{
-  VALUE nn;
-  if (TYPE(n) != T_FIXNUM) nn = INT2FIX(NUM2INT(n));
-  else nn = n;
-  return rb_gsl_sf_eval_e_int(gsl_sf_zetam1_int_e, nn);
-}
-
-static VALUE rb_gsl_sf_zetam1(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval1(gsl_sf_zetam1, x);
-}
-
-static VALUE rb_gsl_sf_zetam1_e(VALUE obj, VALUE x)
-{
-  return rb_gsl_sf_eval_e(gsl_sf_zetam1_e, x);
-}
-
-void Init_gsl_sf_zeta(VALUE module)
-{
-  rb_define_module_function(module, "zeta_int",  rb_gsl_sf_zeta_int, 1);
-  rb_define_module_function(module, "zeta_int_e",  rb_gsl_sf_zeta_int_e, 1);
-  rb_define_module_function(module, "zeta",  rb_gsl_sf_zeta, 1);
-  rb_define_module_function(module, "zeta_e",  rb_gsl_sf_zeta_e, 1);
-
-  rb_define_module_function(module, "hzeta",  rb_gsl_sf_hzeta, 2);
-  rb_define_module_function(module, "hzeta_e",  rb_gsl_sf_hzeta_e, 2);
-  rb_define_module_function(module, "eta_int",  rb_gsl_sf_eta_int, 1);
-  rb_define_module_function(module, "eta_int_e",  rb_gsl_sf_eta_int_e, 1);
-  rb_define_module_function(module, "eta",  rb_gsl_sf_eta, 1);
-  rb_define_module_function(module, "eta_e",  rb_gsl_sf_eta_e, 1);
-
-  rb_define_module_function(module, "zetam1_int",  rb_gsl_sf_zetam1_int, 1);
-  rb_define_module_function(module, "zetam1_int_e",  rb_gsl_sf_zetam1_int_e, 1);
-  rb_define_module_function(module, "zetam1",  rb_gsl_sf_zetam1, 1);
-  rb_define_module_function(module, "zetam1_e",  rb_gsl_sf_zetam1_e, 1);
-}

data/ext/gsl_native/signal.c

@@ -1,303 +0,0 @@
-/*
-  signal.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"
-
-enum FFT_CONV_CORR {
-  RB_GSL_FFT_CONVOLVE = 0,
-  RB_GSL_FFT_CORRELATE = 1,
-  RB_GSL_FFT_REAL = 2,
-  RB_GSL_FFT_HALFCOMPLEX = 3,
-  RB_GSL_FFT_DECONVOLVE = 4,
-};
-
-#ifndef WAVETABLE_P
-#define WAVETABLE_P(x) (rb_obj_is_kind_of(x,cgsl_fft_halfcomplex_wavetable) ? 1 : 0)
-#endif
-
-#ifndef CHECK_WAVETABLE
-#define CHECK_WAVETABLE(x) if(!rb_obj_is_kind_of(x,cgsl_fft_halfcomplex_wavetable)) \
-    rb_raise(rb_eTypeError, "wrong argument type (FFT::HalfComplex::Wavetable expected)");
-#endif
-
-#ifndef WORKSPACE_P
-#define WORKSPACE_P(x) (rb_obj_is_kind_of(x,cgsl_fft_real_workspace) ? 1 : 0)
-#endif
-
-#ifndef CHECK_WORKSPACE
-#define CHECK_WORKSPACE(x) if(!rb_obj_is_kind_of(x,cgsl_fft_real_workspace)) \
-    rb_raise(rb_eTypeError, "wrong argument type (FFT::Real::Workspace expected)");
-#endif
-
-static void complex_mul(double re1, double im1, double re2, double im2,
-                        double *re, double *im)
-{
-  *re = re1*re2 - im1*im2;
-  *im = re1*im2 + im1*re2;
-}
-
-static void complex_conj_mul(double re1, double im1, double re2, double im2,
-                             double *re, double *im)
-{
-  *re = re1*re2 + im1*im2;
-  *im = -re1*im2 + im1*re2;
-}
-
-static void complex_div(double re1, double im1, double re2, double im2,
-                        double *re, double *im)
-{
-  double factor = re2*re2 + im2*im2;
-  complex_conj_mul(re1, im1, re2, im2, re, im);
-  *re /= factor;
-  *im /= factor;
-}
-
-/* data1, data2: FFTed data */
-static void rbgsl_calc_conv_corr_c(const double *data1, const double *data2,
-                                   double *data3, size_t size,
-                                   gsl_fft_halfcomplex_wavetable *table,
-                                   gsl_fft_real_workspace *space, enum FFT_CONV_CORR calcflag)
-{
-  size_t i;
-  double re1, re2, im1, im2;
-  void (*complex_cal)(double, double, double, double, double*, double*);
-
-  switch (calcflag) {
-  case RB_GSL_FFT_CONVOLVE:
-    complex_cal = complex_mul;
-    data3[0] = data1[0]*data2[0];
-    data3[size-1] = data1[size-1]*data2[size-1];
-    break;
-  case RB_GSL_FFT_CORRELATE:
-    data3[0] = data1[0]*data2[0];
-    data3[size-1] = data1[size-1]*data2[size-1];
-    complex_cal = complex_conj_mul;
-    break;
-  case RB_GSL_FFT_DECONVOLVE:
-    complex_cal = complex_div;
-    data3[0] = data1[0]/data2[0];
-    data3[size-1] = data1[size-1]/data2[size-1];
-    break;
-  default:
-    rb_raise(rb_eArgError, "Wrong flag.");
-    break;
-  }
-
-  for (i = 1; i < size-1; i += 2) {
-    re1 = data1[i];    im1 = data1[i+1];
-    re2 = data2[i];    im2 = data2[i+1];
-    (*complex_cal)(re1, im1, re2, im2, &data3[i], &data3[i+1]);
-  }
-}
-
-static VALUE rb_gsl_fft_conv_corr(int argc, VALUE *argv, VALUE obj,
-                                  enum FFT_CONV_CORR flag1,
-                                  enum FFT_CONV_CORR flag2)
-{
-  double *data1, *data2, *data3 = NULL;
-  size_t stride1, stride2, size1, size2;
-  int naflag1, naflag2;
-  gsl_vector *v = NULL;
-  gsl_fft_halfcomplex_wavetable *table = NULL;
-  gsl_fft_real_wavetable *rtable = NULL;
-  gsl_fft_real_workspace *space = NULL, *space2 = NULL;
-  int flagt = 0, flagw = 0;
-  //  size_t i;
-  gsl_vector *vtmp1 = NULL, *vtmp2 = NULL;
-  VALUE ary = Qnil;
-  switch (argc) {
-  case 3:
-    data1 = get_ptr_double3(obj, &size1, &stride1, &naflag1);
-    data2 = get_ptr_double3(argv[0], &size2, &stride2, &naflag2);
-    CHECK_WAVETABLE(argv[1]);
-    Data_Get_Struct(argv[1], gsl_fft_halfcomplex_wavetable, table);
-    CHECK_WORKSPACE(argv[2]);
-    Data_Get_Struct(argv[2], gsl_fft_real_workspace, space);
-    break;
-  case 2:
-    data1 = get_ptr_double3(obj, &size1, &stride1, &naflag1);
-    data2 = get_ptr_double3(argv[0], &size2, &stride2, &naflag2);
-    if (WAVETABLE_P(argv[1])) {
-      Data_Get_Struct(argv[1], gsl_fft_halfcomplex_wavetable, table);
-      space = gsl_fft_real_workspace_alloc(size1);
-      flagw = 1;
-    } else if (WORKSPACE_P(argv[1])) {
-      Data_Get_Struct(argv[1], gsl_fft_real_workspace, space);
-      table = gsl_fft_halfcomplex_wavetable_alloc(size1);
-      flagt = 1;
-    } else {
-      rb_raise(rb_eTypeError,
-               "wrong argument type %s "
-               "(FFT::HalfComplex::Wavetable or FFT::Real::Workspace expected)",
-               rb_class2name(CLASS_OF(argv[2])));
-    }
-    break;
-  case 1:
-    data1 = get_ptr_double3(obj, &size1, &stride1, &naflag1);
-    data2 = get_ptr_double3(argv[0], &size2, &stride2, &naflag2);
-    table = gsl_fft_halfcomplex_wavetable_alloc(size1);
-    space = gsl_fft_real_workspace_alloc(size1);
-    flagt = 1;
-    flagw = 1;
-    break;
-  default:
-    rb_raise(rb_eArgError, "wrong number of arguments (%d for 1-3)", argc);
-  }
-
-  switch (naflag1*naflag2) {
-  case 0:
-    v = gsl_vector_alloc(size1);
-    switch (flag1) {
-    case RB_GSL_FFT_REAL:
-      ary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
-      break;
-    default:
-      ary = Data_Wrap_Struct(cgsl_vector, 0, gsl_vector_free, v);
-      break;
-    }
-    data3 = v->data;
-    break;
-  case 1:
-#ifdef HAVE_NARRAY_H
-    {
-      int shape = (int) size1;
-      ary = na_make_object(NA_DFLOAT, 1, &shape, cNArray);
-      data3 = NA_PTR_TYPE(ary, double*);
-    }
-#endif
-    break;
-  default:
-    break;
-  }
-  switch (flag1) {
-  case RB_GSL_FFT_REAL:  // do FFT
-    vtmp1 = gsl_vector_alloc(size1);
-    vtmp2 = gsl_vector_alloc(size2);
-    memcpy(vtmp1->data, data1, sizeof(double)*size1);
-    memcpy(vtmp2->data, data2, sizeof(double)*size2);
-    data1 = vtmp1->data;
-    data2 = vtmp2->data;
-    rtable = gsl_fft_real_wavetable_alloc(size1);
-    if (size1 == space->n) {
-      gsl_fft_real_transform(data1, stride1, size1, rtable, space);
-    } else {
-      space2 = gsl_fft_real_workspace_alloc(size1);
-      gsl_fft_real_transform(data1, stride1, size1, rtable, space2);
-      /* no freeing space2 here */
-    }
-    if (size1 != size2) {
-      if (rtable) gsl_fft_real_wavetable_free(rtable);
-      rtable = gsl_fft_real_wavetable_alloc(size2);
-    }
-    if (size2 == space->n) {
-      gsl_fft_real_transform(data2, stride2, size2, rtable, space);
-    } else if (size2 == size1) {
-      gsl_fft_real_transform(data2, stride2, size2, rtable, space2);
-      gsl_fft_real_workspace_free(space2);
-    } else {
-      if (space2) gsl_fft_real_workspace_free(space2);
-      space2 = gsl_fft_real_workspace_alloc(size2);
-      gsl_fft_real_transform(data2, stride2, size2, rtable, space2);
-      gsl_fft_real_workspace_free(space2);
-    }
-    gsl_fft_real_wavetable_free(rtable);
-    space2 = NULL;
-    rtable = NULL;
-    break;
-  case RB_GSL_FFT_HALFCOMPLEX:
-    /* do nothing */
-    break;
-  default:
-    /* not occur */
-    break;
-  }
-
-  rbgsl_calc_conv_corr_c(data1, data2, data3, size1, table, space, flag2);
-
-  if (flag1 == RB_GSL_FFT_REAL) {
-    gsl_fft_halfcomplex_inverse(data3, 1, size1, table, space);
-//    for (i = 0; i < size1; i++) data3[i] /= size1;
-  }
-  if (flagt == 1) gsl_fft_halfcomplex_wavetable_free(table);
-  if (flagw == 1) gsl_fft_real_workspace_free(space);
-  if (vtmp1) gsl_vector_free(vtmp1);
-  if (vtmp2) gsl_vector_free(vtmp2);
-  return ary;
-}
-
-/* GSL::Vector#convolve */
-static VALUE rb_gsl_fft_real_convolve(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_fft_conv_corr(argc, argv, obj,
-                              RB_GSL_FFT_REAL,
-                              RB_GSL_FFT_CONVOLVE);
-}
-/* GSL::Vector#deconvolve */
-static VALUE rb_gsl_fft_real_deconvolve(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_fft_conv_corr(argc, argv, obj,
-                              RB_GSL_FFT_REAL,
-                              RB_GSL_FFT_DECONVOLVE);
-}
-
-/* GSL::Vector#correlate */
-static VALUE rb_gsl_fft_real_correlate(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_fft_conv_corr(argc, argv, obj,
-                              RB_GSL_FFT_REAL,
-                              RB_GSL_FFT_CORRELATE);
-}
-
-/* GSL::Vector#halfcomplex_convolve */
-static VALUE rb_gsl_fft_halfcomplex_convolve(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_fft_conv_corr(argc, argv, obj,
-                              RB_GSL_FFT_HALFCOMPLEX,
-                              RB_GSL_FFT_CONVOLVE);
-}
-/* GSL::Vector#halfcomplex_deconvolve */
-static VALUE rb_gsl_fft_halfcomplex_deconvolve(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_fft_conv_corr(argc, argv, obj,
-                              RB_GSL_FFT_HALFCOMPLEX,
-                              RB_GSL_FFT_DECONVOLVE);
-}
-/* GSL::Vector#halfcomplex_correlate */
-static VALUE rb_gsl_fft_halfcomplex_correlate(int argc, VALUE *argv, VALUE obj)
-{
-  return rb_gsl_fft_conv_corr(argc, argv, obj,
-                              RB_GSL_FFT_HALFCOMPLEX,
-                              RB_GSL_FFT_CORRELATE);
-}
-
-void Init_gsl_signal(VALUE module)
-{
-  rb_define_method(cgsl_vector, "real_convolve", rb_gsl_fft_real_convolve, -1);
-  rb_define_method(cgsl_vector, "real_deconvolve", rb_gsl_fft_real_deconvolve, -1);
-  rb_define_method(cgsl_vector, "real_correlate", rb_gsl_fft_real_correlate, -1);
-
-  rb_define_alias(cgsl_vector, "convolve", "real_convolve");
-  rb_define_alias(cgsl_vector, "deconvolve", "real_deconvolve");
-  rb_define_alias(cgsl_vector, "correlate", "real_correlate");
-
-  rb_define_method(cgsl_vector, "halfcomplex_convolve", rb_gsl_fft_halfcomplex_convolve, -1);
-  rb_define_method(cgsl_vector, "halfcomplex_deconvolve", rb_gsl_fft_halfcomplex_deconvolve, -1);
-  rb_define_method(cgsl_vector, "halfcomplex_correlate", rb_gsl_fft_halfcomplex_correlate, -1);
-
-  rb_define_alias(cgsl_vector, "hc_convolve", "halfcomplex_convolve");
-  rb_define_alias(cgsl_vector, "hc_deconvolve", "halfcomplex_deconvolve");
-  rb_define_alias(cgsl_vector, "hc_correlate", "halfcomplex_correlate");
-}
-
-#undef WAVETABLE_P
-#undef CHECK_WAVETABLE
-#undef WORKSPACE_P
-#undef CHECK_WORKSPACE