README: Convolution Lagrangian Perturbation Theory (CLPT) and Gaussian Streaming Redshift-Space Distortions (GSRSD)
This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
- Alejandro Aviles
- Martin White
This package contains two separated programs: CLPT and GSRSD. Those two programs are coded separately, but GSRSD uses the output of CLPT directly. For the overall reference, please refer to:
http://adsabs.harvard.edu/abs/2014MNRAS.437..588W
Compared to the previous version, the
The two packages are now parallelized with OpenMP.
CLPT and GSRSD are both written in standard C++-11 (also works with C++-0x; gcc version > 4.4 will do this correctly)–and nothing else. If you are working under Linux and have g++ and GNUMake in your computer, you can just get into each directory and type “make” respectively.
Both packages contain OpenMP parallelization. If you are using
Clang/LLVM, which does not officialy support OpenMP today, you may
want to remove the ”#include <omp.h>
” and ”omp_set_num_threads(
n_threads );
” statements in both ”main.cpp
” files. Or, instead, you
may want to try this:
The makefiles are written by a human, so you should find it easy to modify them as needed (e.g. if you have gcc 4.4, you may want to use the flag “-std=c++0x” instead of “-std=c++11” therein).
I named the output binary files as “.exe” files–the only purpose of this is for git to recogize them easily and clearly.
This program gives real-space statistics, including correlation functions and pairwise velocity statistics, coming from the CLPT.
For CLPT, please refer to:
http://adsabs.harvard.edu/abs/2013MNRAS.429.1674C
The CLPT code uses a parameter file to specify the parameters. There is an example of parameter file (“par.txt”) in the CLPT folder. You can modify the parameter file (e.g. adjust the values or the sequence of the keys), but please leave the header “CORR_FUNC” and “Init” intact.
Some explanations:
Key | Value |
---|---|
r_max | Maximum value of real-space separation that is calculated |
r_min | Minimum value of real-space separation |
r_bin_num | Number of data points between r_min and r_max |
pow_spec_file | Path to the input file: Linear power spectrum |
xi_file | Path to the output file: real-space |
v12_file | Path to the output file: real-space pairwise infall velocity |
s12_file | Path to the output file: real-space pairwise velocity dispersion |
k_input | If specified, functions of k are read from file; if “none” or missing, they will |
be calculated and stored into the file specified in key “k_output”. | |
k_output | The path for the |
q_input | Similar to “k_input”. Note: Functions of q depend on functions of |
q_output | Similar to “k_output”. |
n_thread | Number of threads that you will be using. If negative or missing, the number |
of threads will be determined by the env variable $OMP_NUM_THREADS |
The file for the velocity dispersions will be attached a suffix: “_p” for the parallel component, and “_v” for the perpandicular component.
The CLPT program also needs a liear power spectrum as input, whose path should be specified in the parameter file.
Definition of the file:
Column | Value |
---|---|
1 | Wave number (in |
2 | Power [in |
The name of the parameter file should be used as the parameter of command-line input. For example, if you are using Linux and your parameter file is named “par.txt”, you may want to type the followings in the terminal:
$ ./clpt.exe par.txt
In our CLPT model, we use
In this section we will use this notation.
Column | Value |
---|---|
1 |
|
2 |
|
3 | |
4 | |
5 | |
6 | |
7 | |
8 |
Veolicties are given in units of
Column | Value |
---|---|
1 |
|
2 |
|
3 | |
4 | |
5 | |
6 | |
7 | |
8 |
Column | Value |
---|---|
1 |
|
2 |
|
3 | |
4 | |
5 | |
6 |
|
7 | |
8 | |
9 |
Other components vanish at the desired (one-loop) order.
This program gives redshift-space statistics, including Legendre expansion of correlation function and wedge statistics, coming from the Gaussian streaming model.
For Gaussian streaming model, please refer to Beth & White (2011):
http://adsabs.harvard.edu/abs/2011MNRAS.417.1913R
For the “wedge” statistcis, you may find it in, such as,
http://adsabs.harvard.edu/abs/2011MNRAS.414..596W
The GSRSD code uses a parameter file to specify the parameters. There is an example of parameter file (“par.txt”) in the GSRSD folder. You can modify the parameter file (e.g. adjust the values or the sequence of the keys), but please leave the header “CORR_FUNC” and “Init” intact.
Some explanations:
Key | Value |
---|---|
fb11b20 | |
fb10b21 | |
fb11b21 | |
fb12b20 | |
fb10b22 | |
sigma_p_100 | The value of |
s_max | Maximum separation in the redshift-space |
s_min | Minimum separation |
s_bin | Number of data points between s_min and s_max |
wedge_bin | Number of data points for the “wedge” statistics |
y_spanning | The integration range of y (see Beth and White 2011) |
dy | The integration step of y |
f_v | Structure growth rate, |
xi_file_name | Path to the input file: real-space |
v_file_name | Path to the input file: real-space |
s_file_name | Path to the input file: real-space velocity dispersion |
out_file_name | Path to the output file |
n_thread | Number of threads that you will be using. If negative or missing, the number |
of threads will be determined by the env variable $OMP_NUM_THREADS |
The wedge statistics will be given by “<out_file_name>_wedge”, 2-D plot by “<out_file_name>_2d”, and some appendices will be given in “<out_file_name>_appendix”.
“fb1*b2*” are separately specified for the convenience of cross-correlation.
Similar to the CLPT code… E.g.:
$ ./xi_stream.exe par.txt
All for biased tracers…
Column | Value |
---|---|
1 |
|
2 |
|
3 |
|
4 |
|
All for biased tracers…
Column | Value |
---|---|
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
7 |
|
All for biased tracers…
Column | Value |
---|---|
1 |
|
Following | “Wedges” |
All for biased tracers…
Column | Value |
---|---|
1 |
|
2 |
|
3 |