diff options
author | Gael Guennebaud <g.gael@free.fr> | 2013-06-12 17:56:15 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2013-06-12 17:56:15 +0200 |
commit | f75419c711b3df73ca1e9b5af94db28d437a0698 (patch) | |
tree | 3f136e0d246bd05c9eb7c074336a0299b18891ad /blas/PackedSelfadjointProduct.h | |
parent | f3a029e9579c7290e1e6737f7239afac51ada097 (diff) |
Add missing changes.
Diffstat (limited to 'blas/PackedSelfadjointProduct.h')
-rw-r--r-- | blas/PackedSelfadjointProduct.h | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/blas/PackedSelfadjointProduct.h b/blas/PackedSelfadjointProduct.h index f7c9b9341..07327a246 100644 --- a/blas/PackedSelfadjointProduct.h +++ b/blas/PackedSelfadjointProduct.h @@ -30,10 +30,9 @@ struct selfadjoint_packed_rank1_update<Scalar,Index,ColMajor,UpLo,ConjLhs,ConjRh for (Index i=0; i<size; ++i) { - Map<Matrix<Scalar,Dynamic,1> >(mat, UpLo==Lower ? size-i : (i+1)) - += alpha * cj(vec[i]) * ConjRhsType(OtherMap(vec+(UpLo==Lower ? i : 0), UpLo==Lower ? size-i : (i+1))); + Map<Matrix<Scalar,Dynamic,1> >(mat, UpLo==Lower ? size-i : (i+1)) += alpha * cj(vec[i]) * ConjRhsType(OtherMap(vec+(UpLo==Lower ? i : 0), UpLo==Lower ? size-i : (i+1))); //FIXME This should be handled outside. - mat[UpLo==Lower ? 0 : i] = real(mat[UpLo==Lower ? 0 : i]); + mat[UpLo==Lower ? 0 : i] = numext::real(mat[UpLo==Lower ? 0 : i]); mat += UpLo==Lower ? size-i : (i+1); } } |