1 /* fft/gsl_fft_complex_float.h
3 * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2007 Brian Gough
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.
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.
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.
20 #ifndef __GSL_FFT_COMPLEX_FLOAT_H__
21 #define __GSL_FFT_COMPLEX_FLOAT_H__
25 #include <gsl/gsl_math.h>
26 #include <gsl/gsl_complex.h>
27 #include <gsl/gsl_fft.h>
32 # define __BEGIN_DECLS extern "C" {
33 # define __END_DECLS }
35 # define __BEGIN_DECLS /* empty */
36 # define __END_DECLS /* empty */
41 /* Power of 2 routines */
44 int gsl_fft_complex_float_radix2_forward (gsl_complex_packed_array_float data,
48 int gsl_fft_complex_float_radix2_backward (gsl_complex_packed_array_float data,
52 int gsl_fft_complex_float_radix2_inverse (gsl_complex_packed_array_float data,
56 int gsl_fft_complex_float_radix2_transform (gsl_complex_packed_array_float data,
59 const gsl_fft_direction sign);
61 int gsl_fft_complex_float_radix2_dif_forward (gsl_complex_packed_array_float data,
65 int gsl_fft_complex_float_radix2_dif_backward (gsl_complex_packed_array_float data,
69 int gsl_fft_complex_float_radix2_dif_inverse (gsl_complex_packed_array_float data,
73 int gsl_fft_complex_float_radix2_dif_transform (gsl_complex_packed_array_float data,
76 const gsl_fft_direction sign);
78 /* Mixed Radix general-N routines */
85 gsl_complex_float *twiddle[64];
86 gsl_complex_float *trig;
88 gsl_fft_complex_wavetable_float;
95 gsl_fft_complex_workspace_float;
98 gsl_fft_complex_wavetable_float *gsl_fft_complex_wavetable_float_alloc (size_t n);
100 void gsl_fft_complex_wavetable_float_free (gsl_fft_complex_wavetable_float * wavetable);
102 gsl_fft_complex_workspace_float *gsl_fft_complex_workspace_float_alloc (size_t n);
104 void gsl_fft_complex_workspace_float_free (gsl_fft_complex_workspace_float * workspace);
107 int gsl_fft_complex_float_memcpy (gsl_fft_complex_wavetable_float * dest,
108 gsl_fft_complex_wavetable_float * src);
111 int gsl_fft_complex_float_forward (gsl_complex_packed_array_float data,
114 const gsl_fft_complex_wavetable_float * wavetable,
115 gsl_fft_complex_workspace_float * work);
117 int gsl_fft_complex_float_backward (gsl_complex_packed_array_float data,
120 const gsl_fft_complex_wavetable_float * wavetable,
121 gsl_fft_complex_workspace_float * work);
123 int gsl_fft_complex_float_inverse (gsl_complex_packed_array_float data,
126 const gsl_fft_complex_wavetable_float * wavetable,
127 gsl_fft_complex_workspace_float * work);
129 int gsl_fft_complex_float_transform (gsl_complex_packed_array_float data,
130 const size_t stride, const size_t n,
131 const gsl_fft_complex_wavetable_float * wavetable,
132 gsl_fft_complex_workspace_float * work,
133 const gsl_fft_direction sign);
137 #endif /* __GSL_FFT_COMPLEX_FLOAT_H__ */