fn cmat_inv(a: &[C64], d: usize) -> Option<Vec<C64>>
Complex matrix inverse via Gauss-Jordan, d×d row-major.