New radial window functions linear_windows() and cubic_windows(), which correspond to linear and cubic spline interpolation of radial functions, respectively. These are overlapping window functions, and it has been difficult to obtain accurate matter power spectra so far.
The partition() function now returns an array of weights to approximate a given function by the window functions. This is necessary to obtain an accurate fit of redshift distributions by overlapping window functions. For example, to get the array of galaxy densities in each shells from dndz, one would now do:
ngal = partition(z, dndz, shells)
A new function getcl() was added to return angular power spectra by index from a list using GLASS ordering.
The gaussian_phz() function now accepts bounds using lower= and upper= keyword parameters.
glass 2023.7 (1 Aug 2023)
New radial window functions
linear_windows()
andcubic_windows()
, which correspond to linear and cubic spline interpolation of radial functions, respectively. These are overlapping window functions, and it has been difficult to obtain accurate matter power spectra so far.The
partition()
function now returns an array of weights to approximate a given function by the window functions. This is necessary to obtain an accurate fit of redshift distributions by overlapping window functions. For example, to get the array of galaxy densities in each shells fromdndz
, one would now do:ngal = partition(z, dndz, shells)
A new function
getcl()
was added to return angular power spectra by index from a list using GLASS ordering.The
gaussian_phz()
function now accepts bounds usinglower=
andupper=
keyword parameters.