example_matrix_algebra.cpp

Programa para el testeo de las funciones.Para compilar o código example_matrix_algebra.cpp:

g++ -static -o example_matrix_algebra example_matrix_algebra.cpp -lpdsramm 

Para executar o programa:

./example_matrix_algebra

Retornando por consola:

A:
1       2       1
1       -1      1
1       1       2
    
ID.T():
1       3       2
    
V=A.Multiindex(ID) :
8
-1
4

Código 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.MultipleMse(Block).Print("A.MultipleMse(Block):\n");
    
    return 0;
}