1#include <unsupported/Eigen/MatrixFunctions>
2#include <iostream>
3
4using namespace Eigen;
5
6int main()
7{
8  using std::sqrt;
9  MatrixXd A(3,3);
10  A << 0.5*sqrt(2), -0.5*sqrt(2), 0,
11       0.5*sqrt(2),  0.5*sqrt(2), 0,
12       0,            0,           1;
13  std::cout << "The matrix A is:\n" << A << "\n\n";
14  std::cout << "The matrix logarithm of A is:\n" << A.log() << "\n";
15}
16