Generalize reduction code so that it can work for higher order statistics.
Add code to compute Reynolds statistics.
Add support offsets in reduction kernels and the before_main_loop hook.