UT/UT_FFT.h

Go to the documentation of this file.

/*
 * PROPRIETARY INFORMATION.  This software is proprietary to
 * Side Effects Software Inc., and is not to be reproduced,
 * transmitted, or disclosed in any way without written permission.
 *
 * Produced by:
 *      Mark Alexander
 *      Side Effects Software Inc
 *      477 Richmond Street West
 *      Toronto, Ontario
 *      Canada   M5V 3E7
 *      416-504-9876
 *
 * NAME:        UT_FFT.h
 *
 * COMMENTS:
 *
 *      Performs a fast discrete fourier transform (or an inverse) on
 *      array data.
 *      FORWARD: Takes real data and produces complex data.
 *      INVERSE: Takes complex data and produces real data.
 */

#ifndef __UT_FFT__
#define __UT_FFT__

#include "UT_API.h"
#include <SYS/SYS_Types.h>

class UT_API UT_FFT
{
public:
    explicit UT_FFT(exint npnts);
    ~UT_FFT();


    void toFrequency(float *source,float *destr,float *desti,exint npnts);
    void toFrequency(float *sourcer, float *sourcei,float *destr,float *desti,exint npnts);
    void toTime(float *sourcer,float *sourcei,float *dest,exint npnts);

    void toTime(float *sourcer,float *sourcei,float *destr,float *desti,exint npnts);

    // two dimensional FFT'ing
    static void toFrequency2D(float *sourcer, float *sourcei, float *destr,float *desti,exint size);
    static void toTime2D(float *sourcer, float *sourcei, float *destr,float *desti,exint size);

private:

    int for_fft(float *yre, float *yim,float *zre,float *zim,exint npts);
    int inv_fft(float *yre, float *yim,float *zre,float *zim,exint npts);

    void dofft(int agamma);

    void fftcopy(float *pr_from, float *pi_from,float *p_to,
                 exint npts,int is_plus);
    
    void ifftcopy(float *p_from, float *pr_to,float *pi_to,
                  int agamma,fpreal factr,fpreal facti);
    void fftfill(float *p_co,exint part,exint all);
    
    void make_w(int neww);
    void butterfly(int fly,int agamma);

    static void init_bitsw();
    static int slow_bitsw(int k,int n);
    static exint bitswitch(exint k, exint n);
        
    float *x_co;
    float *w_co;
    
    static unsigned char tbl[256];
    int gamma;
    exint fft_size;

};
    
#endif

---