diff options
author | Gael Guennebaud <g.gael@free.fr> | 2010-11-19 16:50:15 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2010-11-19 16:50:15 +0100 |
commit | f369b5a7116ae73405d66dee1f2c23f52d40e3d4 (patch) | |
tree | 9d02c839d135d3318d6dc0aee01c77b7570ea23a /Eigen | |
parent | e14f14642d6101f94274459cfd499d2ec28925e4 (diff) |
makes rank 2 update function conformant to BLAS HER2
Diffstat (limited to 'Eigen')
-rw-r--r-- | Eigen/src/Core/products/SelfadjointRank2Update.h | 4 |
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); } }; |