example_matrix_dsp.cpp

 
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix Ans;
    
    Pds::Matrix Kernel( "-2 -1  2\n"
                        "-1  2 -1\n"
                        " 1  0  1");
    Kernel.Print("Kernel:\n");
    
    Pds::Matrix A(  "-1  1 -1 -1\n"
                    " 1  1  1 -1\n"
                    " 1 -1  1 -1\n"
                    " 1 -1 -1 -1");
    A.Print("A:\n");
    
    Ans=A.Conv2D(Kernel,{1,2},Pds::Ra::Full);
    Ans.Print("\nA.Conv2D(Kernel,{1,2},Pds::Ra::Full):\n");
    std::cout<<"Size: "<<Pds::Matrix::Conv2DOutputSize(A.Size(),Kernel.Size(),{1,2},Pds::Ra::Full)<<"\n";
    
    Ans=A.Conv2D(Kernel,{1,2},Pds::Ra::Same);
    Ans.Print("\nA.Conv2D(Kernel,{1,2},Pds::Ra::Same):\n");
    std::cout<<"Size: "<<Pds::Matrix::Conv2DOutputSize(A.Size(),Kernel.Size(),{1,2},Pds::Ra::Same)<<"\n";
    
    Ans=A.Conv2D(Kernel,{1,2},Pds::Ra::Valid);
    Ans.Print("\nA.Conv2D(Kernel,{1,2},Pds::Ra::Valid):\n");
    std::cout<<"Size: "<<Pds::Matrix::Conv2DOutputSize(A.Size(),Kernel.Size(),{1,2},Pds::Ra::Valid)<<"\n";
    
    Ans=Pds::Conv2DCh({A,A},{Kernel,Kernel},{1,1},Pds::Ra::Same);
    Ans.Print("\nPds::Conv2DCh({A},{Kernel},{1,1},Pds::Ra::Valid):\n");
    
    Pds::Matrix B=Pds::RandN(7,7);
    B.Print("B:\n");
    B.MaxPooling(3,3).Print("B.MaxPooling(3,3):\n");
    
    Pds::SampleBlock Block;
    Block=Pds::Matrix::MaxPoolingSampleBlock({B,B},3,3);
    Block[0].Print("Block[0]:\n");
    Block[1].Print("Block[1]:\n");
    
    return 0;
}