What is the recommended (i.e., fast) way to do
folds and scans across one axis of a NumPy array?
(The equivalent of Mma's Fold and FoldList.)
Assume an arbitrary Python function that
produces one array from two arrays, where these
three arrays share the same dimensions.
(So the use of Python's `reduce` is certainly
conceivable for the fold.)