mirror/eigen
2
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2024-12-15 07:10:37 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix infinite loop in sparselu

Browse Source
This commit is contained in:
Gael Guennebaud 2014-02-14 23:10:16 +01:00
parent 0508af4287
commit cd606bbc94
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Eigen/src/SparseLU/SparseLU_Memory.h
View File

@ -70,7 +70,7 @@ Index SparseLUImpl<Scalar,Index>::expand(VectorType& vec, Index& length, Index
if(num_expansions == 0 || keep_prev) if(num_expansions == 0 || keep_prev)
new_len = length ; // First time allocate requested new_len = length ; // First time allocate requested
else else
new_len = Index(alpha * length); new_len = (std::max)(length+1,Index(alpha * length));
VectorType old_vec; // Temporary vector to hold the previous values VectorType old_vec; // Temporary vector to hold the previous values
if (nbElts > 0 ) if (nbElts > 0 )
@ -107,7 +107,7 @@ Index SparseLUImpl<Scalar,Index>::expand(VectorType& vec, Index& length, Index
do do
{ {
alpha = (alpha + 1)/2; alpha = (alpha + 1)/2;
new_len = Index(alpha * length); new_len = (std::max)(length+1,Index(alpha * length));
#ifdef EIGEN_EXCEPTIONS #ifdef EIGEN_EXCEPTIONS
try try
#endif #endif