FFTReal.h Example File

 /*****************************************************************************

         FFTReal.h
         Copyright (c) 2005 Laurent de Soras

 --- Legal stuff ---

 This library is free software; you can redistribute it and/or
 modify it under the terms of the GNU Lesser General Public
 License as published by the Free Software Foundation; either
 version 2.1 of the License, or (at your option) any later version.

 This library is distributed in the hope that it will be useful,
 but WITHOUT ANY WARRANTY; without even the implied warranty of
 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 Lesser General Public License for more details.

 You should have received a copy of the GNU Lesser General Public
 License along with this library; if not, write to the Free Software
 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

 *Tab=3***********************************************************************/

 #if ! defined (FFTReal_HEADER_INCLUDED)
 #define FFTReal_HEADER_INCLUDED

 #if defined (_MSC_VER)
         #pragma once
         #pragma warning (4 : 4250) // "Inherits via dominance."
 #endif

 /*\\\ INCLUDE FILES \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/

 #include        "def.h"
 #include        "DynArray.h"
 #include        "OscSinCos.h"

 template <class DT>
 class FFTReal
 {

 /*\\\ PUBLIC \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/

 public:

         enum {                  MAX_BIT_DEPTH   = 30    };      // So length can be represented as long int

         typedef DT      DataType;

         explicit                        FFTReal (long length);
         virtual                 ~FFTReal () {}

         long                            get_length () const;
         void                            do_fft (DataType f [], const DataType x []) const;
         void                            do_ifft (const DataType f [], DataType x []) const;
         void                            rescale (DataType x []) const;
         DataType *              use_buffer () const;

 /*\\\ PROTECTED \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/

 protected:

 /*\\\ PRIVATE \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/

 private:

    // Over this bit depth, we use direct calculation for sin/cos
    enum {             TRIGO_BD_LIMIT    = 12  };

         typedef OscSinCos <DataType>    OscType;

         void                            init_br_lut ();
         void                            init_trigo_lut ();
         void                            init_trigo_osc ();

         FORCEINLINE const long *
                                                 get_br_ptr () const;
         FORCEINLINE const DataType      *
                                                 get_trigo_ptr (int level) const;
         FORCEINLINE long
                                                 get_trigo_level_index (int level) const;

         inline void             compute_fft_general (DataType f [], const DataType x []) const;
         inline void             compute_direct_pass_1_2 (DataType df [], const DataType x []) const;
         inline void             compute_direct_pass_3 (DataType df [], const DataType sf []) const;
         inline void             compute_direct_pass_n (DataType df [], const DataType sf [], int pass) const;
         inline void             compute_direct_pass_n_lut (DataType df [], const DataType sf [], int pass) const;
         inline void             compute_direct_pass_n_osc (DataType df [], const DataType sf [], int pass) const;

         inline void             compute_ifft_general (const DataType f [], DataType x []) const;
         inline void             compute_inverse_pass_n (DataType df [], const DataType sf [], int pass) const;
         inline void             compute_inverse_pass_n_osc (DataType df [], const DataType sf [], int pass) const;
         inline void             compute_inverse_pass_n_lut (DataType df [], const DataType sf [], int pass) const;
         inline void             compute_inverse_pass_3 (DataType df [], const DataType sf []) const;
         inline void             compute_inverse_pass_1_2 (DataType x [], const DataType sf []) const;

         const long              _length;
         const int               _nbr_bits;
         DynArray <long>
                                                 _br_lut;
         DynArray <DataType>
                                                 _trigo_lut;
         mutable DynArray <DataType>
                                                 _buffer;
    mutable DynArray <OscType>
                                                 _trigo_osc;

 /*\\\ FORBIDDEN MEMBER FUNCTIONS \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/

 private:

                                                 FFTReal ();
                                                 FFTReal (const FFTReal &other);
         FFTReal &               operator = (const FFTReal &other);
         bool                            operator == (const FFTReal &other);
         bool                            operator != (const FFTReal &other);

 };      // class FFTReal

 #include        "FFTReal.hpp"

 #endif  // FFTReal_HEADER_INCLUDED

 /*\\\ EOF \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/