stiff-odes

KineticsLab: implicit solver for the stiff 6-group reactor point-kinetics equations, validated to 6.6e-10 vs a scipy oracle and closed-form inhour solutions, plus protection trip logic with 0 false/0 missed trips on 120 labeled transients.

A research repo exploring how adaptive linear approximation behaves near thin stiff layers. It proposes a threshold law for when layers control node placement, introduces a diagnostic fingerprint for subcritical inversion, and backs both with derivations, simulations, plots, and practical Python examples.

Divide-and-conquer Sobolev gradient solver for stiff ODEs — the original 2005 C++ preserved, plus a modern C++20 port that reproduces the published results

MATLAB benchmark suite for comparing numerical integrators on classical differential-equation problems, including linear decay, harmonic oscillator, Van der Pol, Lorenz, Robertson kinetics, and Kepler two-body dynamics. The repository includes performance plots, accuracy metrics, conservation analysis, and a simulation video available on YouTube.