Struct PicSimulation 

pub struct PicSimulation {
    pub particles: Vec<PicParticle>,
    pub grid: PicGrid,
    pub dt: f32,
    pub time: f32,
}

Fields

particles: Vec<PicParticle>

grid: PicGrid

dt: f32

time: f32

Implementations

impl PicSimulation

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

pub fn deposit_charge(&mut self)

CIC (Cloud-In-Cell) charge deposition: distribute particle charge to grid.

pub fn solve_poisson(&mut self)

Solve Poisson equation: ∇²φ = -ρ/ε₀ using SOR (Successive Over-Relaxation).

pub fn compute_field(&mut self)

Compute electric field from potential: E = -∇φ (central differences).

pub fn push_particles(&mut self)

Push particles using Boris algorithm with interpolated fields.

pub fn step(&mut self)

Full PIC cycle: deposit → solve → field → push.

pub fn initialize_two_stream(&mut self, n_per_species: usize, v_beam: f32)

Initialize a two-stream instability: two counter-propagating beams.

pub fn initialize_landau_damping(&mut self, n: usize, k: f32, amplitude: f32)

Initialize a Landau damping test: uniform plasma with a small density perturbation.

pub fn kinetic_energy(&self) -> f32

Total kinetic energy of all particles.

pub fn field_energy(&self) -> f32

Total field energy.

pub fn total_energy(&self) -> f32

Total energy (kinetic + field).

pub fn total_charge(&self) -> f32

Total charge on the grid (should be conserved).