From 9ccd16609cc100c0749692e90b9ac353b3faf153 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Thu, 27 Jan 2011 14:39:01 +0100
Subject: [PATCH] fix twisted selfadjoint to selfadjoint (conjugation issue)

---
 Eigen/src/Sparse/SparseSelfAdjointView.h | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/Eigen/src/Sparse/SparseSelfAdjointView.h b/Eigen/src/Sparse/SparseSelfAdjointView.h
index 3749c604a..0f76ed581 100644
--- a/Eigen/src/Sparse/SparseSelfAdjointView.h
+++ b/Eigen/src/Sparse/SparseSelfAdjointView.h
@@ -367,6 +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;
   
   Index size = mat.rows();
   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 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] = it.value();
+      dest._valuePtr()[k] = cj(it.value());
     }
   }
 }