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