pub fn matmul(a: &[f64], m: usize, k: usize, b: &[f64], n: usize) -> Vec<f64>
Compute C = A * B where A is m × k and B is k × n (all row-major). Returns a flat m*n vector.
C = A * B
A
m × k
B
k × n
m*n