15 #ifndef INCLUDED_numeric_linear_algebra_rgg_HH
16 #define INCLUDED_numeric_linear_algebra_rgg_HH
19 #include <fem/fem.hpp>
22 namespace linear_algebra {
24 using namespace fem::major_types;
63 arr_cref<double> alfr,
64 arr_cref<double> alfi,
65 arr_cref<double> beta,
void qzhes(int const &nm, int const &n, arr_ref< double, 2 > a, arr_ref< double, 2 > b, bool const &matz, arr_ref< double, 2 > z)
This subroutine is the first step of the QZ algorithm for solving generalized matrix eigenvalue probl...
void qzval(int const &nm, int const &n, arr_ref< double, 2 > a, arr_ref< double, 2 > b, arr_ref< double > alfr, arr_ref< double > alfi, arr_ref< double > beta, bool const &matz, arr_ref< double, 2 > z)
This subroutine is the third step of the QZ algorithm for solving generalized matrix eigenvalue probl...
void qzit(int const &nm, int const &n, arr_ref< double, 2 > a, arr_ref< double, 2 > b, double const &eps1, bool const &matz, arr_ref< double, 2 > z, int &ierr)
This subroutine is the second step of the QZ algorithm for solving generalized matrix eigenvalue prob...
void qzvec(int const &nm, int const &n, arr_cref< double, 2 > a, arr_ref< double, 2 > b, arr_cref< double > alfr, arr_cref< double > alfi, arr_cref< double > beta, arr_ref< double, 2 > z)
This subroutine is the optional fourth step of the QZ algorithm for solving generalized matrix eigenv...
double epslon(double const &x)
Estimate unit roundoff in quantities of size x.
void rgg(int const &nm, int const &n, arr_ref< double, 2 > a, arr_ref< double, 2 > b, arr_ref< double > alfr, arr_ref< double > alfi, arr_ref< double > beta, int const &matz, arr_ref< double, 2 > z, int &ierr)
This subroutine calls the recommended sequence of subroutines from the eigensystem subroutine package...