example_matrix_algebra4.cpp

 
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix A(3,3);
    A.FillId();    
    A.Print("A:\n");
    
    Pds::Matrix X=  " 1  1 -1\n"
                    "-1  2 -1\n"
                    " 2  1 -1\n"
                    " 2 -1 -1\n";
    X.Print("X:\n");
    
    Pds::Matrix X0=" 1 -1\n 2  2";
    Pds::Matrix X1=" 1  2\n 1 -1";
    Pds::Matrix X2="-1 -1\n-1 -1";
    
    Pds::Vector Mu=Pds::Ones(3,1);
    
    A.QuadraticFormWithRows(X,Mu).Print("(x_l-1)^t A (x_l-1):\n");
    
    A.QuadraticFormWithSamples({X0,X1,X2},Mu).Print("(x_l-1)^t A (x_l-1):\n");
    
    return 0;
}