1 #include <iostream> 2 #include <Eigen/Dense> 3 4 using namespace Eigen; 5 using namespace std; main()6int main() 7 { 8 Vector3d v(1,2,3); 9 Vector3d w(0,1,2); 10 11 cout << "Dot product: " << v.dot(w) << endl; 12 double dp = v.adjoint()*w; // automatic conversion of the inner product to a scalar 13 cout << "Dot product via a matrix product: " << dp << endl; 14 cout << "Cross product:\n" << v.cross(w) << endl; 15 } 16