pub fn solve_linear_system(
    a: &CsrMatrix<f64>,
    b: &DVector<f64>,
    method: Method,
    tol: f64,
    max_iter: i32,
    omega: f64
) -> Stat
Expand description

Solves the linear system ax=b and returns a Stat instance:

Arguments:

  • a: matrix
  • b: vector of constant terms
  • method: an istance of the enum Method
  • tol: the tolerance required to stop the routine
  • max_iter: the maximum number of iteration after wich the routine stops
  • omega: the relaxation factor, used only if method is either Jacobi or Gauss-Seidel