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.) Thank you, Alan Isaac