diff --git a/test/test_roots.cpp b/test/test_roots.cpp
index 854ff1eee..a55c12a1e 100644
--- a/test/test_roots.cpp
+++ b/test/test_roots.cpp
@@ -381,15 +381,14 @@ void test_complex_newton()
 // Polynomials which didn't factorize using Newton's method at first:
 void test_daubechies_fails()
 {
+    std::cout << "Testing failures from Daubechies filter computation.\n";
     using std::abs;
     using std::sqrt;
     using boost::math::tools::complex_newton;
     using boost::math::tools::polynomial;
     using boost::math::constants::half;
 
-    std::cout.precision(std::numeric_limits<double>::digits10+3);
-    //std::cout << std::hexfloat;
-    double tol = std::numeric_limits<double>::epsilon();
+    double tol = 500*std::numeric_limits<double>::epsilon();
     polynomial<std::complex<double>> p{{-185961388.136908293,141732493.98435241}, {601080390,0}};
     std::complex<double> guess{1,1};
     polynomial<std::complex<double>> p_prime = p.prime();