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 sparse selfadjoint time dense such that the other triangular part is not used at all

Browse Source
This commit is contained in:
Gael Guennebaud 2011-07-26 09:02:41 +02:00
parent 97ac0fd192
commit 8fa7e92e77
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/Sparse/SparseSelfAdjointView.h
View File

@ -227,10 +227,14 @@ class SparseSelfAdjointTimeDenseProduct
for (Index j=0; j<m_lhs.outerSize(); ++j) for (Index j=0; j<m_lhs.outerSize(); ++j)
{ {
LhsInnerIterator i(m_lhs,j); LhsInnerIterator i(m_lhs,j);
if (ProcessSecondHalf && i && (i.index()==j)) if (ProcessSecondHalf)
{ {
dest.row(j) += i.value() * m_rhs.row(j); while (i.index()<j) ++i;
++i; if(i && i.index()==j)
{
dest.row(j) += i.value() * m_rhs.row(j);
++i;
}
} }
for(; (ProcessFirstHalf ? i && i.index() < j : i) ; ++i) for(; (ProcessFirstHalf ? i && i.index() < j : i) ; ++i)
{ {