1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <unsupported/Eigen/MatrixFunctions> 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <iostream> 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathusing namespace Eigen; 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathint main() 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath const double pi = std::acos(-1.0); 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MatrixXd A(2,2); 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A << cos(pi/3), -sin(pi/3), 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath sin(pi/3), cos(pi/3); 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::cout << "The matrix A is:\n" << A << "\n\n"; 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::cout << "The matrix square root of A is:\n" << A.sqrt() << "\n\n"; 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::cout << "The square of the last matrix is:\n" << A.sqrt() * A.sqrt() << "\n"; 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 17