Struct AxisOperator 

pub struct AxisOperator {
    pub axis: usize,
    pub sub: Vec<f64>,
    pub diag: Vec<f64>,
    pub sup: Vec<f64>,
}

A tridiagonal operator along one axis of a TensorGrid, with per-node coefficients: row i of A u reads sub[i] * u[i - stride] + diag[i] * u[i] + sup[i] * u[i + stride].

sub must be zero on the first plane of the axis and sup on the last (there is no neighbor there); boundary conditions are whatever the boundary rows encode — an all-zero row holds the boundary value fixed through both explicit application and implicit solves.

Fields

axis: usize

sub: Vec<f64>

diag: Vec<f64>

sup: Vec<f64>

Implementations

impl AxisOperator

pub fn zero(grid: &TensorGrid, axis: usize) -> Self

An all-zero operator along axis (a starting point to fill in).

pub fn apply(&self, grid: &TensorGrid, u: &[f64]) -> Vec<f64>

A u.

pub fn solve_shifted(&self, grid: &TensorGrid, c: f64, rhs: &[f64]) -> Vec<f64>

Solve (I - c A) x = rhs, line by line with the Thomas algorithm — the implicit stage of theta and ADI schemes.

Trait Implementations

impl Clone for AxisOperator

fn clone(&self) -> AxisOperator

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for AxisOperator

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.