Skip to content

Commit 11d96a5

Browse files
[opt] Avoid using neutral element during accumulation
1 parent b525707 commit 11d96a5

1 file changed

Lines changed: 2 additions & 2 deletions

File tree

pythran/pythonic/numpy/reduce.hpp

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -269,8 +269,8 @@ namespace numpy
269269
void operator()(E &&e, F &&f, long axis, EIndices &&e_indices)
270270
{
271271
if (axis == std::decay_t<E>::value - 1) {
272-
auto acc = utils::neutral<Op, typename std::decay_t<F>::dtype>::value;
273-
for (long i = 0, n = e.template shape<std::decay_t<E>::value - 1>(); i < n; ++i) {
272+
typename std::decay_t<F>::dtype acc = detail::loader(e, e_indices, 0);
273+
for (long i = 1, n = e.template shape<std::decay_t<E>::value - 1>(); i < n; ++i) {
274274
Op{}(acc, detail::loader(e, e_indices, i));
275275
}
276276
detail::storer(acc, std::forward<F>(f), std::forward<EIndices>(e_indices));

0 commit comments

Comments
 (0)