Struct GpuSphPressureSolver 

pub struct GpuSphPressureSolver {
    pub n_particles: usize,
    pub pressure: Vec<f64>,
    pub density: Vec<f64>,
    pub positions: Vec<[f64; 3]>,
    pub masses: Vec<f64>,
    pub smoothing_h: f64,
    pub rest_density: f64,
    pub stiffness: f64,
}

GPU SPH pressure solver holding all per-particle state.

Fields

n_particles: usize

Total number of simulated particles.

pressure: Vec<f64>

Per-particle pressure values (Pa).

density: Vec<f64>

Per-particle density values (kg/m³).

positions: Vec<[f64; 3]>

Per-particle positions [x, y, z] in metres.

masses: Vec<f64>

Per-particle masses (kg).

smoothing_h: f64

Smoothing length h (m).

rest_density: f64

Rest density ρ₀ (kg/m³).

stiffness: f64

Stiffness constant k for the linearised EOS.

Implementations

impl GpuSphPressureSolver

pub fn new(n: usize, h: f64, rho0: f64, k: f64) -> Self

Create a new solver with n particles, smoothing length h, rest density rho0, and stiffness k.

pub fn particle_count(&self) -> usize

Return the number of particles in this solver.

pub fn set_position(&mut self, i: usize, pos: [f64; 3])

Set the position of particle i.

pub fn set_mass(&mut self, i: usize, mass: f64)

Set the mass of particle i.

pub fn total_mass(&self) -> f64

Compute total mass of all particles.

pub fn density_error(&self) -> f64

Compute the maximum relative density error: max |ρᵢ − ρ₀| / ρ₀.

pub fn compute_stats(&self) -> GpuSphStats

Compute density statistics.

pub fn gpu_compute_density(&mut self)

Mock GPU density computation: ρᵢ = Σⱼ mⱼ · W_poly6(|rᵢ − rⱼ|, h).

pub fn gpu_compute_pressure(&mut self)

Mock GPU pressure computation: pᵢ = k · (ρᵢ − ρ₀).

pub fn gpu_pressure_force(&self, i: usize) -> [f64; 3]

Compute pressure force on particle i: F_press_i = −Σⱼ mⱼ (pᵢ/ρᵢ² + pⱼ/ρⱼ²) ∇W_spiky.

pub fn gpu_viscosity_force( &self, i: usize, velocities: &[[f64; 3]], mu: f64, ) -> [f64; 3]

Compute viscosity force on particle i: F_visc_i = μ Σⱼ mⱼ/ρⱼ (vⱼ − vᵢ) ∇²W_visc.

Trait Implementations

impl Clone for GpuSphPressureSolver

fn clone(&self) -> GpuSphPressureSolver

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for GpuSphPressureSolver

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

Formats the value using the given formatter.