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