GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-14 04:35:57 +08:00
Code Issues Packages Projects Releases Wiki Activity

really fix permute_symm_to_symm for sparse complex matrix

Browse Source
This commit is contained in:
Gael Guennebaud 2011-01-28 15:51:55 +01:00
parent 6ec660ca7e
commit 817d86cbaf
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Eigen/src/Sparse/SparseSelfAdjointView.h
View File

@ -367,7 +367,7 @@ void permute_symm_to_symm(const MatrixType& mat, SparseMatrix<typename MatrixTyp
typedef SparseMatrix<Scalar,DestOrder,Index> Dest;
Dest& dest(_dest.derived());
typedef Matrix<Index,Dynamic,1> VectorI;
internal::conj_if<SrcUpLo!=DstUpLo> cj;
//internal::conj_if<SrcUpLo!=DstUpLo> cj;
Index size = mat.rows();
VectorI count(size);
@ -405,7 +405,11 @@ void permute_symm_to_symm(const MatrixType& mat, SparseMatrix<typename MatrixTyp
Index ip = perm? perm[i] : i;
Index k = count[DstUpLo==Lower ? std::min(ip,jp) : std::max(ip,jp)]++;
dest._innerIndexPtr()[k] = DstUpLo==Lower ? std::max(ip,jp) : std::min(ip,jp);
dest._valuePtr()[k] = cj(it.value());
if((DstUpLo==Lower && ip<jp) || (DstUpLo==Upper && ip>jp))
dest._valuePtr()[k] = conj(it.value());
else
dest._valuePtr()[k] = it.value();
}
}
}