dont use =* might not return a Scalar

This commit is contained in:
janos 2020-10-02 14:36:51 +02:00
parent 9078f47cd6
commit f7b185a8b1

View File

@ -136,7 +136,9 @@ EIGEN_DONT_INLINE void triangular_solve_matrix<Scalar,Index,OnTheLeft,Mode,Conju
} }
else else
{ {
Scalar b = (other(i,j) *= a); Scalar& otherij = other(i,j);
otherij *= a;
Scalar b = otherij;
typename OtherMapper::LinearMapper r = other.getLinearMapper(s,j); typename OtherMapper::LinearMapper r = other.getLinearMapper(s,j);
typename TriMapper::LinearMapper l = tri.getLinearMapper(s,i); typename TriMapper::LinearMapper l = tri.getLinearMapper(s,i);
for (Index i3=0;i3<rs;++i3) for (Index i3=0;i3<rs;++i3)