Struct StencilOperators 

pub struct StencilOperators { /* private fields */ }

Stencil operators

Implementations

impl StencilOperators

pub fn new() -> Self

pub fn register_operator(&mut self, name: &str, stencil: StencilOperator)

Register a named stencil operator.

pub fn register_boundary_stencil( &mut self, name: &str, stencil: BoundaryStencil, )

Register a named boundary stencil.

pub fn apply_derivative( &self, name: &str, field: &[f64], dx: f64, index: usize, ) -> Result<f64, PhysicsError>

Apply a registered stencil to compute the spatial derivative at index.

The derivative is computed as:

  sum_i( coefficients[i] * field[index + offset_i] ) / dx

where offset_i is taken from stencil_points[i].relative_position[0]. The coefficients are expected to already include the normalisation factor (e.g. [-0.5, 0.0, 0.5] for a 2nd-order central difference).

pub fn central_difference_2nd_order() -> StencilOperator

Create a 3-point 2nd-order central difference stencil.

Coefficients [-0.5, 0.0, 0.5] at offsets [-1, 0, +1] give (field[i+1] - field[i-1]) / (2*dx).

pub fn forward_difference_1st_order() -> StencilOperator

Create a 2-point 1st-order forward difference stencil.

Coefficients [-1.0, 1.0] at offsets [0, +1] give (field[i+1] - field[i]) / dx.