aboutsummaryrefslogtreecommitdiffhomepage
path: root/blas/PackedSelfadjointProduct.h
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2013-06-12 17:56:15 +0200
committerGravatar Gael Guennebaud <g.gael@free.fr>2013-06-12 17:56:15 +0200
commitf75419c711b3df73ca1e9b5af94db28d437a0698 (patch)
tree3f136e0d246bd05c9eb7c074336a0299b18891ad /blas/PackedSelfadjointProduct.h
parentf3a029e9579c7290e1e6737f7239afac51ada097 (diff)
Add missing changes.
Diffstat (limited to 'blas/PackedSelfadjointProduct.h')
-rw-r--r--blas/PackedSelfadjointProduct.h5
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);
}
}