FFT.cxx

Go to the documentation of this file.

/*
 * Copyright (c) 2001-2004 MUSIC TECHNOLOGY GROUP (MTG)
 *                         UNIVERSITAT POMPEU FABRA
 *
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program 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 General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */


#include "FFT.hxx"
#include "Audio.hxx"
#include "Spectrum.hxx"
#include "SpectrumConfig.hxx"
#include "ProcessingFactory.hxx"

namespace CLAM 
{

namespace Hidden
{
        static const char* metadata[] = {
                "key", "FFT",
                "category", "Analysis",
                "description", "FFT",
                0
        };
        static FactoryRegistrator<ProcessingFactory, FFT> reg = metadata;
}

        SpecTypeFlags FFT_base::mComplexflags;

        FFT_base::FFT_base() :
                mSize(0),
                mInput("Audio Input",this),
                mOutput("Spectrum Output",this),
                fftbuffer( NULL )
        {
        };

        FFT_base::~FFT_base()
        {
                if (fftbuffer) delete [] fftbuffer;
        }

        void FFT_base::ChangeSize(int n)
        {
                CLAM_ASSERT(n>=0,"Wrong (negative) Size in control input.");

                CLAM_ASSERT(false, "Controls not yet implemented.");

        }

        void FFT_base::CheckTypes(const Audio& in, const Spectrum &out) const
        {

                CLAM_BEGIN_CHECK
                // Input object checking
                if (in.GetSize()!=mSize) { 
                        std::stringstream ss;
                        ss << "FFT::Do: Wrong size in FFT Audio input\n"
                           << "  Expected: " << mSize << ", used " << in.GetSize();
                        CLAM_ASSERT(0,ss.str().c_str());
                }
                if (!in.HasBuffer())
                        CLAM_ASSERT(0,"FFT Do: Float attribute required for Audio object.");
                if (out.GetSize() != mSize/2+1 ) { // ALGORITHM DEPENDENT CHECKING
                        std::stringstream ss;
                        ss << "FFT::Do: wrong size  Spectrum.\n"
                           << "  Expected: " << mSize/2+1 << ", used " << out.GetSize();
                        CLAM_ASSERT(0,ss.str().c_str());
                }

                CLAM_END_CHECK
        }

        bool FFT_base::SetPrototypes(const Audio& in,const Spectrum &out)
        {
                CheckTypes(in,out);

                SpecTypeFlags flags;
                out.GetType(flags);

                if (flags.bComplex)
                        if (flags.bPolar || flags.bMagPhase || flags.bMagPhaseBPF)
                                mState=sComplexSync;
                        else
                                mState=sComplex;
                else
                        if (flags.bPolar || flags.bMagPhase || flags.bMagPhaseBPF)
                                mState=sOther;
                        else
                                CLAM_ASSERT(false,"FFT_numrec: SetPrototypes(): Spectrum with no attributes!");

                return true;
        }

        bool FFT_base::ConcreteConfigure(const ProcessingConfig& c)
        {
                int oldSize = mSize;
                
                CopyAsConcreteConfig(mConfig, c);
                if (mConfig.HasAudioSize()) {
                        CLAM_ASSERT(mSize>=0, "Negative Size in FFT configuration");
                        mSize = mConfig.GetAudioSize();
                        if(mSize>0)
                        {
                                mInput.SetSize( mSize );
                                mInput.SetHop( mSize );
                        }
                }

                CLAM_ASSERT(mSize>=0, "Negative Size in FFT configuration");

                mState=sOther;
                mComplexflags.bComplex=1;
                mComplexflags.bMagPhase=0;
                if (mSize == 0) 
                {
                        fftbuffer = 0;
                        return false;
                }
                if (mSize == oldSize)
                        return true;

                delete [] fftbuffer;
                fftbuffer = new TData[mSize];   

                SpectrumConfig cfg;
                SpecTypeFlags fl;
                fl.bMagPhase=0;
                fl.bComplex=1;
                cfg.SetType(fl);
                cfg.SetSize(mSize);
                mComplexSpectrum.Configure(cfg);                

                return true;
        }

        bool FFT_base::UnsetPrototypes()
        {
                mState=sOther;
                return true;
        }

        void FFT_base::ToOther(Spectrum &out)
        {
                if(out.HasComplexArray()) {
                        ToComplex(out);
                        out.SynchronizeTo(mComplexflags);
                }
                else {
                        ToComplex(mComplexSpectrum);
                        out.SynchronizeTo(mComplexSpectrum);
                }
        }       
};//namespace CLAM

---