HessenbergDecomposition_packedMatrix.cpp revision c981c48f5bc9aefeffc0bcb0cc3934c2fae179dd
124943d2ee8bfaa7cf5893e4709143924157a5c1eChris LattnerMatrix4d A = Matrix4d::Random(4,4); 224943d2ee8bfaa7cf5893e4709143924157a5c1eChris Lattnercout << "Here is a random 4x4 matrix:" << endl << A << endl; 324943d2ee8bfaa7cf5893e4709143924157a5c1eChris LattnerHessenbergDecomposition<Matrix4d> hessOfA(A); 424943d2ee8bfaa7cf5893e4709143924157a5c1eChris LattnerMatrix4d pm = hessOfA.packedMatrix(); 524943d2ee8bfaa7cf5893e4709143924157a5c1eChris Lattnercout << "The packed matrix M is:" << endl << pm << endl; 624943d2ee8bfaa7cf5893e4709143924157a5c1eChris Lattnercout << "The upper Hessenberg part corresponds to the matrix H, which is:" 724943d2ee8bfaa7cf5893e4709143924157a5c1eChris Lattner << endl << hessOfA.matrixH() << endl; 824943d2ee8bfaa7cf5893e4709143924157a5c1eChris LattnerVector3d hc = hessOfA.householderCoefficients(); 924943d2ee8bfaa7cf5893e4709143924157a5c1eChris Lattnercout << "The vector of Householder coefficients is:" << endl << hc << endl; 1024943d2ee8bfaa7cf5893e4709143924157a5c1eChris Lattner