example_fir_create.cpp

 
#include <iomanip>
#include <Pds/Ra>
#include <Pds/Sp>
    
int main(void)
{
    double result;
    Pds::Fir FIR(Pds::Vector(100,0.01));
 
    double Tmax=10;
    double fs=1.0;
    unsigned int N=10000000;
    
    Pds::Vector T=Pds::LinSpace(0,Tmax,N);
    Pds::Vector X=Pds::Sin((2*Pds::Ra::Pi*fs)*T);
 
    Pds::Ra::Tic();
    for(unsigned int i=0;i<N;i++)
        result=FIR.Evaluate(X.Get(i));
    Pds::Ra::Toc();
    
    std::cout<<"Final:"<<result<<"\n";
}