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,7,8),Pds::LinSpace(-2,2,5),X,Y);    
    X.Print("\nX:\n");
    
    Pds::Vector B("1 1 1 0 0 0 1 0");
    
    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");
    
    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;
}