13 #ifndef eli_mutil_poly_root_radius_hpp
14 #define eli_mutil_poly_root_radius_hpp
18 #include "eli/code_eli.hpp"
30 template<
typename data__>
33 data__ term, max_term(0), radius(1);
43 term=std::abs(coef(i)/coef(deg));
52 template<
typename data__>
55 data__ term1, term2, radius(1);
63 term2=std::pow(std::abs(coef[0]/coef[deg]), 1/static_cast<data__>(deg));
66 term1=deg*std::abs(coef[0]/coef[1]);
67 return std::min(term1, term2);
data__ max_radius(const polynomial< data__ > &f)
Definition: radius.hpp:31
Definition: polynomial.hpp:31
Eigen::Matrix< data_type, Eigen::Dynamic, 1 > coefficient_type
Definition: polynomial.hpp:35
index_type degree() const
Definition: polynomial.hpp:72
void get_coefficients(coefficient_type &aout) const
Definition: polynomial.hpp:86
coefficient_type::Index index_type
Definition: polynomial.hpp:36
data__ at_least_radius(const polynomial< data__ > &f)
Definition: radius.hpp:53