University of Cambridge > > Institute of Astronomy Seminars > oscode: fast solutions of oscillatory ODEs in cosmology

oscode: fast solutions of oscillatory ODEs in cosmology

Add to your list(s) Download to your calendar using vCal

If you have a question about this talk, please contact Catrina Diener.

Bayesian analysis of the cosmic microwave background (CMB) allows us to constrain cosmological parameters and infer the physics of the early universe. The forward modelling phase of Bayesian inference involves calculation of the likelihood function of the model parameters, which requires one to numerically solve linear, ordinary differential equations with oscillatory solutions efficiently. I’ll present a novel numerical routine, oscode, with associated open-source software, for such purposes. The method makes use of the Wentzel—Kramers—Brillouin approximation used in quantum mechanics to rapidly traverse highly oscillatory regions of the solution, giving a significant speedup compared to Runge—Kutta-like methods. Its efficiency is demonstrated on physical examples of the Schrödinger equation and the Mukhanov—Sasaki equation; the latter is solved to generate primordial power spectra of curvature perturbations in different cosmological scenarios.

This talk is part of the Institute of Astronomy Seminars series.

Tell a friend about this talk:

This talk is included in these lists:

Note that ex-directory lists are not shown.


© 2006-2023, University of Cambridge. Contact Us | Help and Documentation | Privacy and Publicity