fn solve_linear_system(
a: &Array2<f64>,
b: &Array1<f64>,
) -> Result<Vec<f64>, RError>Expand description
Solve a linear system A * x = b using Gaussian elimination with partial pivoting. A must be square (ncol x ncol), b must have length ncol.