GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-12 11:49:02 +08:00
Code Issues Packages Projects Releases Wiki Activity

Bug fix in MatrixExponential.h

Browse Source 
Initialize matrices for intermediate results to correct dimension
This commit is contained in:
Jitse Niesen 2009-09-01 10:50:54 +01:00
parent 6e4e94ff32
commit 32f95ec267
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -292,7 +292,10 @@ namespace MatrixExponentialInternal {
template <typename MatrixType> template <typename MatrixType>
void compute(const MatrixType &M, MatrixType* result) void compute(const MatrixType &M, MatrixType* result)
{ {
MatrixType num, den, U, V; MatrixType num(M.rows(), M.cols());
MatrixType den(M.rows(), M.cols());
MatrixType U(M.rows(), M.cols());
MatrixType V(M.rows(), M.cols());
MatrixType Id = MatrixType::Identity(M.rows(), M.cols()); MatrixType Id = MatrixType::Identity(M.rows(), M.cols());
float l1norm = static_cast<float>(M.cwise().abs().colwise().sum().maxCoeff()); float l1norm = static_cast<float>(M.cwise().abs().colwise().sum().maxCoeff());
int squarings; int squarings;