License CC BY-NC-SA

Closure-phases

Package to generate an ionisation field in 2D or 3D and then compute its triangle correlations (ref https://arxiv.org/abs/1903.11402). Requires FFTW and OPENMP libraries.

First generate a bubble field with python (for a periodic 2D field with 512x512 pixels, with bubbles of radius 20 pixels, to reach a given filling fraction, allowed to overlap):

from Random_bubbles import *
cube=RandomBubbles(DIM=50, fillingfraction=0.05, radius=2, NDIM = 2, nooverlap = False, periodic=True) 
RandomBubbles.write_ionisation_field(cube) #saves real field in a text file

The output text file will be: Field_10binary_bubbles_nooverlap=False_radius=2_xhII0.05_N50_2D.txt

Then modify the header file SC.h with the correct information: box side length, computation parallelised on 20 cores, range of correlation scales considered (in Mpc)... Compile make 2d for 2D TCF make 3d for 3D TCF
Run the executable: ./SC.o - it will take a while.

The output will be a textfile with the r considered (hard coded, in Mpc), the real part and the imaginary part of the TCF for each r, and the number of modes used to compute the TCF for each r.

Note that you might have to modify the Makefile to include paths to your libraries. If you are using a Mac, we recommend using homebrew and installing both fftw and openmp (llvm) through homebrew. The paths to Homebrew libraries are commented out in the Makefile and given as an example.