look into optimizing Stridedrolling

Paths to explore:

• optimize series slicing for functions which requires one (or more) pd.Series as input
• C-implementation (or using cython memoryview of arrays)
• numba

caveats; as we sort-of inject lambda-like functions, and support a large variability of array types, how feasible is this?