test_matrix_static_array.cpp

 
#include <iostream>
#include <fstream>
#include <Pds/Ra>
#include <cmath>
    
int main(void)
{
    std::string str;
    Pds::Matrix A=Pds::Eye(3);
    
    unsigned int Nlin=2;
    unsigned int Ncol=2;
    double **array=NULL;
    
    char filepath[]="filedat.txt";
    
    
    
    array=Pds::Array<double>::ArrayAllocate(Nlin,Ncol);
    if(array==NULL)   return 0; 
    
    array[0][0]=0.0/0.0;    array[0][1]=0.0/1.0;
    array[1][0]=1.0/0.0;    array[1][1]=1.0/1.0;
    std::cout<<"\n";
    
    
    Pds::Array<double>::ArraySave(filepath,array,Nlin,Ncol);
    
    std::cout   <<"ArrayAllocate(Nlin,Ncol):\n"
                <<Pds::Array<double>::ArrayToString(array,Nlin,Ncol);
    
    Pds::Array<double>::ArrayRelease(array,Nlin);
    if(array==NULL) std::cout<<"Array released.\n\n\n";
    
    A.Print("A:\n");
    
    array=Pds::Array<double>::ArrayAllocate(A);
    if(array==NULL)   {std::cout<<"Error Alocate.\n";return 0;}
    else              std::cout<<"New array allocated.\n";
    
    std::cout   <<"ArrayAllocate(A):\n"
                <<Pds::Array<double>::ArrayToString(array,A.Nlin(),A.Ncol());
    
    Pds::Array<double>::ArrayRelease(array,A.Nlin());
    if(array==NULL) std::cout<<"Array released.\n\n\n";
    
    
    
    array=Pds::Matrix::ArrayAllocate(sin,A);
    if(array==NULL)   {std::cout<<"Error Alocate.\n";return 0;}
    else              std::cout<<"New array allocated.\n";
    
    std::cout   <<"ArrayAllocate(sin,A):\n"
                <<Pds::Array<double>::ArrayToString(array,A.Nlin(),A.Ncol());
    
    Pds::Array<double>::ArrayRelease(array,A.Nlin());
    if(array==NULL) std::cout<<"Array released.\n\n\n";
    
    
    
    array=Pds::Array<double>::ArrayLoad(filepath,Nlin,Ncol);
    if(array==NULL)   {std::cout<<"Error Alocate.\n";return 0;}
    else              std::cout<<"New array allocated.\n";
    
    std::cout   <<"ArrayLoad(filepath,Nlin,Ncol):\n"
                <<Pds::Array<double>::ArrayToString(array,Nlin,Ncol);
    
    Pds::Array<double>::ArrayRelease(array,Nlin);
    if(array==NULL) std::cout<<"Array released.\n\n\n";
    
    
    
    array=Pds::Array<double>::ArrayFromString("1 2 3\n4 5 6\n",Nlin,Ncol);
    if(array==NULL)   {std::cout<<"Error Alocate.\n";return 0;}
    else              std::cout<<"New array allocated.\n";
    
    std::cout   <<"ArrayFromString(\"1 2 3\\n4 5 6\\n\",Nlin,Ncol):\n"
                <<Pds::Array<double>::ArrayToString(array,Nlin,Ncol);
    
    Pds::Array<double>::ArrayRelease(array,Nlin);
    if(array==NULL) std::cout<<"Array released.\n\n\n";
    
    array=Pds::Array<double>::ArrayColFromString("1 2 3\n4 5 6\n",Nlin,Ncol);
    if(array==NULL)   {std::cout<<"Error Alocate.\n";return 0;}
    else              std::cout<<"New array allocated.\n";
    
    std::cout   <<"ArrayColFromString(\"1 2 3\\n4 5 6\\n\",Nlin,Ncol):\n"
                <<Pds::Array<double>::ArrayToString(array,Nlin,Ncol);
                
    std::cout<<"\n";
                
    unsigned int Nel;
    Pds::Ra::ElementsInFile(filepath,Nel);
    std::cout<<"Nel: "<<Nel<<"\n\n";
                        
    std::ifstream ifs;
    ifs.open (filepath);
    
    if(ifs.is_open())
    while(!ifs.eof())
    { 
        Nlin=Ncol=0;
        
        array=Pds::Array<double>::ArrayLoadLineFromStream(ifs,Nlin,Ncol);
        if(array!=NULL)
        std::cout<<"["<<Ncol<<"]:"<<Pds::Array<double>::ArrayToString(array,Nlin,Ncol);
        
        Pds::Array<double>::ArrayRelease(array,Nlin);
    }
    ifs.close();
    
    return 0;
}