aboutsummaryrefslogtreecommitdiffhomepage
path: root/Eigen/src/Core/products/SelfadjointRank2Update.h
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2010-11-19 16:50:15 +0100
committerGravatar Gael Guennebaud <g.gael@free.fr>2010-11-19 16:50:15 +0100
commitf369b5a7116ae73405d66dee1f2c23f52d40e3d4 (patch)
tree9d02c839d135d3318d6dc0aee01c77b7570ea23a /Eigen/src/Core/products/SelfadjointRank2Update.h
parente14f14642d6101f94274459cfd499d2ec28925e4 (diff)
makes rank 2 update function conformant to BLAS HER2
Diffstat (limited to 'Eigen/src/Core/products/SelfadjointRank2Update.h')
-rw-r--r--Eigen/src/Core/products/SelfadjointRank2Update.h4
1 files 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);
}
};