Struct SphericalRepr 

pub struct SphericalRepr {
    pub data: Vec<f64>,
    pub shape: [usize; 3],
    pub r_range: (f64, f64),
    pub v_range: (f64, f64),
    pub l_range: (f64, f64),
    pub dr: f64,
    pub dv: f64,
    pub dl: f64,
    pub domain: Domain,
}

Spherically symmetric phase-space representation on a (r, v_r, L) grid.

Fields

data: Vec<f64>

Distribution function values: f[ir * nv * nl + iv * nl + il].

shape: [usize; 3]

Grid sizes [nr, nv, nl].

r_range: (f64, f64)

Radial grid: r ∈ [r_min, r_max].

v_range: (f64, f64)

Radial velocity grid: v_r ∈ [-v_max, v_max].

l_range: (f64, f64)

Angular momentum grid: L ∈ [0, L_max].

dr: f64

Cell spacings.

dv: f64

dl: f64

domain: Domain

Domain (for PhaseSpaceRepr trait compatibility).

Implementations

impl SphericalRepr

pub fn new( domain: Domain, nr: usize, nv: usize, nl: usize, r_max: f64, v_max: f64, l_max: f64, ) -> Self

Create a new spherical representation with zero-initialized f.

r_min is set to r_max / nr to avoid the r=0 coordinate singularity.

pub fn r_at(&self, ir: usize) -> f64

Get radial coordinate for cell index.

pub fn vr_at(&self, iv: usize) -> f64

Get radial velocity for cell index.

pub fn l_at(&self, il: usize) -> f64

Get angular momentum for cell index.

pub fn index(&self, ir: usize, iv: usize, il: usize) -> usize

Flat index into the data array for grid point (ir, iv, il).

Trait Implementations

impl PhaseSpaceRepr for SphericalRepr

fn compute_density(&self) -> DensityField

Compute spherically averaged density rho(r) by integrating f over v_r and L.

fn advect_x(&mut self, _displacement: &DisplacementField, dt: f64)

Radial drift sub-step: semi-Lagrangian shift along r with dr/dt = v_r.

fn advect_v(&mut self, acceleration: &AccelerationField, dt: f64)

Velocity kick sub-step: semi-Lagrangian shift along v_r with dv_r/dt = -dPhi/dr + L^2/r^3 (gravitational + centrifugal acceleration).

fn moment(&self, _position: &[f64; 3], _order: usize) -> Tensor

Velocity moment (stub – returns an empty tensor).

fn total_mass(&self) -> f64

Total mass by integrating f over the (r, v_r, L) grid with 4pir^2 Jacobian.

fn casimir_c2(&self) -> f64

Casimir C2 = integral of f^2 over the reduced phase space.

fn entropy(&self) -> f64

Entropy S = -integral of f ln f over the reduced phase space.

fn stream_count(&self) -> StreamCountField

Stream count (stub – returns zeros; not applicable in spherical coords).

fn velocity_distribution(&self, _position: &[f64; 3]) -> Vec<f64>

Local velocity distribution (stub – returns empty; use radial profiles instead).

fn total_kinetic_energy(&self) -> Option<f64>

Total kinetic energy T = 0.5 * integral of f*(v_r^2 + L^2/r^2) over phase space.

fn to_snapshot(&self, time: f64) -> Option<PhaseSpaceSnapshot>

Serialize the 3D grid into a PhaseSpaceSnapshot for checkpointing.

fn load_snapshot( &mut self, snap: PhaseSpaceSnapshot, ) -> Result<(), CausticError>

Restore the grid data from a previously saved snapshot.

fn as_any(&self) -> &dyn Any

Downcast to &dyn Any for runtime type queries.

fn as_any_mut(&mut self) -> &mut dyn Any

Downcast to &mut dyn Any for runtime type queries.

fn memory_bytes(&self) -> usize

Heap memory used by the distribution function data array.

fn set_progress(&mut self, _p: Arc<StepProgress>)

Attach shared progress state for intra-phase cell-level reporting. Implementations can use this to report progress via StepProgress::set_intra_progress.

fn can_materialize(&self) -> bool

Whether full 6D materialization fits in available memory. Default: true. Compressed representations (HT, TT) should override to check shape vs a memory threshold.