pub fn lu_decompose(n: usize, data: &[f64]) -> Result<Lu, SolversError>Expand description
LU-decompose a row-major n×n matrix with partial pivoting. The reusable primitive
behind determinant (and a building block for solves / condition estimates). O(n³).
Returns SolversError::InvalidDimension for an empty or non-square input.