GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-10-10 15:16:31 +08:00
Code Issues Packages Projects Releases Wiki Activity

CREDIT Hauke Heibel, windows compatibility fixes in MatrixExponential

Browse Source
This commit is contained in:
Benoit Jacob 2009-05-07 18:11:49 +00:00
parent 7080282a6d
commit 159ab4a043
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h
View File

@ -25,6 +25,10 @@
#ifndef EIGEN_MATRIX_EXPONENTIAL #ifndef EIGEN_MATRIX_EXPONENTIAL
#define EIGEN_MATRIX_EXPONENTIAL #define EIGEN_MATRIX_EXPONENTIAL
#ifdef _MSC_VER
template <typename Scalar> Scalar log2(Scalar v) { return std::log(v)/std::log(Scalar(2)); }
#endif
/** Compute the matrix exponential. /** Compute the matrix exponential.
* *
* \param M matrix whose exponential is to be computed. * \param M matrix whose exponential is to be computed.
@ -131,7 +135,7 @@ void ei_matrix_exponential(const MatrixBase<Derived> &M, typename ei_plain_matri
squarings = std::max(0, (int)ceil(log2(l1norm / maxnorm))); squarings = std::max(0, (int)ceil(log2(l1norm / maxnorm)));
PlainMatrixType A, A2, A4, A6; PlainMatrixType A, A2, A4, A6;
A = M / pow(2, squarings); A = M / pow(Scalar(2), squarings);
A2 = (A * A).lazy(); A2 = (A * A).lazy();
A4 = (A2 * A2).lazy(); A4 = (A2 * A2).lazy();
A6 = (A4 * A2).lazy(); A6 = (A4 * A2).lazy();