# Benjamin L’Huillier(벤자민)

Cosmologist and Astrophysicist, Research Fellow at KASI, Daejeon, South Korea

# ComputePk (ascl:1403.015)

If you use the code for your research, please cite it as L'Huillier (2014). You may use the following bibtex reference:

@MISC{2014ascl.soft03015L,
author = {{L'Huillier}, {B.}},
title = "{computePk: Power spectrum computation}",
note = {in Astrophysics Source Code Library, record ascl:1403.015},
year = 2014,
archivePrefix = "ascl",
eprint = {1403.015},
month = mar,
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}


## Description

Code to compute the density power spectrum in cosmological simulations. It uses MPI and FFTW-2.15. It is strongly inspired by N-GenIC for the parallelisation and from GenPK for the calculation of the power spectrum. So far it can read gadget format, but more will be available. It has been tested on $N_\text{g} = 2048^3$ simulation on 24 nodes and performed well. You may chose NGP, CIC or TSC for the mass assignment. I personally use $N_\text{g}^3 = 8\,N_\text{p}^3$ for a simulation using $N_\text{p}^3$ particles. Should you find any bug, please contact me.

The code has only be tested using Intel compiler (icc).

## Release notes

• Update 2016-06-13: the code was previously not accessible, it is now fixed
• 2015-02-03: v 2.1.0. Fixed bugs related to Gadget with baryons
• 2014-05-12: The code was upgraded to version 2.0. The main change is that it now can deal with baryons (Gadget simulations only). $\Omega_\text{b}$ should be defined in the header (cf allvars.h and the README).