From 817d86cbaf04a6f92a7945b8cd5f3a8865c2dfbd Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Fri, 28 Jan 2011 15:51:55 +0100
Subject: [PATCH] really fix permute_symm_to_symm for sparse complex matrix

---
 Eigen/src/Sparse/SparseSelfAdjointView.h | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/Eigen/src/Sparse/SparseSelfAdjointView.h b/Eigen/src/Sparse/SparseSelfAdjointView.h
index 0f76ed581..8be358db8 100644
--- a/Eigen/src/Sparse/SparseSelfAdjointView.h
+++ b/Eigen/src/Sparse/SparseSelfAdjointView.h
@@ -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();
     }
   }
 }