HessenbergDecomposition_packedMatrix.cpp 482 B

123456789
  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;