Struct MacroMicroRepr 

pub struct MacroMicroRepr {
    pub density: Vec<f64>,
    pub mean_velocity: Vec<f64>,
    pub temperature: Vec<f64>,
    pub inner: Box<dyn PhaseSpaceRepr>,
    pub spatial_shape: [usize; 3],
    pub velocity_shape: [usize; 3],
    pub domain: Domain,
    /* private fields */
}

Macro-micro representation: f = f_M + g.

f_M is a local Maxwellian parameterized by (rho, u, T). g is the kinetic deviation, stored by an inner PhaseSpaceRepr.

Density is O(N^3) from the macro fields – no velocity integration needed. Conservation of mass, momentum, and energy is exact by construction.

Fields

density: Vec<f64>

Macro density rho(x) – flat row-major [nx*ny*nz].

mean_velocity: Vec<f64>

Macro mean velocity u(x) – flat [nxnynz * 3] (ux, uy, uz interleaved).

temperature: Vec<f64>

Macro temperature T(x) – flat [nx*ny*nz] (scalar temperature).

inner: Box<dyn PhaseSpaceRepr>

Micro deviation g = f - f_M, stored by inner representation.

spatial_shape: [usize; 3]

Spatial grid shape [nx, ny, nz].

velocity_shape: [usize; 3]

Velocity grid shape [nv1, nv2, nv3].

domain: Domain

Domain specification.

Implementations

impl MacroMicroRepr

pub fn from_repr(inner: Box<dyn PhaseSpaceRepr>, domain: &Domain) -> Self

Create from an existing PhaseSpaceRepr by extracting macro moments and computing g = f - f_M.

pub fn reproject_moments(&mut self)

Re-project: recompute macro fields (rho, u, T) from the current inner representation state, ensuring exact conservation of macroscopic moments.

Trait Implementations

impl PhaseSpaceRepr for MacroMicroRepr

fn compute_density(&self) -> DensityField

Returns density directly from the cached macro field (O(N^3), no velocity integration).

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

Delegates spatial advection to the inner representation. Macro fields are re-synced via reproject_moments() after a full step.

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

Delegates velocity advection to the inner representation. Macro momentum update is applied via reproject_moments() after a full step.

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

Computes velocity moments by delegating to the inner kinetic representation.

fn total_mass(&self) -> f64

Returns the cached total mass from macro density integration.

fn casimir_c2(&self) -> f64

Approximates C2 = integral(f^2) using the fluid-level density: sum(rho^2 * dx^3).

fn entropy(&self) -> f64

Fluid-level entropy: configurational (-integral rho ln rho) plus thermal (3/2 rho ln(2 pi e T)).

fn stream_count(&self) -> StreamCountField

Delegates stream counting to the inner kinetic representation.

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

Delegates local velocity distribution extraction to the inner representation.

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

Kinetic energy from macro fields: bulk (1/2 rho |u|^2) plus thermal (3/2 rho T).

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

Produces a full 6D snapshot by delegating to the inner representation.

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

Loads a snapshot into the inner representation and re-syncs macro fields.

fn can_materialize(&self) -> bool

Reports whether the inner representation supports dense materialization.

fn memory_bytes(&self) -> usize

Total memory: macro fields (density + velocity + temperature) plus inner representation.

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

Attaches a shared progress reporter to both this wrapper and the inner representation.

fn as_any(&self) -> &dyn Any

Downcast to concrete type for implementation-specific queries (e.g. HT rank data).

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

Mutable downcast for in-place modification (e.g. BUG integrator leaf updates). Returns self as &mut dyn Any. Only HtTensor overrides this; the default returns self but concrete downcasts will yield None.