GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-07-26 06:44:27 +08:00
Code Issues Packages Projects Releases Wiki Activity

makes rank 2 update function conformant to BLAS HER2

Browse Source
This commit is contained in:
Gael Guennebaud 2010-11-19 16:50:15 +01:00
parent e14f14642d
commit f369b5a711
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Eigen/src/Core/products/SelfadjointRank2Update.h
View File

@ -43,7 +43,7 @@ struct selfadjoint_rank2_update_selector<Scalar,Index,UType,VType,Lower>
for (Index i=0; i<size; ++i) for (Index i=0; i<size; ++i)
{ {
Map<Matrix<Scalar,Dynamic,1> >(mat+stride*i+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); + (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(); const Index size = u.size();
for (Index i=0; i<size; ++i) for (Index i=0; i<size; ++i)
Map<Matrix<Scalar,Dynamic,1> >(mat+stride*i, i+1) += 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); + (alpha * conj(v.coeff(i))) * u.head(i+1);
} }
}; };