Function gen_krange 

pub fn gen_krange<T>(k_mesh: &Array1<usize>) -> Array3<T>
where
    T: usefloat + Div<Output = T>,