pub fn linear_solve(matrix: &[Vec<f64>], rhs: &[f64]) -> Option<Vec<f64>>
小型矩阵使用显式主元和融合加乘消元，避免 BLAS 选择不同舍入顺序。