Struct LoMaC 

pub struct LoMaC {
    pub kfvs: KfvsSolver,
    pub spatial_shape: [usize; 3],
    pub velocity_shape: [usize; 3],
    pub dv: [f64; 3],
    pub v_min: [f64; 3],
    pub active: bool,
}

LoMaC conservation manager.

Maintains the macroscopic state alongside the kinetic solution, and provides the projection step after each truncation.

Fields

kfvs: KfvsSolver

KFVS macroscopic solver for evolving (ρ, J, e).

spatial_shape: [usize; 3]

Spatial grid shape [nx, ny, nz].

velocity_shape: [usize; 3]

Velocity grid shape [nv1, nv2, nv3].

dv: [f64; 3]

Velocity cell spacings.

v_min: [f64; 3]

Minimum velocity coordinates.

active: bool

Whether LoMaC is active. When false, acts as passthrough.

Implementations

impl LoMaC

pub fn new( spatial_shape: [usize; 3], velocity_shape: [usize; 3], dx: [f64; 3], dv: [f64; 3], v_min: [f64; 3], ) -> Self

Create a new LoMaC manager.

pub fn initialize_from_kinetic(&mut self, f: &[f64])

Initialize from a kinetic distribution function.

Extracts macroscopic moments from f and initializes the KFVS solver.

pub fn advance_macroscopic(&mut self, dt: f64, acceleration: &[[f64; 3]])

Step 2: Advance macroscopic state by dt using KFVS.

Call this BEFORE advancing the kinetic equation.

pub fn project(&self, f_truncated: &[f64]) -> Vec<f64>

Step 4: Project the truncated kinetic solution to restore moments.

Call this AFTER advancing and truncating the kinetic equation. Returns the corrected distribution function with exact conservation.

pub fn apply( &mut self, dt: f64, acceleration: &[[f64; 3]], f_truncated: &[f64], ) -> Vec<f64>

Combined step: advance macroscopic, then project kinetic solution.

This is the main LoMaC entry point for use after each time step:

1. Advance KFVS: (ρ,J,e)^n → (ρ,J,e)^{n+1}
2. Project: f̃^{n+1} → f^{n+1} matching (ρ,J,e)^{n+1}

pub fn total_mass(&self) -> f64

Get the current total mass from the macroscopic solver.

pub fn total_momentum(&self) -> [f64; 3]

Get the current total momentum from the macroscopic solver.

pub fn total_energy(&self) -> f64

Get the current total energy from the macroscopic solver.

pub fn conservation_error(&self, f: &[f64]) -> (f64, f64, f64)

Conservation error: difference between kinetic and macroscopic moments.

Returns (|Δρ/ρ|, |ΔJ/J|, |Δe/e|) where Δ = kinetic − macroscopic.