diff options
Diffstat (limited to 'Eigen/src/misc/SparseSolve.h')
-rw-r--r-- | Eigen/src/misc/SparseSolve.h | 6 |
1 files changed, 4 insertions, 2 deletions
diff --git a/Eigen/src/misc/SparseSolve.h b/Eigen/src/misc/SparseSolve.h index 244bb8ec7..05caa9266 100644 --- a/Eigen/src/misc/SparseSolve.h +++ b/Eigen/src/misc/SparseSolve.h @@ -54,8 +54,10 @@ template<typename _DecompositionType, typename Rhs> struct sparse_solve_retval_b static const int NbColsAtOnce = 4; int rhsCols = m_rhs.cols(); int size = m_rhs.rows(); - Eigen::Matrix<DestScalar,Dynamic,Dynamic> tmp(size,rhsCols); - Eigen::Matrix<DestScalar,Dynamic,Dynamic> tmpX(size,rhsCols); + // the temporary matrices do not need more columns than NbColsAtOnce: + int tmpCols = (std::min)(rhsCols, NbColsAtOnce); + Eigen::Matrix<DestScalar,Dynamic,Dynamic> tmp(size,tmpCols); + Eigen::Matrix<DestScalar,Dynamic,Dynamic> tmpX(size,tmpCols); for(int k=0; k<rhsCols; k+=NbColsAtOnce) { int actualCols = std::min<int>(rhsCols-k, NbColsAtOnce); |