1 Matrix4d A
= Matrix4d::Random(4,4);
2 cout
<< "Here is a random 4x4 matrix:" << endl
<< A
<< endl
;
3 HessenbergDecomposition
<Matrix4d
> hessOfA(A
);
4 Matrix4d pm
= hessOfA
.packedMatrix();
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
;
8 Vector3d hc
= hessOfA
.householderCoefficients();
9 cout
<< "The vector of Householder coefficients is:" << endl
<< hc
<< endl
;