25 #ifndef EIGEN_SELFADJOINTRANK2UPTADE_H
26 #define EIGEN_SELFADJOINTRANK2UPTADE_H
36 template<
typename Scalar,
typename Index,
typename UType,
typename VType,
int UpLo>
37 struct selfadjoint_rank2_update_selector;
39 template<
typename Scalar,
typename Index,
typename UType,
typename VType>
40 struct selfadjoint_rank2_update_selector<Scalar,Index,UType,VType,
Lower>
42 static void run(Scalar* mat, Index stride,
const UType& u,
const VType& v, Scalar alpha)
44 const Index size = u.size();
45 for (Index i=0; i<size; ++i)
47 Map<Matrix<Scalar,Dynamic,1> >(mat+stride*i+i, size-i) +=
48 (
conj(alpha) *
conj(u.coeff(i))) * v.tail(size-i)
49 + (alpha *
conj(v.coeff(i))) * u.tail(size-i);
54 template<
typename Scalar,
typename Index,
typename UType,
typename VType>
55 struct selfadjoint_rank2_update_selector<Scalar,Index,UType,VType,
Upper>
57 static void run(Scalar* mat, Index stride,
const UType& u,
const VType& v, Scalar alpha)
59 const Index size = u.size();
60 for (Index i=0; i<size; ++i)
61 Map<Matrix<Scalar,Dynamic,1> >(mat+stride*i, i+1) +=
62 (
conj(alpha) *
conj(u.coeff(i))) * v.head(i+1)
63 + (alpha *
conj(v.coeff(i))) * u.head(i+1);
67 template<
bool Cond,
typename T>
struct conj_expr_if
68 : conditional<!Cond, const T&,
69 CwiseUnaryOp<scalar_conjugate_op<typename traits<T>::Scalar>,T> > {};
73 template<
typename MatrixType,
unsigned int UpLo>
74 template<
typename DerivedU,
typename DerivedV>
78 typedef internal::blas_traits<DerivedU> UBlasTraits;
79 typedef typename UBlasTraits::DirectLinearAccessType ActualUType;
80 typedef typename internal::remove_all<ActualUType>::type _ActualUType;
81 typename internal::add_const_on_value_type<ActualUType>::type actualU = UBlasTraits::extract(u.derived());
83 typedef internal::blas_traits<DerivedV> VBlasTraits;
84 typedef typename VBlasTraits::DirectLinearAccessType ActualVType;
85 typedef typename internal::remove_all<ActualVType>::type _ActualVType;
86 typename internal::add_const_on_value_type<ActualVType>::type actualV = VBlasTraits::extract(v.derived());
91 enum { IsRowMajor = (internal::traits<MatrixType>::Flags&
RowMajorBit) ? 1 : 0 };
92 Scalar actualAlpha = alpha * UBlasTraits::extractScalarFactor(u.derived())
97 internal::selfadjoint_rank2_update_selector<
Scalar,
Index,
98 typename internal::remove_all<typename internal::conj_expr_if<IsRowMajor ^ UBlasTraits::NeedToConjugate,_ActualUType>::type>::type,
99 typename internal::remove_all<typename internal::conj_expr_if<IsRowMajor ^ VBlasTraits::NeedToConjugate,_ActualVType>::type>::type,
101 ::run(_expression().const_cast_derived().data(),_expression().outerStride(),actualU,actualV,actualAlpha);
108 #endif // EIGEN_SELFADJOINTRANK2UPTADE_H