1 2 3 4 5 6 7 8 9 10 11
pub fn gemm_f32(a: &[f32], b: &[f32], c: &mut [f32], m: usize, n: usize, k: usize) { for i in 0..m { for j in 0..n { let mut sum = 0.0; for l in 0..k { sum += a[i * k + l] * b[j * k + l]; } c[i * n + j] = sum; } } }