example_matrix_algebra.cpp

 
#include <cmath>
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix A("1 2 1\n1 -2 1\n1 1 2");
    std::vector<Pds::Matrix> Block(3);
    Block[0]=Pds::Ones(2,2);
    Block[1]=Pds::Ones(2,2);
    Block[2]=Pds::Ones(2,2);
    A.Print("\nA:\n");
    
    Pds::Vector ID({1,3,2});    
    ID.T().Print("\nID.T():\n");
 
    Pds::Vector V=A.MultiIndex(ID);
 
    V.Print("\nV=A.MultiIndex(ID) :\n");
    
    std::cout<<"A.Det():\t"<<A.Det()<<std::endl;
 
    A.Print("\nA:\n");
    A.IdInMultipleMse(Block).Print("A.IdInMultipleMse(Block):\n");
    
    A.Print("\nA:\n");
    A.NormalizeRows();
    A.Print("\nA:\n");
    
    Pds::Matrix B("0.00109512 -0.00027612\n-0.00027612 6.962e-05");
    B.Print("\nB:\n");
    std::cout<<"B.Det():\t"<<B.Det()<<std::endl;
    std::cout<<"B.Det():\t"<<B.Get(0)*B.Get(3)-B.Get(1)*B.Get(2)<<std::endl;
    
    return 0;
}