N-D array reduction by arbitrary function

Implements a fast reduction of N-D array that computes arbitrary function to sub-blocks of A

This function reduces a d-dimensional array into a d-dimensional array of smaller size by applying an arbitrary function to any sub-block extracted from A of given size. It includes testing scripts to demonstrate how it can be used.