Function decimate 

pub fn decimate(
    x: &[f32],
    dims: &[usize],
    stride: usize,
) -> (Vec<f32>, Vec<usize>)

Every-stride-th sample along each axis longer than 4·stride (row-major), with the new shape.