From f369b5a7116ae73405d66dee1f2c23f52d40e3d4 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Fri, 19 Nov 2010 16:50:15 +0100
Subject: [PATCH] makes rank 2 update function conformant to BLAS HER2

---
 Eigen/src/Core/products/SelfadjointRank2Update.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/Eigen/src/Core/products/SelfadjointRank2Update.h b/Eigen/src/Core/products/SelfadjointRank2Update.h
index 41968a596..583537a8a 100644
--- a/Eigen/src/Core/products/SelfadjointRank2Update.h
+++ b/Eigen/src/Core/products/SelfadjointRank2Update.h
@@ -43,7 +43,7 @@ struct selfadjoint_rank2_update_selector<Scalar,Index,UType,VType,Lower>
     for (Index i=0; i<size; ++i)
     {
       Map<Matrix<Scalar,Dynamic,1> >(mat+stride*i+i, size-i) +=
-                        (alpha * conj(u.coeff(i))) * v.tail(size-i)
+                        (conj(alpha)  * conj(u.coeff(i))) * v.tail(size-i)
                       + (alpha * conj(v.coeff(i))) * u.tail(size-i);
     }
   }
@@ -57,7 +57,7 @@ struct selfadjoint_rank2_update_selector<Scalar,Index,UType,VType,Upper>
     const Index size = u.size();
     for (Index i=0; i<size; ++i)
       Map<Matrix<Scalar,Dynamic,1> >(mat+stride*i, i+1) +=
-                        (alpha * conj(u.coeff(i))) * v.head(i+1)
+                        (conj(alpha)  * conj(u.coeff(i))) * v.head(i+1)
                       + (alpha * conj(v.coeff(i))) * u.head(i+1);
   }
 };