1 #include <unsupported/Eigen/MatrixFunctions>
2 #include <iostream>
3 
4 using namespace Eigen;
5 
main()6 int 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