GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-06-04 18:54:00 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix typo in MatrixExponential noticed by Markos.

Browse Source
This commit is contained in:
Jitse Niesen 2014-07-25 13:34:03 +01:00
parent a0a87410d0
commit 5f3d542b8a
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -176,8 +176,8 @@ void matrix_exp_pade17(const MatrixType &A, MatrixType &U, MatrixType &V)
const MatrixType A4 = A2 * A2; const MatrixType A4 = A2 * A2;
const MatrixType A6 = A4 * A2; const MatrixType A6 = A4 * A2;
const MatrixType A8 = A4 * A4; const MatrixType A8 = A4 * A4;
V = b[17] * m_tmp1 + b[15] * A6 + b[13] * A4 + b[11] * A2; // used for temporary storage V = b[17] * A8 + b[15] * A6 + b[13] * A4 + b[11] * A2; // used for temporary storage
matrixType tmp = A8 * V; MatrixType tmp = A8 * V;
tmp += b[9] * A8 + b[7] * A6 + b[5] * A4 + b[3] * A2 tmp += b[9] * A8 + b[7] * A6 + b[5] * A4 + b[3] * A2
+ b[1] * MatrixType::Identity(A.rows(), A.cols()); + b[1] * MatrixType::Identity(A.rows(), A.cols());
U.noalias() = A * tmp; U.noalias() = A * tmp;