example_matrix_algebra3.cpp

 
#include <Pds/Ra>
    
int main(void)
{
    Pds::Matrix A(3,3);
    A.FillId();    
    A.Print("A:\n");
    
    Pds::Vector x="1 1 -1";
    x.Print("x:\n");
    
    std::cout<<"x^t A x: "<<A.QuadraticForm(x)<<"\n";
    
    Pds::Matrix X=  " 1  1 -1\n"
                    "-1  2 -1\n"
                    " 2  1 -1\n"
                    " 2 -1 -1\n";
    X.Print("X:\n");
    A.QuadraticFormWithRows(X).Print("x_l^t A x_l:\n");
    
    Pds::Vector Mu=Pds::Ones(3,1);
    
    A.QuadraticFormWithRows(X,Mu).Print("(x_l-1)^t A (x_l-1):\n");
    
    X.QuadraticDotWithRows(Mu).Print("(x_l-1)^t (x_l-1):\n");
    
    A.FillRandU(-4,4);
    A.Print("A:\n");
    
    std::cout<<"              A.MinOfDiag(): "<<A.MinOfDiag()<<"\n";
    std::cout<<"           A.MinAbsOfDiag(): "<<A.MinAbsOfDiag()<<"\n";
    std::cout<<"           A.MaxAbsOfDiag(): "<<A.MaxAbsOfDiag()<<"\n";
    std::cout<<"A.MaxAbsOfLowerTriangular(): "<<A.MaxAbsOfLowerTriangular()<<"\n";
    
    return 0;
}