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

fix twisted selfadjoint to selfadjoint (conjugation issue)

Browse Source
This commit is contained in:
Gael Guennebaud 2011-01-27 14:39:01 +01:00
parent f5d0f115b4
commit 9ccd16609c
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

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