Hamevol1.0: a C++ code for differential equations based on Runge-Kutta algorithm. An application to matter enhanced neutrino oscillation

TL;DR

Hamevol1.0 is a C++ software implementing a fifth order semi-implicit Runge-Kutta algorithm, efficiently solving differential equations, specifically applied to neutrino oscillation probabilities in matter, demonstrating competitive performance and accuracy.

Contribution

The paper introduces Hamevol1.0, a novel C++ implementation of a high-order Runge-Kutta algorithm tailored for differential equations, with specific application to neutrino physics.

Findings

Competitive performance with existing methods

High accuracy in neutrino oscillation calculations

Versatile application to various differential equations

Abstract

We present a C++ implementation of a fifth order semi-implicit Runge-Kutta algorithm for solving Ordinary Differential Equations. This algorithm can be used for studying many different problems and in particular it can be applied for computing the evolution of any system whose Hamiltonian is known. We consider in particular the problem of calculating the neutrino oscillation probabilities in presence of matter interactions. The time performance and the accuracy of this implementation is competitive with respect to the other analytical and numerical techniques used in literature. The algorithm design and the salient features of the code are presented and discussed and some explicit examples of code application are given.