test_matrix_operators_binary.cpp

 
#include <iostream>
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix A(2,2);
    Pds::Matrix B(2,2);
    Pds::Matrix B1(1,2,3);
    Pds::Matrix C;
    Pds::Vector V(2);
    Pds::Matrix D(3,3);
 
    
    A.Fill(3.0);    std::cout<<"A:\n"<<A;
    B.Fill(2.0);    std::cout<<"B:\n"<<B;
    V.Fill(1.0);    std::cout<<"V:\n"<<V;
    D.Fill(2.0);    std::cout<<"D:\n"<<D;    
    
    // + operator
    std::cout<<std::endl;
    
    C=A.Add(B);
    C.Print("A.Add(B)\n");
    
    C=A+B;
    C.Print("A+B\n");
    
    C=2-A+5;
    C.Print("2-A+5\n");
    
    // * operator
    std::cout<<std::endl;
    
    C=A.Mul(B);
    C.Print("A.Mul(B)\n");
    
    C=A.Mul(V);
    C.Print("A.Mul(V)\n");
    
    Pds::Matrix F(4,3); F.Fill(1.0); std::cout<<"F:\n"<<F;
    Pds::Matrix G(4,2); G.Fill(1.0); std::cout<<"G:\n"<<G;
    C=F.TMul(G);
    C.Print("F.TMul(G)\n");
    
    C=D.MulComp(-1,V);
    C.Print("D.MulComp(-1,V)\n");
        
    C=A*B;
    C.Print("A*B\n");
    
    // / operator
    std::cout<<std::endl;
    
    C=A.Div(2);
    C.Print("A.Div(2)\n");
    
    C=A/2;
    C.Print("A/2\n");
    
    // Geq y Leq
    C=Pds::MatrixId(4,4);
    C.Print("\nC\n");
 
    C.Geq(8).Print("\nC.Geq(8)\n");
 
    C.Leq(8).Print("\nC.Leq(8)\n");
 
    A.Print("\nA\n");
    B1.Print("B1\n");
 
    C=A.AddRowMatrix(B1);
    C.Print("A.AddRowMatrix(B1)\n");
 
    C=A.SubRowMatrix(B1);
    C.Print("A.SubRowMatrix(B1)\n");
    
    // (*) operator
    std::cout<<std::endl;
    
    A.Print("\nA\n");
    B.Print("B\n");
    
    C=B.Product(A.Product(B));
    C.Print("B.Product(A.Product(B))\n");
        
    C=B&A&B;
    C.Print("B&A&B\n");
    
    C=A.Div(B);
    C.Print("A.Division(B)\n");
        
    C=A/B;
    C.Print("A/B\n");
   
    C=A.Pow(B);
    C.Print("A.Pow(B)\n");
        
    return 0;
}