Function gemm 

pub fn gemm(
    a: &[f128],
    b: &[f128],
    c: &mut [f128],
    m: usize,
    n: usize,
    k: usize,
)

C = A * B, row-major. A is m×k, B is k×n, C is m×n.