HessenbergDecomposition_compute.cpp revision c981c48f5bc9aefeffc0bcb0cc3934c2fae179dd
1MatrixXcf A = MatrixXcf::Random(4,4); 2HessenbergDecomposition<MatrixXcf> hd(4); 3hd.compute(A); 4cout << "The matrix H in the decomposition of A is:" << endl << hd.matrixH() << endl; 5hd.compute(2*A); // re-use hd to compute and store decomposition of 2A 6cout << "The matrix H in the decomposition of 2A is:" << endl << hd.matrixH() << endl; 7