Function gen_kmesh 

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