reduce-operations.hpp

#ifndef DIY_REDUCE_OPERATIONS_HPP
#define DIY_REDUCE_OPERATIONS_HPP

#include "reduce.hpp"
#include "partners/swap.hpp"
#include "detail/reduce/all-to-all.hpp"

namespace diy
{

template<class Op>
void
all_to_all(Master&              master,     
           const Assigner&      assigner,   
           const Op&            op,         
           int                  k = 2       
          )
{
  auto scoped = master.prof.scoped("all_to_all");
  RegularDecomposer<DiscreteBounds> decomposer(1, interval(0,assigner.nblocks()-1), assigner.nblocks());
  RegularSwapPartners  partners(decomposer, k, false);
  reduce(master, assigner, partners, detail::AllToAllReduce<Op>(op, assigner), detail::SkipIntermediate(partners.rounds()));
}

}

#endif