xref: /aosp_15_r20/external/eigen/unsupported/doc/examples/MatrixPower.cpp (revision bf2c37156dfe67e5dfebd6d394bad8b2ab5804d4)
1*bf2c3715SXin Li #include <unsupported/Eigen/MatrixFunctions>
2*bf2c3715SXin Li #include <iostream>
3*bf2c3715SXin Li 
4*bf2c3715SXin Li using namespace Eigen;
5*bf2c3715SXin Li 
main()6*bf2c3715SXin Li int main()
7*bf2c3715SXin Li {
8*bf2c3715SXin Li   const double pi = std::acos(-1.0);
9*bf2c3715SXin Li   Matrix3d A;
10*bf2c3715SXin Li   A << cos(1), -sin(1), 0,
11*bf2c3715SXin Li        sin(1),  cos(1), 0,
12*bf2c3715SXin Li 	   0 ,      0 , 1;
13*bf2c3715SXin Li   std::cout << "The matrix A is:\n" << A << "\n\n"
14*bf2c3715SXin Li 	       "The matrix power A^(pi/4) is:\n" << A.pow(pi/4) << std::endl;
15*bf2c3715SXin Li   return 0;
16*bf2c3715SXin Li }
17