Convolver.h

Go to the documentation of this file.

#ifndef INCLUDE_CONVOLUTION_H
#define INCLUDE_CONVOLUTION_H

#include "CSL_Core.h"
#include "Spectral.h"           // fftw includes and types


#define CHECK_PTR(ptr)                              \
    if (!ptr)    throw MemoryError("Can't allocate buffer")

#define initvec( name, size, type )                 \
    if ((name = (type *) calloc(size, sizeof(type))) == NULL)  \
        throw MemoryError("Can't allocate buffer");

#define cmac(in1, in2, out)                         \
    out[0] += in1[0] * in2[0] - in1[1] * in2[1];    \
    out[1] += in1[0] * in2[1] + in1[1] * in2[0];

#define cbinc(buf, size)    if (++buf > size-1)     buf = 0;
#define cbdec(buf, size)    if (--buf < 0)          buf = size - 1;

#define cbarb(buf, size, amt)                       \
    buf += amt;                                     \
    if (buf >= size) buf -= size;                   \
    else if (buf < 0) buf += size;

namespace csl {


class Convolver : public Effect {

public:                                     
    Convolver() : Effect() { };
    Convolver(char *IRfilename);            
    Convolver(unsigned len, Buffer & impulseResp);
    Convolver(unsigned len, char * IRfilename);
    Convolver(Buffer & inbuf);
    ~Convolver();
    void initialize(Buffer & buf) throw (CException);       
    void setFilters(fftwf_complex ** filterFFTs);           
    void setInputf(fftwf_complex * inFFT);                  
    CSL_FFTW_sample * mSampleBuffer;                        

    void nextBuffer(Buffer &outputBuffer, unsigned outBufNum) throw (CException);

protected:  
    CSL_FFTW_cmplx ** mFilterFFTs;                          
    CSL_FFTW_cmplx ** mInputFFTs;                           
    CSL_FFTW_cmplx * mSpectrumBuffer;                       
    CSL_FFTW_plan mForwardPlan, mInversePlan;               
    unsigned mFFTSize, mWindowCount, mNumBufs;              
    bool mMyBuffers;                                        
    bool mMyInput;                                          
    
    void initialize (const char *filename) throw (CException);  
    
    void complex_multiply_accumulate(fftwf_complex* left, fftwf_complex* right, fftwf_complex * output);

};

}

#endif

---