PolynomialUtils1.cpp
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include <unsupported/Eigen/Polynomials>
	#include <iostream>

	using namespace Eigen;
	using namespace std;

	int main()
	{
	Vector4d roots = Vector4d::Random();
	cout << "Roots: " << roots.transpose() << endl;
	Eigen::Matrix<double,5,1> polynomial;
	roots_to_monicPolynomial( roots, polynomial );
	cout << "Polynomial: ";
	for( int i=0; i<4; ++i ){ cout << polynomial[i] << ".x^" << i << "+ "; }
	cout << polynomial[4] << ".x^4" << endl;
	Vector4d evaluation;
	for( int i=0; i<4; ++i ){
	evaluation[i] = poly_eval( polynomial, roots[i] ); }
	cout << "Evaluation of the polynomial at the roots: " << evaluation.transpose();
	}