example_ra_indices.cpp

 
#include <Pds/Ra>
#include <ctime>
 
int main(void)
{
    Pds::Matrix A(5,5);
    A.FillId();
    A.Print("\nA:\n");
    //(A/Pds::Ra::Pi-4).Print("\nA:\n");
 
    std::srand(std::time(0));
    
    // Getting index vector
    std::vector<unsigned int> Vec=Pds::Ra::IndicesRandVector(0,4,3);
    std::cout<<"Vec:\t"<<Pds::Ra::IndicesToString(Vec,",")<<"\n";
    
    // Getting index list
    std::list<unsigned int> List=Pds::Ra::IndicesRandList(0,4,3);
    std::cout<<"List:\t"<<Pds::Ra::IndicesToString(Vec,",")<<"\n";
 
    // Getting rows from matrix and index list
    Pds::Matrix B=A.GetRows(List);
    B.Print("\nA.GetRows(List):\n");
 
    // Getting rows from matrix and index list
    B.Fill(-1.0);
    A.GetRows(3,5,B);
    B.Print("\nB:\n");
 
    // Getting columns from matrix and index list
    B=A.GetCols(List);
    B.Print("\nA.GetCols(List):\n");
    
    //               0 1 2 3 4 5 6 7 8 9 10 11 12 13
    //Pds::Vector V="1 0 1 1 1 0 0 0 1 0  1  1  0  0";
    //Pds::Vector V="1 1 1 1 1 1 1 1 1 1  1  1  1  1";
    Pds::Vector   V="0 0 0 0 0 0 0 0 0 0  0  0  0  0";
    V.T().Print("\nV.T():\n");
    std::vector<unsigned int> ID;
    ID=Pds::Ra::IndicesRandVector(V,0.5,8);
    std::cout<<Pds::Ra::IndicesToString(ID,",")<<std::endl;
    
    std::cout<<"ID:\n";
    std::string str=Pds::Ra::IndicesXmlToString({1, 2 ,5,0},"Indices");
    std::cout<<str;
    ID=Pds::Ra::IndicesXmlFromString(str,"Indices");
    std::cout<<Pds::Ra::IndicesToString(ID,",")<<"\n";
    
    unsigned int lin=0; 
    unsigned int col=0;
    std::cout<<Pds::Ra::Ind2Sub(Pds::Size(3,2),5,lin,col)<<"\n";
    std::cout<<lin<<"\n";
    std::cout<<col<<"\n";
    
    return 0;
}