Efficient numerical integration of neutrino oscillations in matter
For neutrino physics researchers, this method accelerates simulations needed for experimental data analysis, but it is an incremental improvement in numerical methods.
The authors propose a Magnus expansion-based solver for three-neutrino oscillations in matter, achieving up to 100x speedup over general integrators, which could enable massive numerical integration for experimental data analysis.
A special purpose solver, based on the Magnus expansion, well suited for the integration of the linear three neutrino oscillations equations in matter is proposed. The computations are speeded up to two orders of magnitude with respect to a general numerical integrator, a fact that could smooth the way for massive numerical integration concomitant with experimental data analyses. Detailed illustrations about numerical procedure and computer time costs are provided.