example_matrix_dsp_conv.cpp

 
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix Ans;
    
    Pds::Matrix K2( "-2  1\n"
                    "-1  1\n");
    K2.Print("K2:\n");
    
    Pds::Matrix K3( "-2 -1  2\n"
                    "-1  2 -1\n"
                    " 1  0  1");
    K3.Print("K3:\n");
    
    Pds::Matrix A(  "-1  1 -1 -1  1\n"
                    " 1  1  1 -1 -1\n"
                    " 1 -1  1 -1  1\n"
                    "-1  1  1 -1 -1\n"
                    " 1 -1 -1 -1  1\n");
    A.Print("A:\n");
    
    std::cout<<"\n";
    
    Ans=A.Conv(K2,Pds::Ra::Full);
    Ans.Print("A*^{F}K2:\n");
    
    Ans=A.Conv(K2,Pds::Ra::Same);
    Ans.Print("A*^{S}K2:\n");
    
    Ans=A.Conv(K2,Pds::Ra::Valid);
    Ans.Print("A*^{V}K2:\n");
    
    std::cout<<"\n";
    
    Ans=A.Conv(K3,Pds::Ra::Full);
    Ans.Print("A*^{F}K3:\n");
    
    Ans=A.Conv(K3,Pds::Ra::Same);
    Ans.Print("A*^{S}K3:\n");
    
    Ans=A.Conv(K3,Pds::Ra::Valid);
    Ans.Print("A*^{V}K3:\n");
    
    return 0;
}