unsupported/doc/examples/MatrixExponential.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

#include <unsupported/Eigen/MatrixFunctions>

using namespace Eigen;

int main()
{
  const double pi = std::acos(-1.0);

  MatrixXd A(3,3);
  A << 0,    -pi/4, 0,
       pi/4, 0,     0,
       0,    0,     0;
  std::cout << "The matrix A is:\n" << A << "\n\n";

  MatrixXd B = ei_matrix_exponential(A);
  std::cout << "The matrix exponential of A is:\n" << B << "\n\n";
}