Radio Galaxy Detection in the Visibility Domain
Bayesian model-fitting method working in the Fourier domain adopting a single SF galaxy model (the exponential profile, i.e. Sersic of index 1). The resulting multimodal posterior distribution is sampled using a multimodal nested sampling algorithm. See the related paper by A. Malyali, M. Rivi, F.B. Abdalla, J.D. McEwen, 2019, MNRAS 486(2) 2695–2704.
This is the C version of the python code. MPI + OpenMP parallelization. where each MPI task will read a different spectral window of the dataset that must be split in independent MS files. MPI can be enabled/disabled in the Makefile.
MultiNest output files must be processed to remove fake modes using the python code Modes_selection.py, specifying the flux/arcsec threshold.
Update the Makefile with the path to the dependencies, the run make.
- casacore library for reading Measurement Set file format.
- MultiNest library.
- MPI library for multi-node parallelization
pyMultiNest package for reading MultiNest output files.