| 65 | |
| 66 | template<typename _Scalar, int _Deg> |
| 67 | void CauchyBounds(int deg) |
| 68 | { |
| 69 | typedef internal::increment_if_fixed_size<_Deg> Dim; |
| 70 | typedef Matrix<_Scalar,Dim::ret,1> PolynomialType; |
| 71 | typedef Matrix<_Scalar,_Deg,1> EvalRootsType; |
| 72 | |
| 73 | PolynomialType pols(deg+1); |
| 74 | EvalRootsType roots = EvalRootsType::Random(deg); |
| 75 | roots_to_monicPolynomial( roots, pols ); |
| 76 | _Scalar M = cauchy_max_bound( pols ); |
| 77 | _Scalar m = cauchy_min_bound( pols ); |
| 78 | _Scalar Max = roots.array().abs().maxCoeff(); |
| 79 | _Scalar min = roots.array().abs().minCoeff(); |
| 80 | bool eval = (M >= Max) && (m <= min); |
| 81 | if( !eval ) |
| 82 | { |
| 83 | cerr << "Roots: " << roots << endl; |
| 84 | cerr << "Bounds: (" << m << ", " << M << ")" << endl; |
| 85 | cerr << "Min,Max: (" << min << ", " << Max << ")" << endl; |
| 86 | } |
| 87 | VERIFY( eval ); |
| 88 | } |
| 89 | |
| 90 | template<typename _Scalar> void CauchyBounds_scalar() |
| 91 | { |
nothing calls this directly
no test coverage detected