00001 MatrixXcf A = MatrixXcf::Random(4,4); 00002 HessenbergDecomposition<MatrixXcf> hd(4); 00003 hd.compute(A); 00004 cout << "The matrix H in the decomposition of A is:" << endl << hd.matrixH() << endl; 00005 hd.compute(2*A); // re-use hd to compute and store decomposition of 2A 00006 cout << "The matrix H in the decomposition of 2A is:" << endl << hd.matrixH() << endl;