itrharrison/PyCosmo — pull request #8 by polyphant1 (2013-07-17 10:16:01+00:00)
Title: Class updates + void radius function

-----------------------------------
Progress as of 17th July
-----------------------------------

- Classes updated so that the Cosmology and PowSpec classes integrate effectively. When instantiating a Cosmology object, during initialisation a PowSpec instance is also created by passing the (half completed) Cosmology instance to the PowSpec instance. This allows the PowSpec instance to reference the methods in the Cosmology class (e.g. the growth function method).

-----------------------------------

- The void_distribution script has been updated to calculate the void radii distribution. The distribution of void sizes does not appear to match that shown in the Jennings, Li & Hu (2013) paper, and from looking at the elements of the number density formula (equation 10) I have narrowed it down to the multiplicity function f_sigma. There are 3 methods included to calculate f_sigma from JLH paper:

the exact infinite series (equation 6)
the JLH approximation (equation 8)
the SVDW approximation (equation B2)

I have included methods for each in the void_distribution.py script.

The exact and JLH approximations appear to underestimate the number of very small and very large voids, but show the correct spread of number densities when changing delta_c (the collapse limit, see figure 1 in the JLH paper). The SVDW approximation appears to show the correct distribution of voids, but a greater spread in large voids when changing delta_c.

I can't see any problems with my code, but a cursory glance might identify some silly mistakes I have missed. The only other area I can identify as a source of error is my calculation of sigma in the powspec.py script.

-----------------------------------