test_matrix_statistics.cpp

 
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix A(200,200);
    A.FillRandN();
    
    Pds::Matrix B=A+2.0*Pds::RandU(200,200);
    
    
    std::cout<<"A.Max()  : "<<A.Max() <<std::endl;
    std::cout<<"A.Min()  : "<<A.Min() <<std::endl;
    std::cout<<"A.Sum()  : "<<A.Sum() <<std::endl;
    std::cout<<"A.Mean() : "<<A.Mean()<<std::endl;
    std::cout<<"A.Std()  : "<<A.Std() <<std::endl;
    std::cout<<"A.Corr(B): "<<A.Corr(B)<<std::endl;
    
    Pds::Matrix C("-2 2 2 \n-1 2 -1\n 1 0 3");
 
    Pds::Matrix D("1 1\n1 -1");
    
    C.Print("C:\n");
    D.Print("D:\n");
    
    C.XCorr(D).Print("C.XCorr(D):\n");
    C.XCorr(D,true).Print("C.XCorr(D,true):\n");
    
    C.Print("\nC:\n");
    Pds::Array<unsigned int> Lin;
    C.MinInCols(Lin).Print("C.MinInCols(Lin):\n");
    Lin.Print("Lin:\n");
    
    return 0;
}