12   Vector5d roots = Vector5d::Random();  in main()  local
13   cout << "Roots: " << roots.transpose() << endl;  in main()
15   roots_to_monicPolynomial( roots, polynomial );  in main()
18   cout << "Complex roots: " << psolve.roots().transpose() << endl;  in main()
23   cout << "Real roots: " << mapRR.transpose() << endl;  in main()
33   cout << "Complex roots: " << psolvef.roots().transpose() << endl;  in main()
35 …for( int i=0; i<6; ++i ){ evals[i] = std::abs( poly_eval( hardCase_polynomial, psolvef.roots()[i] …  in main()
36 …cout << "Norms of the evaluations of the polynomial at the roots: " << evals.transpose() << endl <…  in main()
41   cout << "Complex roots: " << psolve6d.roots().transpose() << endl;  in main()
44     std::complex<float> castedRoot( psolve6d.roots()[i].real(), psolve6d.roots()[i].imag() );  in main()
47 …cout << "Norms of the evaluations of the polynomial at the roots: " << evals.transpose() << endl <…  in main()
50 …cout << "The last root in float then in double: " << psolvef.roots()[5] << "\t" << psolve6d.roots(…  in main()
51   std::complex<float> castedRoot( psolve6d.roots()[5].real(), psolve6d.roots()[5].imag() );  in main()
52   cout << "Norm of the difference: " << std::abs( psolvef.roots()[5] - castedRoot ) << endl;  in main()