1 | // FFT (Fast Fourier Transform) Server based on FFTPack
|
---|
2 | // R. Ansari 1999-2000
|
---|
3 | // DAPNIA/SPP (Saclay) / CEA LAL - IN2P3/CNRS (Orsay)
|
---|
4 |
|
---|
5 | #ifndef FFTPServer_H_SEEN
|
---|
6 | #define FFTPServer_H_SEEN
|
---|
7 |
|
---|
8 | #include "fftservintf.h"
|
---|
9 |
|
---|
10 | // implementation de FFTServerInterface en utilisant FFTPack
|
---|
11 |
|
---|
12 | namespace SOPHYA {
|
---|
13 |
|
---|
14 | //! An FFT server based on fftpack.
|
---|
15 | class FFTPackServer : public FFTServerInterface {
|
---|
16 | public:
|
---|
17 | FFTPackServer();
|
---|
18 | virtual ~FFTPackServer();
|
---|
19 |
|
---|
20 | // Implementation de l'interface FFTServerInterface
|
---|
21 |
|
---|
22 | virtual FFTServerInterface * Clone();
|
---|
23 |
|
---|
24 | // Transforme unidimensionnel sur des double
|
---|
25 | virtual void FFTForward(TArray< complex<r_8> > const & in, TArray< complex<r_8> > & out);
|
---|
26 | virtual void FFTBackward(TArray< complex<r_8> > const & in, TArray< complex<r_8> > & out);
|
---|
27 | virtual void FFTForward(TArray< r_8 > const & in, TArray< complex<r_8> > & out);
|
---|
28 | virtual void FFTBackward(TArray< complex<r_8> > const & in, TArray< r_8 > & out);
|
---|
29 |
|
---|
30 | // Transforme unidimensionnel sur des float
|
---|
31 | virtual void FFTForward(TArray< complex<r_4> > const & in, TArray< complex<r_4> > & out);
|
---|
32 | virtual void FFTBackward(TArray< complex<r_4> > const & in, TArray< complex<r_4> > & out);
|
---|
33 | virtual void FFTForward(TArray< r_4 > const & in, TArray< complex<r_4> > & out);
|
---|
34 | virtual void FFTBackward(TArray< complex<r_4> > const & in, TArray< r_4 > & out);
|
---|
35 |
|
---|
36 | // Methodes statiques pour reordonner les donnees
|
---|
37 | static void ReShapetoReal(TArray< complex<r_8> > const & in, TArray< r_8 > & out);
|
---|
38 | static void ReShapetoCompl(TArray< r_8 > const & in, TArray< complex<r_8> > & out);
|
---|
39 | static void ReShapetoReal(TArray< complex<r_4> > const & in, TArray< r_4 > & out);
|
---|
40 | static void ReShapetoCompl(TArray< r_4 > const & in, TArray< complex<r_4> > & out);
|
---|
41 |
|
---|
42 | // Methodes propres a cette classe
|
---|
43 | virtual void fftf(int_4 l, float* inout);
|
---|
44 | virtual void fftb(int_4 l, float* inout);
|
---|
45 | virtual void fftf(int_4 l, double* inout);
|
---|
46 | virtual void fftb(int_4 l, double* inout);
|
---|
47 | virtual void fftf(int_4 l, complex<float>* inout);
|
---|
48 | virtual void fftb(int_4 l, complex<float>* inout);
|
---|
49 | virtual void fftf(int_4 l, complex<double>* inout);
|
---|
50 | virtual void fftb(int_4 l, complex<double>* inout);
|
---|
51 |
|
---|
52 | protected:
|
---|
53 | virtual void checkint_rfft(int_4 l);
|
---|
54 | virtual void checkint_dfft(int_4 l);
|
---|
55 | virtual void checkint_cfft(int_4 l);
|
---|
56 | virtual void checkint_cdfft(int_4 l);
|
---|
57 |
|
---|
58 | int sz_rfft;
|
---|
59 | r_4* ws_rfft;
|
---|
60 |
|
---|
61 | int sz_cfft;
|
---|
62 | r_4* ws_cfft;
|
---|
63 |
|
---|
64 | int sz_dfft;
|
---|
65 | r_8* ws_dfft;
|
---|
66 |
|
---|
67 | int sz_cdfft;
|
---|
68 | r_8* ws_cdfft;
|
---|
69 |
|
---|
70 | FFTArrayChecker<r_4> ckR4;
|
---|
71 | FFTArrayChecker<r_8> ckR8;
|
---|
72 | };
|
---|
73 |
|
---|
74 | } // Fin du namespace
|
---|
75 |
|
---|
76 |
|
---|
77 | #endif
|
---|