Normalized to: Auld, T.
[1]
oai:arXiv.org:astro-ph/0608174 [pdf] - 84051
Fast cosmological parameter estimation using neural networks
Submitted: 2006-08-08, last modified: 2007-09-17
We present a method for accelerating the calculation of CMB power spectra,
matter power spectra and likelihood functions for use in cosmological parameter
estimation. The algorithm, called CosmoNet, is based on training a multilayer
perceptron neural network and shares all the advantages of the recently
released Pico algorithm of Fendt & Wandelt, but has several additional benefits
in terms of simplicity, computational speed, memory requirements and ease of
training. We demonstrate the capabilities of CosmoNet by computing CMB power
spectra over a box in the parameter space of flat \Lambda CDM models containing
the 3\sigma WMAP1 confidence region. We also use CosmoNet to compute the WMAP3
likelihood for flat \Lambda CDM models and show that marginalised posteriors on
parameters derived are very similar to those obtained using CAMB and the WMAP3
code. We find that the average error in the power spectra is typically 2-3% of
cosmic variance, and that CosmoNet is \sim 7 \times 10^4 faster than CAMB (for
flat models) and \sim 6 \times 10^6 times faster than the official WMAP3
likelihood code. CosmoNet and an interface to CosmoMC are publically available
at www.mrao.cam.ac.uk/software/cosmonet.
[2]
oai:arXiv.org:astro-ph/0703445 [pdf] - 90254
{\sc CosmoNet}: fast cosmological parameter estimation in non-flat
models using neural networks
Submitted: 2007-03-16
We present a further development of a method for accelerating the calculation
of CMB power spectra, matter power spectra and likelihood functions for use in
cosmological Bayesian inference. The algorithm, called {\sc CosmoNet}, is based
on training a multilayer perceptron neural network. We compute CMB power
spectra (up to $\ell=2000$) and matter transfer functions over a hypercube in
parameter space encompassing the $4\sigma$ confidence region of a selection of
CMB (WMAP + high resolution experiments) and large scale structure surveys (2dF
and SDSS). We work in the framework of a generic 7 parameter non-flat
cosmology. Additionally we use {\sc CosmoNet} to compute the WMAP 3-year, 2dF
and SDSS likelihoods over the same region. We find that the average error in
the power spectra is typically well below cosmic variance for spectra, and
experimental likelihoods calculated to within a fraction of a log unit. We
demonstrate that marginalised posteriors generated with {\sc CosmoNet} spectra
agree to within a few percent of those generated by {\sc CAMB} parallelised
over 4 CPUs, but are obtained 2-3 times faster on just a \emph{single}
processor. Furthermore posteriors generated directly via {\sc CosmoNet}
likelihoods can be obtained in less than 30 minutes on a single processor,
corresponding to a speed up of a factor of $\sim 32$. We also demonstrate the
capabilities of {\sc CosmoNet} by extending the CMB power spectra and matter
transfer function training to a more generic 10 parameter cosmological model,
including tensor modes, a varying equation of state of dark energy and massive
neutrinos. {\sc CosmoNet} and interfaces to both {\sc CosmoMC} and {\sc
Bayesys} are publically available at {\tt
www.mrao.cam.ac.uk/software/cosmonet}.