fn cmat_mul(a: &[C64], b: &[C64], d: usize) -> Vec<C64>
Complex matrix multiply: C = A * B, all d×d row-major.