pub fn lu_solve(n: usize, a: &[f64], b: &[f64]) -> Option<Vec<f64>>Expand description
Solve a general row-major n×n system A x = b via LU with partial pivoting. None
on a shape mismatch or a singular matrix. The canonical dense solve for the engine.