| Brian Silverman | 72890c2 | 2015-09-19 14:37:37 -0400 | [diff] [blame^] | 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; |