Struct Diagnostics 

pub struct Diagnostics {
    pub history: Vec<GlobalDiagnostics>,
}

Accumulates time series of GlobalDiagnostics.

Fields

history: Vec<GlobalDiagnostics>

Implementations

impl Diagnostics

pub fn compute( &mut self, repr: &dyn PhaseSpaceRepr, potential: &PotentialField, time: f64, dx3: f64, ) -> GlobalDiagnostics

Compute all global diagnostics from the current representation and potential. dx3 is the spatial cell volume dx1dx2dx3.

pub fn compute_with_density( &mut self, repr: &dyn PhaseSpaceRepr, density: &DensityField, potential: &PotentialField, time: f64, dx3: f64, ) -> GlobalDiagnostics

Compute all global diagnostics using a pre-computed density field.

This avoids redundant compute_density() calls when the caller already has the density (e.g. from a Poisson solve in the same step).

pub fn kinetic_energy(repr: &dyn PhaseSpaceRepr) -> f64

Total kinetic energy T = ½∫fv² dx³dv³.

Returns 0.0 if the representation does not support kinetic energy computation.

pub fn potential_energy( density: &DensityField, potential: &PotentialField, dx3: f64, ) -> f64

Total potential energy W = ½∫ρΦ dx³. dx3 is the spatial cell volume.

pub fn virial_ratio(t: f64, w: f64) -> f64

Virial ratio 2T/|W|. Equals 1.0 at equilibrium.

pub fn max_relative_drift<F: Fn(&GlobalDiagnostics) -> f64>( &self, quantity: F, ) -> f64

Maximum relative drift of a conserved quantity over the full history.

Returns |Q(t) - Q(0)| / |Q(0)| for the worst timestep, or 0 if history is empty.

pub fn conservation_summary(&self) -> ConservationSummary

Conservation summary: max relative drifts in mass, energy, and Casimir C₂.

pub fn density_profile(density: &DensityField, dx: [f64; 3]) -> Vec<(f64, f64)>

Spherically averaged density profile ρ(r) at current timestep. Bins density cells by radius from domain centre, returns (r_bin, ρ_avg) pairs.