gsl_fft_real.h 2.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081
  1. /* fft/gsl_fft_real.h
  2. *
  3. * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 Brian Gough
  4. *
  5. * This program is free software; you can redistribute it and/or modify
  6. * it under the terms of the GNU General Public License as published by
  7. * the Free Software Foundation; either version 3 of the License, or (at
  8. * your option) any later version.
  9. *
  10. * This program is distributed in the hope that it will be useful, but
  11. * WITHOUT ANY WARRANTY; without even the implied warranty of
  12. * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  13. * General Public License for more details.
  14. *
  15. * You should have received a copy of the GNU General Public License
  16. * along with this program; if not, write to the Free Software
  17. * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  18. */
  19. #ifndef __GSL_FFT_REAL_H__
  20. #define __GSL_FFT_REAL_H__
  21. #include <stddef.h>
  22. #include "gsl_math.h"
  23. #include "gsl_complex.h"
  24. #include "gsl_fft.h"
  25. #undef __BEGIN_DECLS
  26. #undef __END_DECLS
  27. #ifdef __cplusplus
  28. # define __BEGIN_DECLS extern "C" {
  29. # define __END_DECLS }
  30. #else
  31. # define __BEGIN_DECLS /* empty */
  32. # define __END_DECLS /* empty */
  33. #endif
  34. __BEGIN_DECLS
  35. int gsl_fft_real_radix2_transform (double data[], const size_t stride, const size_t n) ;
  36. typedef struct
  37. {
  38. size_t n;
  39. size_t nf;
  40. size_t factor[64];
  41. gsl_complex *twiddle[64];
  42. gsl_complex *trig;
  43. }
  44. gsl_fft_real_wavetable;
  45. typedef struct
  46. {
  47. size_t n;
  48. double *scratch;
  49. }
  50. gsl_fft_real_workspace;
  51. gsl_fft_real_wavetable * gsl_fft_real_wavetable_alloc (size_t n);
  52. void gsl_fft_real_wavetable_free (gsl_fft_real_wavetable * wavetable);
  53. gsl_fft_real_workspace * gsl_fft_real_workspace_alloc (size_t n);
  54. void gsl_fft_real_workspace_free (gsl_fft_real_workspace * workspace);
  55. int gsl_fft_real_transform (double data[], const size_t stride, const size_t n,
  56. const gsl_fft_real_wavetable * wavetable,
  57. gsl_fft_real_workspace * work);
  58. int gsl_fft_real_unpack (const double real_coefficient[],
  59. double complex_coefficient[],
  60. const size_t stride, const size_t n);
  61. __END_DECLS
  62. #endif /* __GSL_FFT_REAL_H__ */