test_matrix_apply2.cpp

 
#include <cmath>
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix A(3,2);
 
    A.FillId();
    A.Print("\nA:\n");
    
    auto func=[](double x, const std::vector<double> &a)->double
    {
        if(a.size()>=2) return a[0]*x+a[1];
        return 0;
    };
    
    A.Apply(func,{2,-1});
    A.Print("\nA:\n");
    
    A.Apply(2,A);
    A.Print("\nA:\n");
    
    
    return 0;
}