eigen-2.1.1: eigen3/doc/snippets/MatrixBase_array.cpp

Vector3d v(1,2,3);
v.array() += 3;
v.array() -= 2;
cout << v << endl;