Enum IncidentField 

pub enum IncidentField {
    PlaneWave {
        direction: [f64; 3],
        amplitude: Complex64,
    },
    PointSource {
        position: [f64; 3],
        strength: Complex64,
    },
    MultiplePlaneWaves(Vec<([f64; 3], Complex64)>),
    MultiplePointSources(Vec<([f64; 3], Complex64)>),
}

Incident field source type

Variants

PlaneWave

Plane wave: p = A * exp(i k·x)

Fields

direction: [f64; 3]

Direction of propagation (unit vector)

amplitude: Complex64

Complex amplitude

PointSource

Point source (monopole): p = A * exp(ikr) / (4πr)

Fields

position: [f64; 3]

Source position

strength: Complex64

Source strength (volume velocity amplitude)

MultiplePlaneWaves(Vec<([f64; 3], Complex64)>)

Multiple plane waves

MultiplePointSources(Vec<([f64; 3], Complex64)>)

Multiple point sources

Implementations

impl IncidentField

pub fn plane_wave_z() -> Self

Create a plane wave with unit amplitude traveling in +z direction

This matches the standard convention for Mie scattering theory where the incident wave travels toward +z (positive z direction).

pub fn plane_wave_neg_z() -> Self

Create a plane wave with unit amplitude traveling in -z direction

pub fn plane_wave(direction: [f64; 3], amplitude: f64) -> Self

Create a plane wave with specified direction and amplitude

pub fn point_source(position: [f64; 3], strength: f64) -> Self

Create a point source at given position

pub fn evaluate_pressure( &self, points: &Array2<f64>, physics: &PhysicsParams, ) -> Array1<Complex64>

Evaluate incident pressure at given points

Arguments

• points - Evaluation points (N × 3 array)
• physics - Physical parameters (contains wave number k)

Returns

Complex pressure values at each point

pub fn evaluate_normal_derivative( &self, points: &Array2<f64>, normals: &Array2<f64>, physics: &PhysicsParams, ) -> Array1<Complex64>

Evaluate incident velocity (∂p/∂n) at given points with normals

Arguments

• points - Evaluation points (N × 3 array)
• normals - Unit normal vectors at each point (N × 3 array)
• physics - Physical parameters

Returns

Complex normal velocity values (actually ∂p/∂n, not v_n)

pub fn compute_rhs( &self, element_centers: &Array2<f64>, element_normals: &Array2<f64>, physics: &PhysicsParams, use_burton_miller: bool, ) -> Array1<Complex64>

Compute the right-hand side vector for BEM

For exterior Neumann problem (rigid scatterer) using DIRECT formulation where the unknown is the total surface pressure p:

The RHS comes from the incident field contribution to the integral equation. From NumCalc NC_IncidentWaveRHS: RHS = -(γp_inc + βτ*∂p_inc/∂n)

For exterior problems: τ = +1, γ = 1 For interior problems: τ = -1

pub fn compute_rhs_with_beta( &self, element_centers: &Array2<f64>, element_normals: &Array2<f64>, physics: &PhysicsParams, beta: Complex64, ) -> Array1<Complex64>

Compute RHS with custom Burton-Miller coupling parameter

From NumCalc NC_IncidentWaveRHS: RHS = -(γp_inc + βτ*∂p_inc/∂n)

This formula matches the C++ implementation for proper BEM assembly.

Trait Implementations

impl Clone for IncidentField

fn clone(&self) -> IncidentField

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for IncidentField

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

Formats the value using the given formatter.