Struct Transport 

pub struct Transport<const N: usize> {
    pub plus: [f64; N],
    pub minus: [f64; N],
    pub source: f64,
}

Per-dimension transport coefficients for a finite-difference operator.

The operator is decomposed as

L V = sum_i [ plus_i (V_{i,+} - V) + minus_i (V_{i,-} - V) ] + source

plus_i and minus_i are the forward and backward transport coefficients documented in docs/src/solver/solver_pde.md. source collects the purely local terms.

The transport coefficients and the source have clean semantics: transport goes into the operator and the source goes into the right-hand side. There is no hidden cancellation term.

Fields

plus: [f64; N]

Forward transport coefficient in each dimension.

minus: [f64; N]

Backward transport coefficient in each dimension.

source: f64

Purely local source term.

Implementations

impl<const N: usize> Transport<N>

pub fn new(plus: [f64; N], minus: [f64; N], source: f64) -> Self

Builds a transport bundle from raw coefficients.

Examples

use solver::numeric::finite_difference::discretization::{Transport, DimensionKind};
let t = Transport::<2>::new([0.5, 0.1], [0.25, 0.0], 1.0);
assert_eq!(t.plus, [0.5, 0.1]);
assert_eq!(t.minus, [0.25, 0.0]);
assert_eq!(t.source, 1.0);
let _kind = DimensionKind::Diffusion;

pub fn zero() -> Self

An empty transport bundle with a zero source.

Trait Implementations

impl<const N: usize> Clone for Transport<N>

fn clone(&self) -> Transport<N>

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl<const N: usize> Copy for Transport<N>

impl<const N: usize> Debug for Transport<N>

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

Formats the value using the given formatter.