example_fft_create.cpp

 
#include <Pds/Ca>
#include <Pds/Ra>
    
int main(void)
{
    unsigned int N=512;
    double fc=4;
    
    Pds::Vector T=Pds::LinSpace(0,1,N);
    Pds::Vector X=Pds::Sin(2*Pds::Ra::Pi*fc*T);
    
    Pds::Fft FT(N);
    if(FT.IsEmpty())    
    {
        pds_print_error_message("Error creando Pds::Fft.");
        return 0;
    }
    
    Pds::CVector Y(N);
    
    FT.Evaluate(X,Y);
    
    for(unsigned int k=0;k<8;k++)
        std::cout<<"|Y["<<k<<"]|: "<<Y[k].Module()<<std::endl;
    
}