code-archive

AniCosmo: Bayesian analysis of anisotropic cosmologies

AniCosmo provides functionality to perform a Bayesian analysis of anisotropic cosmologies. It has been used to study Bianchi models of universal rotation in a number of papers. AniCosmo is not yet publicly available but will be made available in future. Please get in touch in you’re interested in the code.

BIANCHI: Bianchi ${VII}_h$ simulations

The BIANCHI code provides functionality to support the simulation of Bianchi Type $\text{VII}_h$ induced temperature fluctuations in CMB maps of a universe with shear and rotation. The implementation is based on the solutions to the Bianchi models derived by Barrow et al. (1985), which do not incorporate any dark energy component. Functionality is provided to compute the induced fluctuations on the sphere directly in either real or harmonic space. For a description of the equations implemented to perform the Bianchi simulations see the Appendix of our paper: Non-Gaussianity detections in the Bianchi $\text{VII}_h$ corrected WMAP 1-year data made with directional spherical wavelets.

BIANCHI2: Dark Bianchi ${VII}_h$ simulations

The BIANCHI2 code provides functionality to support the simulation of Bianchi Type $\text{VII}_h$ induced temperature fluctuations in CMB maps of a universe with shear and rotation. The implementation is based on the solutions to the Bianchi models derived by Anthony Lasenby (not published) that incorporate a cosmological constant. Functionality is provided to compute the induced fluctuations on the sphere directly in either real or harmonic space. Functionality to support functions defined on the sphere is provided by the S2 code.

COMB: Compact embedded object simulations

COMB provides functionality to support the simulation on the sphere of compact objects embedded in a stochastic background process of specified power spectrum. Support is provided to add additional white noise and convolve with beam functions. Functionality to support functions defined on the sphere is provided by the S2 code. COMB is typically used with the S2FIL code to generate template functions for the construction of optimal filters, such as the directional matched filter (MF) and scale adaptive filter (SAF) on the sphere, and to general simulated data with which to test object decection algorithms.

FastCSWT: Fast directional continuous spherical wavelet transform

FastCSWT provides functionality to perform a fast directional continuous wavelet transform on the sphere. The transform is based on the construction of the continuous spherical wavelet transform (CSWT) developed by Antoine and Vandergheynst (1999). Since the transform is based on a continuous framework, while it is possible to compute a wavelet transform of a signal of interest, it is not possible to reconstruct that signal from its wavelet coefficients. It order to consider exact synthesis, please see our S2LET code, which is based on an alternative scale-discretised wavelet framework on the sphere.

FLAG: Exact Fourier-Laguerre transform on the ball

FLAG provides a fast implementation of the Fourier-Laguerre transform, a novel 3D transform exploiting an exact quadrature rule on the ball to construct an exact harmonic transform in 3D spherical coordinates. The angular part of the Fourier-Laguerre transform uses the MW sampling theorem and the exact spherical harmonic transform implemented in the SSHT code. The radial sampling scheme arises from an exact quadrature on the radial half-line using damped Laguerre polynomials.

FLAGLET: Exact wavelets on the ball

FLAGLET provides efficient routines for the wavelet analysis of signals on the ball, the space formed by augmenting the sphere with the radial half-line. FLAGLET exploits the S2LET, FLAG and SSHT codes to assist with wavelet construction and support harmonic transforms on the ball and sphere, respectively. The flaglet transform is theoretically exact, i.e. the original signal can be synthesised from its wavelet coefficients exactly since the wavelet coefficients capture all the information of band-limited signals.

GLaSS: Generalised Lensing and Shear Spectra

The Generalised Lensing and Shear Spectra (GLaSS) code provides functionality to compute spherical Bessel, tomographic and generalised lensing spectra for weak gravitational lensing (i.e. cosmic shear) studies. GLaSS provides a module for Cosmosis and is called by Cosmosis. Run-mode options are specified in a Cosmosis ini file. GLaSS can either be run like any other Cosmosis module or interactively in Python. See the jupyter notebooks to see how to run interactively.

massmappy: Mapping dark matter on the celestial sphere

massmappy provides functionality to recover convergence mass maps on the celestial sphere from weak lensing cosmic shear observations. At present the spherical Kaiser-Squires estimator is implemented. Additional estimators will be added in future. massmappy also supports numerous projection operator to project the sphere to planar regions for comparison with the spherical Kaiser-Squires estimator. massmappy relies on the SSHT (built on the MW sampling theorem) and HEALPIX codes to handle sampled data on the sphere.

NSHT: Novel optimal sampling spherical harmonic transforms

NSHT implements a novel sampling scheme to accurately sample a signal band-limited at $L$ in $L^2$ samples. Whereas sampling theorems on the sphere require at best $\sim 2 L^2$ samples (see A novel sampling theorem on the sphere; implemented in SSHT), NSHT sampling achieves the optimal sampling rate specified by the dimensionality of band-limited signals in harmonic space. However, a theoretically exact sampling theorem in not obtained. Nevertheless, NSHT provides functionality to compute fast spherical harmonic transforms for signals sampled according to the optimal sampling scheme developed.