2 cout <<
"Here is a random 4x4 matrix:" << endl <<
A << endl;
3 HessenbergDecomposition<Matrix4d>
hessOfA(
A);
5 cout <<
"The packed matrix M is:" << endl <<
pm << endl;
6 cout <<
"The upper Hessenberg part corresponds to the matrix H, which is:"
7 << endl <<
hessOfA.matrixH() << endl;
9 cout <<
"The vector of Householder coefficients is:" << endl <<
hc << endl;
cout<< "Here is a random 4x4 matrix:"<< endl<< A<< endl;HessenbergDecomposition< Matrix4d > hessOfA(A)
Matrix< double, 4, 4 > Matrix4d
4×4 matrix of type double.
Matrix< double, 3, 1 > Vector3d
3×1 vector of type double.