fn mat_inv_small(a: &[f64], n: usize) -> Vec<f64>
Small matrix inverse (Gauss-Jordan) for nf x nf matrices.