Parallel scan primitives compute element-wise inclusive or exclusive prefix sums of input vectors contributed by p consecutively ranked processors under an associative, possibly expensive, binary operator ⊕. In message-passing systems with bounded, one-ported communication capabilities, at least ⌈log2p⌉ or ⌈log2(p−1)⌉ send-receive communication rounds are required to perform the scans. While there are well-known, simple algorithms for the inclusive scan that solve the problem in ⌈log2p⌉ send-receive communication rounds with ⌈log2p⌉ applications of the ⊕ operator, the exclusive scan is different and has been much less addressed.
By considering natural invariants for the exclusive prefix sums problem, we present two different algorithms that are efficient in the number of communication rounds and in the number of applications of the ⊕ operator. The first algorithm consists of an inclusive scan phase and an exclusive scan phase and trades the number of communication rounds against the number of applications of the ⊕ operator. The smallest number of inclusive scan rounds with q=⌈log2p⌉ rounds in total is q′≥q−log2(2q−p+1). The other algorithm is a modification of a round-optimal all-reduce algorithm, and the number of additional applications of the ⊕ operator is dependent on the number of bits set (popcount of) in p−1. Both algorithms are relevant for small(er) input vectors where performance is dominated by the number of communication rounds. For large input vectors, other (pipelined, fixed-degree tree) algorithms must be used.Parallel scan primitives are included in MPI as the MPI Scan and MPI Exscan collectives.
An experimental comparison with implementations in MPI of the MPI Exscan collective
operation indicates that the considerations here are of practical relevance.
CC BY 4.0
