1. Seager et al., 1999,  ApJ, 523, L1–L5

  2. Chluba & Thomas, 2010, arXiv:1010.3631

  3. Rubino-Martin et al., 2010, MNRAS, 403, 439-452

  4. Chluba, 2010, MNRAS, 402, 1195–1207

  5. Chluba, Vasil & Dursi, 2010, MNRAS, 407, 599–612

as well as:

  1. Switzer & Hirata, 2008, Phys.Rev.D. 77, 083006

  2. Grin & Hirata, 2010, Phys.Rev.D, 81, 083005

  3. Ali-Haimoud & Hirata, 2010, Phys.Rev.D, 82, 063521

This page introduces my C++ version of Recfast. This code is based on the C version by Seager et al. 1999 (the original version can be found here), but was strongly rewritten and updated to:

  1. perform the normal Recfast computation (Recfast v1.4.2 with normal helium rec)

  2. include the recombination corrections given by Chluba & Thomas 2010

  3. allow running a simple Dark matter annihilation module

These tasks are achieved using the simple 3-level approach of Recfast. To solve the system of ODE's we use a simplified version of the stiff ODE-solver developed by Chluba, Vasil & Dursi, 2010. This solver allows us to avoid the switches in the original Recfast code.

To download the code, please click here. This tar-ball contains a makefile and detailed explanation about how to compile and run the code. Please feel free to use, modify and distribute it, however, when using the code please refer to:

Fig. 1: Dependence of the cosmological ionization history on the efficiency of dark matter annihilation. The curves were obtained running Recfast++ with the dark matter module. Recombination corrections were not included.

weebly analytics