example_fft_evaluate2.cpp

 
#include <Pds/Ca>
#include <Pds/Ra>
    
int main(void)
{
    unsigned int N=16;
    unsigned int Nlin=4;
    
    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::Matrix  Xmat(Nlin,N);
    Pds::CMatrix Ymat(Nlin,N);
    
    for(unsigned int lin=0;lin<Xmat.Nlin();lin++)
        Xmat.SetRowVector(lin,X);
    
    
    FT.Evaluate(Xmat,Ymat);
    
    Ymat.Module().T().Print("Y.Module().T():\n");
    
}