Author: | Svetlin Tassev |
---|---|
Version: | |release| |
Date: | June 26, 2014 |
Homepage: | pyCOLA Homepage |
Documentation: | PDF Documentation |
License: | GPLv3+ License |
pyCOLA is a multithreaded Python/Cython N-body code, implementing the Comoving Lagrangian Acceleration (COLA) method in the temporal and spatial domains. pyCOLA also implements a novel method to compute second-order cosmological initial conditions for given initial conditions at first-order for arbitrary initial particle configurations (including glass initial conditions, as well as initial conditions having refined subregions).
pyCOLA is based on the following two papers: [temporalCOLA], [spatialCOLA]. We kindly ask you [1] to acknowledge them and their authors in any program or publication in which you use the COLA method in the temporal and/or spatial domains.
The new method for calculating second-order cosmological initial conditions is based on the following paper: (todo: Daniel, let me know what to cite). Again, we kindly ask you to acknowledge that paper and its authors in any program or publication in which you use that method.
Footnotes
[1] | We do not require you, however, as we want pyCOLA to be GPLv3 compatible. |
pyCOLA requires NumPy, SciPy, pyFFTW, h5py, as well as their respective dependencies. Note that pyFFTW v0.9.2 does not support large arrays, so one needs to install the development version from github, where the bug has been fixed.
Note
All lengthscales are in units of comoving \mathrm{Mpc}/h, unless otherwise specified.
.. todo:: If there is interest in the code (i.e. not only the algorithm), it should be converted to use classes as that will enormously reduce the amount of arguments to be passed around, will make the code more readable, and reduce the chances for introducing bugs. Some of the functions are already converted to using classes in a separate branch, but converting the whole code will take some time.
[MUSIC] | Multi-scale initial conditions for cosmological simulations, O. Hahn, T. Abel, Monthly Notices of the Royal Astronomical Society, 415, 2101 (2011), arXiv:1103.6031. The code can be found on this website. |
[spatialCOLA] | Extending the N-body Comoving Lagrangian Acceleration Method to the Spatial Domain, S. Tassev, D. J. Eisenstein, B. D. Wandelt, M. Zaldarriaga, (2014), arXiv:14??.???? |
[temporalCOLA] | Solving Large Scale Structure in Ten Easy Steps with COLA, S. Tassev, M. Zaldarriaga, D. J. Eisenstein, Journal of Cosmology and Astroparticle Physics, 06, 036 (2013), arXiv:1301.0322 |