-
Notifications
You must be signed in to change notification settings - Fork 7
/
README
47 lines (28 loc) · 1.21 KB
/
README
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
==GenPK==
This is GenPK, a program to calculate the 3D matter power spectrum.
To build, first make sure FFTW3 is installed on your system
and that the headers are in your PATH.
To build, type:
make
For the test suite (requires Boost::Test > 1.34), do:
make test
To run, do:
./gen-pk -i filename -o outputdir
filename can be any Gadget snapshot, with or without the .o, in format I or II.
For the doxygen-generated documentation, do:
make doc
To delete compiler output:
make clean
Further documentation may be found in the doxygen generated files in doc/html/index.html
(type "make doc" to make this file), by typing ./gen-pk --help or in gen-pk.cpp
Powerspectrum output is in Fourier units; to convert to, eg, Mpc/h, you should multiply
k by 2 \pi / Boxsize and P(k) by (Boxsize/2\pi)^3,
where Boxsize is the size of the box in Mpc/h.
In addition, CAMB uses a different Fourier convention to FFTW (which GenPK is based on),
so for comparison to CAMB you need to also multiply P(k) by (2\pi) ^3.
The total conversion factors are then
k *= 2\pi / Boxsize
P(k) *= Boxsize ^3
A script to do this can be found in plot.py
If you use this software in a paper, you can cite it via the ASCL:
http://ascl.net/1706.006