Support `OrdinalRange`s in `ndrange` argument

[MasonProtter]

I think it'd be nice if one could specify e.g.

kernel(u, v, ndrange=(2:2:size(u, 1), 2:2:size(u,2)))

to get kernel indices that don't necessarily start at 1 or iterate in unit steps. This'd also fit better with the name ndrange