Struct CoupledSolver 

pub struct CoupledSolver {
    pub tolerance: f64,
    pub max_iter: usize,
    pub relaxation: f64,
}

Multiphysics coupled solver strategies.

Provides operator splitting, staggered scheme, and monolithic Newton iteration.

Fields

tolerance: f64

Convergence tolerance.

max_iter: usize

Maximum iterations.

relaxation: f64

Relaxation parameter ω for successive substitution (0 < ω ≤ 1).

Implementations

impl CoupledSolver

pub fn new(tolerance: f64, max_iter: usize, relaxation: f64) -> Self

Creates a new coupled solver.

pub fn operator_split<FA, FB>( &self, a0: f64, b0: f64, solve_a: FA, solve_b: FB, ) -> (f64, f64, usize)

where FA: Fn(f64) -> f64, FB: Fn(f64) -> f64,

Operator splitting: alternately solve field A then field B.

solve_a(b) -> a and solve_b(a) -> b are the individual field solvers. Returns the converged (a, b) pair and number of iterations.

pub fn staggered<F>(&self, state0: &[f64], step_fn: F) -> (Vec<f64>, usize)

where F: Fn(&[f64]) -> Vec<f64>,

Staggered scheme: sequential solve with fixed-point iteration.

Returns converged state vector and iteration count.

pub fn monolithic_newton<F>(&self, x0: &[f64], residual: F) -> (Vec<f64>, usize)

where F: Fn(&[f64]) -> Vec<f64>,

Monolithic Newton iteration for a coupled residual R(x) = 0.

Uses finite-difference Jacobian and direct solve via Gaussian elimination.