mirror/eigen
2
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-01-06 14:14:46 +08:00
Code Issues Packages Projects Releases Wiki Activity

Avoid division by zero in nonZerosEstimate() for empty blocks.

Browse Source
This commit is contained in:
Rasmus Munk Larsen 2020-07-22 01:38:30 +00:00
parent ac2eca6b11
commit acab22c205
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
Eigen/src/SparseCore/SparseBlock.h
View File

@ -446,9 +446,13 @@ struct unary_evaluator<Block<ArgType,BlockRows,BlockCols,InnerPanel>, IteratorBa
{} {}
inline Index nonZerosEstimate() const { inline Index nonZerosEstimate() const {
Index nnz = m_block.nonZeros(); const Index nnz = m_block.nonZeros();
if(nnz<0) if(nnz < 0) {
return m_argImpl.nonZerosEstimate() * m_block.size() / m_block.nestedExpression().size(); // Scale the non-zero estimate for the underlying expression linearly with block size.
// Return zero if the underlying block is empty.
const Index nested_sz = m_block.nestedExpression().size();
return nested_sz == 0 ? 0 : m_argImpl.nonZerosEstimate() * m_block.size() / nested_sz;
}
return nnz; return nnz;
} }