example_matrix_operators_binary2.cpp

 
#include <iostream>
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix A(4,2);
    Pds::Vector V(2);
    Pds::Matrix D(2,3);
    
    Pds::Matrix C;
    Pds::Vector R;
 
    
    A.FillId();         std::cout<<"A:\n"<<A;
    V.Fill(1.0);        std::cout<<"V:\n"<<V;
    D.FillRandC(0.8);   std::cout<<"D:\n"<<D;    
 
    C=D.CompMul(-1,A);
    C.Print("D.CompMul(-1,A)\n");
    
    R=D.MulComp(-1,V);
    R.Print("D.MulComp(-1,V)\n");
    
    R.Fill(0);
    D.MulComp(-1,V,R);
    R.Print("D.MulComp(-1,V,R)\n");
    
    A="1 -1\n 1 2";
    std::cout<<"A:\n"<<A;
    
    C=A;C.Fill(0);
    
    D.MulComp(-1,A,C);
    C.Print("D.MulComp(-1,A,C)\n");
    
    return 0;
}