example_goertzel_create.cpp

 
#include <Pds/Ca>
#include <Pds/Ra>
    
int main(void)
{
    unsigned int N=512;
    unsigned int K=8;
    unsigned int k;
    double fc=4;
    
    Pds::Vector T=Pds::LinSpace(0,1,N);
    Pds::Vector X=Pds::Sin(2*Pds::Ra::Pi*fc*T);
    
    Pds::Goertzel FT[K];
    Pds::Complex Out[K];
    
    for(k=0;k<K;k++)
    FT[k].Configure(k,N);
    
    for(k=0;k<K;k++)
    Out[k]=FT[k].Evaluate(X);
    
    for(k=0;k<K;k++)
    std::cout<<"FT{X}["<<k<<"]: "<<Out[k].Module()<<std::endl;
    
}