Multi-GPU maximum entropy image synthesis for radio astronomy

The maximum entropy method (MEM) is a well known deconvolution technique in radio-interferometry. This method solves a non-linear optimization problem with an entropy regularization term. Other heuristics such as CLEAN are faster but highly user dependent. Nevertheless, MEM has the following advantages: it is unsupervised, it has an statistical basis, it has a better resolution and better image quality under certain conditions. This work presents a high performance GPU version of non-gridded MEM, which is tested using interferometric and simulated data. We propose a single-GPU and a multi-GPU implementation for single and multi-spectral data, respectively. We also make use of the Peer-to-Peer and Unified Virtual Addressing features of newer GPUs which allows to exploit transparently and efficiently multiple GPUs. Several ALMA data sets are used to demonstrate the effectiveness in imaging and to evaluate GPU performance. The results show that a speedup from 1000 to 5000 times faster than a sequential version can be achieved, depending on data and image size. This has allowed us to reconstruct the HD142527 CO(6-5) short baseline data set in 2.1 minutes, instead of the 2.5 days that takes on CPU.

Comments: 11 pages, 13 figures

Similar Publications

Pulsar Wind Nebulae (PWNe) constitute an ideal astrophysical environment to test our current understanding of relativistic plasma processes. It is well known that magnetic fields play a crucial role in their dynamics and emission properties. At present, one of the main issues concerns the level of magnetic turbulence present in these systems, which in the absence of space resolved X-ray polarization measures cannot be directly constrained. Read More

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures to map the faint Cherenkov light emitted in extensive air showers onto their image sensors. Segmented reflectors fulfill these needs using mass produced and light weight mirror facets. However, as the overall image is the sum of the individual mirror facet images, alignment is important. Read More

For many years, lunar laser ranging (LLR) observations using a green wavelength have suffered an inhomogeneity problem both temporally and spatially. This paper reports on the implementation of a new infrared detection at the Grasse LLR station and describes how infrared telemetry improves this situation. Our first results show that infrared detection permits us to densify the observations and allows measurements during the new and the full Moon periods. Read More

The machine learning techniques are widely applied in many modern optical sky surveys, i.e. Pan-STARRS1, PTF/iPTF and Subaru/Hyper Suprime-Cam survey, to reduce the human intervention for data verification. Read More

Different combinations of input parameters to filament identification algorithms, such as Disperse and FilFinder, produce numerous different output skeletons. The skeletons are a one pixel wide representation of the filamentary structure in the original input image. However, these output skeletons may not necessarily be a good representation of that structure. Read More

Determining the velocity distribution of halo stars is essential for estimating the mass of the Milky Way and for inferring its formation history. Since the stellar halo is a dynamically hot system, the velocity distribution of halo stars is well described by the 3-dimensional velocity dispersions $(\sigma_r, \sigma_\theta, \sigma_\phi)$, or by the velocity anisotropy parameter $\beta=1-(\sigma_\theta^2+\sigma_\phi^2)/(2\sigma_r^2)$. Direct measurements of $(\sigma_r, \sigma_\theta, \sigma_\phi)$ consistently suggest $\beta =0. Read More

POLAR is space-borne detector designed for a precise measurement of gamma-ray polarization of the prompt emissions of Gamma-Ray Bursts in the energy range 50 keV - 500 keV. POLAR is a compact Compton polarimeter consisting of 40$\times$ 40 plastic scintillator bars read out by 25 multi-anode PMTs. In May 2015, we performed a series of tests of the POLAR flight model with 100\% polarized x-rays beams at the European Synchrotron Radiation Facility beam-line ID11 aming to study thresholds, crosstalk between channels and responses of POLAR flight model to polarized X-ray beams. Read More

A future compact and modular X and gamma-ray spectrometer (XGS) has been designed and a series of prototypes have been developed and tested. The experiment envisages the use of CsI scintillator bars read out at both ends by single-cell 25 mm2 Silicon Drift Detectors. Digital algorithms are used to discriminate between events absorbed in the Silicon layer (lower energy X rays) and events absorbed in the scintillator crystal (higher energy X rays and gamma-rays). Read More

Extensive air showers still are our only access to the highest-energy particles in the universe, namely cosmic-ray nuclei with energies up to several 100 EeV. Studying open questions in cosmic-ray physics, like their yet unknown origin requires the reconstruction of the energy and mass of the primary particles from the air-shower measurements. Great progress has been achieved lately in the development of the radio detection technique for this purpose. Read More