test_matrix_operate.cpp

 
#include <cmath>
#include <Pds/Ra>
    
double func_vec(const Pds::Vector &V)
{
    return V.Sum();
}
 
double func_vec_par(const Pds::Vector &V,const Pds::Vector &C)
{
    return V.Dot(C);
}
 
Pds::Vector func_vec_trans(const Pds::Vector &V,const Pds::Vector &C)
{
    return Pds::MergeVer({V,C});
}
 
int main(void)
{
    Pds::Matrix A(3,2);
    Pds::Matrix B;
    Pds::Vector C;
 
    A.FillId();
 
    A.Print("\nA:\n");
 
    B=A.OperateCols(Pds::SumSquare);
    B.Print("\nA.OperateCols(Pds::SumSquare):\n");
 
    B=A.OperateCols(func_vec);
    B.Print("\nA.OperateCols(func_vec):\n");
 
    C="2 1 1";
    B=A.OperateCols(func_vec_par,C);
    B.Print("\nA.OperateCols(func_vec_par,C):\n");
 
    A.Print("\nA:\n");
    
    B=A.OperateRows(Pds::SumSquare);
    B.Print("\nA.OperateRows(Pds::SumSquare):\n");
 
    B=A.OperateRows(Pds::RMS);
    B.Print("\nA.OperateRows(Pds::RMS):\n");
 
    B=A.OperateRows(func_vec);
    B.Print("\nA.OperateRows(func_vec):\n");
 
    C="2 1";
    B=A.OperateRows(func_vec_par,C);
    B.Print("\nA.OperateRows(func_vec_par,C):\n");
 
    A.Print("\nA:\n");
 
    C="2 1";
    B=A.TransformRows(func_vec_trans,C);
    B.Print("\nA.TransformRows(func_vec_trans,C):\n");
 
    return 0;
}