Struct SimulationGrid3D 

pub struct SimulationGrid3D {
    pub width: usize,
    pub height: usize,
    pub depth: usize,
    pub pressure: Vec<f32>,
    pub pressure_prev: Vec<f32>,
    pub cell_types: Vec<CellType>,
    pub reflection_coeff: Vec<f32>,
    pub params: AcousticParams3D,
    pub step: u64,
    pub time: f32,
    /* private fields */
}

3D simulation grid using Structure of Arrays (SoA) layout.

Memory layout is z-major (z varies slowest): index = z * (width * height) + y * width + x

Fields

width: usize

Grid dimensions

height: usize

depth: usize

pressure: Vec<f32>

Current pressure field (row-major, z-major ordering)

pressure_prev: Vec<f32>

Previous pressure field (for time-stepping)

cell_types: Vec<CellType>

Cell types (boundary conditions)

reflection_coeff: Vec<f32>

Reflection coefficients per cell (0 = full absorption, 1 = full reflection)

params: AcousticParams3D

Acoustic parameters

step: u64

Current simulation step

time: f32

Accumulated simulation time

Implementations

impl SimulationGrid3D

pub fn new( width: usize, height: usize, depth: usize, params: AcousticParams3D, ) -> Self

Create a new 3D simulation grid.

Arguments

• width - Number of cells in X direction
• height - Number of cells in Y direction
• depth - Number of cells in Z direction
• params - Acoustic simulation parameters

pub fn index(&self, x: usize, y: usize, z: usize) -> usize

Convert (x, y, z) coordinates to linear index.

pub fn coords(&self, idx: usize) -> (usize, usize, usize)

Convert linear index to (x, y, z) coordinates.

pub fn in_bounds(&self, x: i32, y: i32, z: i32) -> bool

Check if coordinates are within bounds.

pub fn inject_impulse(&mut self, x: usize, y: usize, z: usize, amplitude: f32)

Inject a pressure impulse at a specific position.

pub fn inject_impulse_at(&mut self, pos: Position3D, amplitude: f32)

Inject a pressure impulse at a 3D position (in meters).

pub fn inject_spherical_impulse( &mut self, center: Position3D, amplitude: f32, radius_cells: f32, )

Inject a spherical impulse (Gaussian distribution).

pub fn set_cell_type( &mut self, x: usize, y: usize, z: usize, cell_type: CellType, )

Set cell type at a specific location.

pub fn create_obstacle( &mut self, min: Position3D, max: Position3D, cell_type: CellType, )

Create a rectangular obstacle.

pub fn step_sequential(&mut self)

Perform one FDTD time step (CPU, sequential).

pub fn step_parallel(&mut self)

Perform one FDTD time step (CPU, parallel with Rayon).

pub fn sample_pressure(&self, pos: Position3D) -> f32

Get pressure value at a specific position (interpolated).

pub fn get_xy_slice(&self, z: usize) -> Vec<f32>

Get an XY slice of the pressure field at a given Z.

pub fn get_xz_slice(&self, y: usize) -> Vec<f32>

Get an XZ slice of the pressure field at a given Y.

pub fn get_yz_slice(&self, x: usize) -> Vec<f32>

Get a YZ slice of the pressure field at a given X.

pub fn total_energy(&self) -> f32

Get the total energy in the simulation.

pub fn max_pressure(&self) -> f32

Get the maximum absolute pressure.

pub fn reset(&mut self)

Reset the simulation to initial state.

pub fn dimensions(&self) -> (usize, usize, usize)

Get grid dimensions as (width, height, depth).

pub fn physical_size(&self) -> (f32, f32, f32)

Get the physical size of the simulation domain in meters.

Trait Implementations

impl From<&SimulationGrid3D> for GridParams

fn from(grid: &SimulationGrid3D) -> Self

Converts to this type from the input type.