test_cmatrix_binary_operator.cpp

 
#include <Pds/Ca>
    
int main(void)
{
    Pds::CMatrix A(2,2,3,4);
    Pds::CMatrix B(2,2,1,1);
    Pds::Complex c(1,1);
    Pds::CMatrix R;
    
    A.Print("A:\n");
    B.Print("B:\n");
 
    A.Add(B).Print("A.Add(B):\n");
    (A+B).Print("A+B:\n");
    A.Add(2).Print("A.Add(2):\n");
    (A+2).Print("A+2:\n");
    (2+A).Print("2+A:\n");
    A.Add(c).Print("A.Add(c):\n");
 
    std::cout<<"-----------------"<<std::endl;
    
    A.Print("A:\n");
    B.Print("B:\n");
 
    A.Sub(B).Print("A.Sub(B):\n");
    (A-B).Print("A-B:\n");
    A.Sub(2).Print("A.Sub(2):\n");
    (A-2).Print("A-2:\n");
    (2-A).Print("2-A:\n");
    A.Sub(c).Print("A.Sub(c):\n");
 
    std::cout<<"-----------------"<<std::endl;
    
    A.Print("A:\n");
    B.Print("B:\n");
 
    A.Mul(B).Print("A.Mul(B):\n");
    (A*B).Print("A*B:\n");
    A.Mul(2).Print("A.Mul(2):\n");
    (A*2).Print("A*2:\n");
    (2*A).Print("2*A:\n");
    A.Mul(c).Print("A.Mul(c):\n");
 
    std::cout<<"-----------------"<<std::endl;
    
    A.Print("A:\n");
    B.Print("B:\n");
 
    A.Div(B).Print("A.Div(B):\n");
    (A/B).Print("A/B:\n");
    A.Div(2).Print("A.Div(2):\n");
    (A/2).Print("A/2:\n");
    (2/A).Print("2/A:\n");
    (c/A).Print("c/A:\n");
 
    std::cout<<"-----------------"<<std::endl;
    
    A.Print("A:\n");
    B.Print("B:\n");
 
    A.Pow(B).Print("A.Pow(B):\n");
    A.Pow(2).Print("A.Pow(2):\n");
 
    return 0;
}