pub fn solve_linear_system<F>( a: &ArrayView2<'_, F>, b: &ArrayView1<'_, F>, ) -> IntegrateResult<Array1<F>>where F: Float + FromPrimitive + Debug + AddAssign + SubAssign + MulAssign,
Solve a linear system Ax = b using Gaussian elimination with partial pivoting
a
b
Result<Array1<F>, IntegrateError>