1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <iostream>
2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <Eigen/Dense>
3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathusing namespace Eigen;
5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath
6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathint main()
7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{
8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  Matrix2d a;
9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  a << 1, 2,
10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath       3, 4;
11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  MatrixXd b(2,2);
12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  b << 2, 3,
13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath       1, 4;
14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  std::cout << "a + b =\n" << a + b << std::endl;
15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  std::cout << "a - b =\n" << a - b << std::endl;
16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  std::cout << "Doing a += b;" << std::endl;
17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  a += b;
18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  std::cout << "Now a =\n" << a << std::endl;
19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  Vector3d v(1,2,3);
20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  Vector3d w(1,0,0);
21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath  std::cout << "-v + w - v =\n" << -v + w - v << std::endl;
22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}
23