test_matrix_find.cpp

 
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix A = Pds::Eye(5);
    
    A.Print("\nA:\n");
    std::cout<<"Pds::Find(A): "<<Pds::Ra::IndicesToString(Pds::Find(A),", ")<<std::endl;
    
    Pds::Matrix X;
    Pds::Matrix Y;
    Pds::Meshgrid(Pds::LinSpace(0,6,7),Pds::LinSpace(-2,2,5),X,Y);    
    X.Print("\nX:\n");
    
    Pds::Vector B("1 1 1 0 0 0 1");
    
    Pds::Matrix X0=X.FindRows(B);
    X0.Print("\nX0=X.FindRows(B):\n");
    
    Pds::Matrix X1=X.FindRows(B,0.5);
    X1.Print("\nX1=X.FindRows(B,0.5):\n");
    
    std::vector<Pds::Matrix> D;
    D=X.FindRowsBlock(B);
    D[0].Print("\nD[0]:\n");
    D[1].Print("\nD[1]:\n");
    
    X0.Geq(2.5).Print("\nX0.Geq(2.5)\n");
 
    std::cout<<"X0.BalancedBinaryWeight(2.5):\n"<<X0.BalancedBinaryWeight(2.5);    
    std::cout<<"X0.P1Probability(2.5): "<<X0.P1Probability(2.5)<<std::endl;    
    std::cout<<"X0.BinaryEntropy(2.5): "<<X0.BinaryEntropy(2.5)<<std::endl;
    
    return 0;
}