12#ifndef TLAPACK_LAHQR_SCHUR22_HH
13#define TLAPACK_LAHQR_SCHUR22_HH
43template <TLAPACK_REAL T>
85 const T
bcmax = max(abs(
b), abs(c));
153 const T
sab = sqrt(abs(
b));
154 const T
sac = sqrt(abs(c));
173 const T
temp = sqrt(abs(
b)) * sqrt(abs(c));
real_type< TX, TY > lapy2(const TX &x, const TY &y)
Finds , taking care not to cause unnecessary overflow.
Definition lapy2.hpp:34
void lahqr_schur22(T &a, T &b, T &c, T &d, complex_type< T > &s1, complex_type< T > &s2, T &cs, T &sn)
Computes the Schur factorization of a 2x2 matrix A.
Definition lahqr_schur22.hpp:44
Sort the numbers in D in increasing order (if ID = 'I') or in decreasing order (if ID = 'D' ).
Definition arrayTraits.hpp:15
typename traits::real_type_traits< Types..., int >::type real_type
The common real type of the list of types.
Definition scalar_type_traits.hpp:113
constexpr int sgn(const T &val)
Type-safe sgn function.
Definition utils.hpp:109