diff --git a/Eigen/src/SparseCore/SparseDenseProduct.h b/Eigen/src/SparseCore/SparseDenseProduct.h
index a9084192e..08846b16a 100644
--- a/Eigen/src/SparseCore/SparseDenseProduct.h
+++ b/Eigen/src/SparseCore/SparseDenseProduct.h
@@ -177,7 +177,7 @@ struct sparse_time_dense_product_impl<SparseLhsType,DenseRhsType,DenseResType, R
       Index n = lhs.outerSize();
       for(Index j=0; j<n; ++j)
       {
-        typename Res::Scalar tmp(0);
+        typename Res::Scalar tmp(res.coeffRef(j,c));
         for(LhsInnerIterator it(lhs,j); it ;++it)
           tmp += it.value() * rhs.coeff(it.index(),c);
         res.coeffRef(j,c) = alpha * tmp;