Rasmus Munk LarsenandRasmus Munk Larsen
8a3d61c790
libeigen/eigen!2509 Closes #2868 Co-authored-by: Rasmus Munk Larsen <rmlarsen@gmail.com>
20 lines
513 B
C++
20 lines
513 B
C++
// SPDX-FileCopyrightText: The Eigen Authors
|
|
// SPDX-License-Identifier: MPL-2.0
|
|
|
|
#include <unsupported/Eigen/MatrixFunctions>
|
|
#include <iostream>
|
|
|
|
using namespace Eigen;
|
|
|
|
std::complex<double> expfn(std::complex<double> x, int) { return std::exp(x); }
|
|
|
|
int main() {
|
|
const double pi = std::acos(-1.0);
|
|
|
|
MatrixXd A(3, 3);
|
|
A << 0, -pi / 4, 0, pi / 4, 0, 0, 0, 0, 0;
|
|
|
|
std::cout << "The matrix A is:\n" << A << "\n\n";
|
|
std::cout << "The matrix exponential of A is:\n" << A.matrixFunction(expfn) << "\n\n";
|
|
}
|