oscode: Oscillatory ordinary differential equation solver
Abstract
oscode solves oscillatory ordinary differential equations efficiently. It is designed to deal with equations of the form d^{2}x/dt^{2}+2γ(t)xË(t)+ω2(t)x(t)=0, where γ(t) and ω(t) can be given as explicit functions or sequence containers (Eigen::Vectors, arrays, std::vectors, lists) in C++ or as numpy.arrays in Python. oscode makes use of an analytic approximation of x(t) embedded in a stepping procedure to skip over long regions of oscillations, giving a reduction in computing time. The approximation is valid when the frequency changes slowly relative to the timescales of integration, it is therefore worth applying when this condition holds for at least some part of the integration range.
 Publication:

Astrophysics Source Code Library
 Pub Date:
 August 2019
 Bibcode:
 2019ascl.soft08012A
 Keywords:

 Software