Tridiagonalization_Tridiagonalization_MatrixType.cpp revision c981c48f5bc9aefeffc0bcb0cc3934c2fae179dd
1544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack RusinMatrixXd X = MatrixXd::Random(5,5); 2544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack RusinMatrixXd A = X + X.transpose(); 3544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack Rusincout << "Here is a random symmetric 5x5 matrix:" << endl << A << endl << endl; 4544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack RusinTridiagonalization<MatrixXd> triOfA(A); 5544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack RusinMatrixXd Q = triOfA.matrixQ(); 6544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack Rusincout << "The orthogonal matrix Q is:" << endl << Q << endl; 7544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack RusinMatrixXd T = triOfA.matrixT(); 8544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack Rusincout << "The tridiagonal matrix T is:" << endl << T << endl << endl; 9544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack Rusincout << "Q * T * Q^T = " << endl << Q * T * Q.transpose() << endl; 10544dd4b11f7be76bb00fe29a60eaf2772dcc69caZack Rusin