Function coordinates 

pub fn coordinates(
    flat: usize,
    shape: &[usize],
) -> impl Iterator<Item = (usize, usize)>

Each axis’s coordinate of flat position flat, from the last axis to the first. flat must be within the shape.