Function rgsl::linear_algebra::LU_svx

pub fn LU_svx(lu: &MatrixF64, p: &Permutation, x: &mut VectorF64) -> Value

This function solves the square system A x = b in-place using the precomputed LU decomposition of A into (LU,p). On input x should contain the right-hand side b, which is replaced by the solution on output.