Struct LinearSystem

pub struct LinearSystem { /* private fields */ }

A linear system Ax = b with a matrix A and a right-hand side vector b.

Implementations

impl LinearSystem

pub fn new(mat: DMatrix<f64>, rhs: DVector<f64>) -> IterSolverResult<Self>

Create a new LinearSystem with a matrix A and a right-hand side vector b.

Arguments

• mat - The matrix A.
• rhs - The right-hand side vector b.

Errors

Returns an error if the matrix is not square or if the matrix and the right-hand side vector do not match.

Examples

use nalgebra::{DMatrix, DVector};
use iterative_solvers::LinearProblem;

let mat = DMatrix::from_row_slice(2, 2, &[4.0, 1.0, 1.0, 3.0]);
let rhs = DVector::from_vec(vec![1.0, 2.0]);
let problem = LinearProblem::new(mat, rhs).unwrap();

pub fn solve( &self, solver: impl IterativeSolver, ) -> IterSolverResult<DVector<f64>>

Solve the linear problem Ax = b using the given solver.

Arguments

• solver - The solver to use.

Errors

Returns an error if the solver fails to solve the linear problem.

Examples

use nalgebra::{DMatrix, DVector};
use iterative_solvers::{LinearProblem, CG};

let mat = DMatrix::from_row_slice(2, 2, &[4.0, 1.0, 1.0, 3.0]);
let rhs = DVector::from_vec(vec![1.0, 2.0]);
let problem = LinearProblem::new(mat, rhs).unwrap();
let solver = CG::new(&mat, &rhs, 1e-6).unwrap();
let solution = problem.solve(solver).unwrap();