Struct AcousticCoreInverse2d 

pub struct AcousticCoreInverse2d<R>
where
    R: CfdScalar + ConjugateScalar<Real = R>,
{ /* private fields */ }

Closed-form inverse of the 2-D constant-coefficient acoustic core A₀ = I − β·∇² on a periodic 2^lx × 2^ly grid, via ADI dimensional splitting: A₀⁻¹ ≈ (I − β·∂ₓ²)⁻¹·(I − β·∂ᵧ²)⁻¹, each factor the 1-D closed-form inverse acting along one axis. The splitting error is the O(β²·∂ₓ²∂ᵧ²) cross term; free-stream exactness is preserved exactly (each 1-D factor maps a uniform field to itself). The per-axis stiffness is s = β/Δx², β/Δy².

Implementations

impl<R> AcousticCoreInverse2d<R>

where R: CfdScalar + ConjugateScalar<Real = R>,

pub fn new( lx: usize, ly: usize, dx: R, dy: R, beta: R, trunc: Truncation<R>, ) -> Result<Self, PhysicsError>

Build the 2-D ADI inverse of A₀ = I − β·∇² on a periodic 2^lx × 2^ly grid (cell sizes dx, dy). The shift powers S±^{2^j} are the existing 1-D shifts lifted onto each axis of the serial x-then-y mode layout — no new operator.

Errors

PhysicsError::DimensionMismatch if lx == 0 or ly == 0; PhysicsError::NumericalInstability if β, dx, or dy is not finite and positive.

pub fn apply( &self, b: &CausalTensorTrain<R>, ) -> Result<CausalTensorTrain<R>, PhysicsError>

Apply the 2-D inverse A₀⁻¹·b = (I−β∂ₓ²)⁻¹·(I−β∂ᵧ²)⁻¹·b.

Errors

Propagates the per-axis apply / rounding errors.