diff options
author | Gael Guennebaud <g.gael@free.fr> | 2018-12-09 23:05:46 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2018-12-09 23:05:46 +0100 |
commit | cf697272e1326e5ceecfd876ea67943b88425521 (patch) | |
tree | 6cef7b4af29ab363455361583622650d3a153c82 /unsupported | |
parent | 450dc97c6b14cd738def377d8b04c12427c6449a (diff) |
Remove debug code.
Diffstat (limited to 'unsupported')
-rw-r--r-- | unsupported/Eigen/src/Polynomials/PolynomialSolver.h | 4 |
1 files changed, 1 insertions, 3 deletions
diff --git a/unsupported/Eigen/src/Polynomials/PolynomialSolver.h b/unsupported/Eigen/src/Polynomials/PolynomialSolver.h index d4f737134..5e0ecbb43 100644 --- a/unsupported/Eigen/src/Polynomials/PolynomialSolver.h +++ b/unsupported/Eigen/src/Polynomials/PolynomialSolver.h @@ -355,13 +355,11 @@ class PolynomialSolver : public PolynomialSolverBase<_Scalar,_Deg> companion.balance(); m_eigenSolver.compute( companion.denseMatrix() ); m_roots = m_eigenSolver.eigenvalues(); - MatrixXcd A = companion.denseMatrix(); // cleanup noise in imaginary part of real roots: // if the imaginary part is rather small compared to the real part // and that cancelling the imaginary part yield a smaller evaluation, // then it's safe to keep the real part only. - RealScalar coarse_prec = std::pow(4,poly.size()+1)*NumTraits<RealScalar>::epsilon(); - std::cout << coarse_prec << "\n"; + RealScalar coarse_prec = RealScalar(std::pow(4,poly.size()+1))*NumTraits<RealScalar>::epsilon(); for(Index i = 0; i<m_roots.size(); ++i) { if( internal::isMuchSmallerThan(numext::abs(numext::imag(m_roots[i])), |