pub fn rank_k_matrix(m: usize, n: usize, k: usize) -> DTensor<f64, 2>
Generate a rank-k matrix by multiplying m×k and k×n matrices