rb-gsl 1.15.3.1

data/ext/sf_trigonometric.c

@@ -0,0 +1,210 @@
+/*
+  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 "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/sf_zeta.c

@@ -0,0 +1,119 @@
+/*
+  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 "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);
+}
+
+#ifdef GSL_1_4_9_LATER
+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);
+}
+#endif
+
+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);
+
+#ifdef GSL_1_4_9_LATER
+  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);
+#endif
+}

data/ext/signal.c

@@ -0,0 +1,308 @@
+/*
+  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 "rb_gsl_config.h"
+#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;
+  size_t stride1, stride2, size1, size2;
+#ifdef HAVE_NARRAY_H
+  int naflag1, naflag2, shape;
+#else
+  int naflag1, naflag2;
+#endif
+  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;
+  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
+    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