Module ode_solver 

ODE Solver for coupled differential equations

Implements Runge-Kutta 4th order (RK4) solver for systems of ordinary differential equations, with GPU acceleration via PlatformGpuIntegrator.

Architecture

• Iterative Workload: Chains thousands of RK4 steps through the production queue
• Kahan Accumulation: f64 precision for error propagation control
• GPU Acceleration: Uses PlatformGpuIntegrator for batch step computation
• WAL Persistence: Each step persisted to NVMe WAL for fault tolerance

Usage

use qualia_core_db::modalities::calculus::ode_solver::{ExponentialDecay, Rk4Solver};

let system = ExponentialDecay::new(0.5);
let mut solver = Rk4Solver::new(system, 0.001);
let solution = solver.solve(0.0, 10.0, 1.0);

Structs

ChaoitonProfile

Chaoiton radial profile β(r) for astrophysics applications

CoupledBoltzmann

Coupled Boltzmann equations for step-size sensitivity analysis

ExponentialDecay

Exponential decay: dy/dt = -λy

HarmonicOscillator

Simple harmonic oscillator: d²x/dt² + ω²x = 0

LinearDecayBvp

Simple linear BVP for testing: dy/dt = -y

QuantizationMapper

Canonical quantization equivalence mapper

Rk4Solver

Runge-Kutta 4th order ODE solver with GPU acceleration

ShootingMethod

Boundary Value Problem solver using the Shooting Method

StandardModelMasses

Standard Model particle mass constants for validation

StepSizeAnalyzer

Step-size sensitivity analyzer for ODE solvers

Traits

BvpSystem

Represents a Boundary Value Problem system

OdeSystem

Represents a system of first-order ODEs: dy/dt = f(t, y)

Functions

create_ode_step_quin

Creates a Quin for an ODE solver step

extract_ode_state

Extracts ODE state from a Quin

pack_ode_state

Packs ODE state into a Quin