Struct FdtdGrid2D 

pub struct FdtdGrid2D {
    pub ez: Vec<f32>,
    pub hx: Vec<f32>,
    pub hy: Vec<f32>,
    pub nx: usize,
    pub ny: usize,
    pub dx: f32,
    pub dt: f32,
    pub permittivity: Vec<f32>,
    pub permeability: Vec<f32>,
    pub conductivity: Vec<f32>,
    /* private fields */
}

2D FDTD grid for TM (Transverse Magnetic) mode: Ez, Hx, Hy.

Fields

ez: Vec<f32>

hx: Vec<f32>

hy: Vec<f32>

nx: usize

ny: usize

dx: f32

dt: f32

permittivity: Vec<f32>

permeability: Vec<f32>

conductivity: Vec<f32>

Implementations

impl FdtdGrid2D

pub fn new(nx: usize, ny: usize, dx: f32, dt: f32) -> Self

pub fn step(&mut self)

Advance one timestep using the leapfrog scheme. Updates H first (half-step), then E (full step).

pub fn add_source(&mut self, x: usize, y: usize, value: f32)

Inject a soft source at point (x, y) into Ez.

pub fn field_energy(&self) -> f32

Total electromagnetic energy in the 2D grid.

pub fn courant_number(&self) -> f32

Courant number for 2D: dt * c / dx. Must be < 1/sqrt(2) for stability.

pub fn apply_pml_boundary(&mut self, layers: usize)

Apply Perfectly Matched Layer absorbing boundary conditions. layers is the number of PML cells on each side.

pub fn apply_mur_abc(&mut self)

Apply first-order Mur absorbing boundary condition.

pub fn set_material_region( &mut self, x0: usize, y0: usize, x1: usize, y1: usize, eps: f32, mu: f32, sigma: f32, )

Set material properties for a rectangular region.