oscode: fast solutions of oscillatory ODEs in cosmology

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

• Fruzsina Agocs
• Wednesday 29 January 2020, 13:45-14:15
• Sackler Lecture Theatre, IoA.

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.

This talk is included in these lists:

• Combined External Astrophysics Talks DAMTP
• Cosmology lists
• Cosmology, Astrophysics and General Relativity
• Institute of Astronomy Seminars
• Institute of Astronomy Talk Lists
• Kavli Institute for Cosmology Talk Lists
• Priscilla
• Sackler Lecture Theatre, IoA

Note that ex-directory lists are not shown.