For a while now, I had been working on a lightweight code to compute neutrino oscillation probabilities for arbitrary oscillation scenarios.
NuOscProbExact uses the method developed by Ohlsson and Snellman to compute probabilities exactly without the need of diagonalizing the Hamiltonian.
NuOscProbExact is fully written in Python 3.7, open source, and publicly available at this GitHub repository: github.com/mbustama/NuOscProbExact .
In the repository you will also find documentation and examples; the code itself is distributed with example files ready to run.
For details, see the accompanying paper: