example_vector_dsp.cpp

 
#include <Pds/Ra>
#include <cmath>
    
int main(void)
{
    Pds::Vector R;
    Pds::Vector A("1 -1 5 -1 1 2 -1");
    Pds::Vector B("1 2 3");
        
    A.T().Print("\nA.T():\n");  // Transpose of A
    B.T().Print("\nB.T():\n");  // Transpose of B
 
    
    R=A.Conv(B);        // Conv(A,B)
    R.T().Print("\nA.Conv(B).T():\n");
    
    R=A.Conv1D(B,1,Pds::Ra::Full);   // Conv(A,B)
    R.T().Print("\nA.Conv1D(B,2,Pds::Ra::Full).T():\n");
    
    R=A.Conv(B,Pds::Ra::Same);   // Conv(A,B) with size of A
    R.T().Print("\nA.Conv(B,Pds::Ra::Same).T():\n");
    
    R=A.Conv1D(B,1,Pds::Ra::Same);   // Conv(A,B) with size of A
    R.T().Print("\nA.Conv1D(B,1,Pds::Ra::Same).T():\n");
    std::cout<<Pds::Vector::Conv1DOutputSize(A.Nlin(),B.Nlin(),1,Pds::Ra::Same)<<"\n";
    
    R=A.Conv1D(B,1,Pds::Ra::Full);   // Conv(A,B) with size of A
    R.T().Print("\nA.Conv1D(B,1,Pds::Ra::Full).T():\n");
    std::cout<<Pds::Vector::Conv1DOutputSize(A.Nlin(),B.Nlin(),1,Pds::Ra::Full)<<"\n";
 
    R=A.Conv1D(B,1,Pds::Ra::Valid);   // Conv(A,B) with size of A
    R.T().Print("\nA.Conv1D(B,1,Pds::Ra::Valid).T():\n");
    std::cout<<Pds::Vector::Conv1DOutputSize(A.Nlin(),B.Nlin(),1,Pds::Ra::Valid)<<"\n";
    
    R=A.XCorr(B);       // XConv(A,B)
    R.T().Print("\nA.XCorr(B).T():\n");
    
    R=A.XCorr(B,true);  // XConv(A,B) with size of A
    R.T().Print("\nA.XCorr(B,true).T():\n");
    
    std::vector<Pds::Vector> Block;
    Block=Pds::Vector::MaxPoolingVectorBlock({A,A},2);
    A.T().Print("\nA.T():\n");  // Transpose of A
    Block[0].T().Print("Block[0].T():\t");
    Block[1].T().Print("Block[1].T():\t");
    
    Pds::Vector::FlattenVectorBlock(Block).T().Print("Pds::Vector::FlattenVectorBlock(Block).T():\t");
    return 0;
}