Multi-GPU maximum entropy image synthesis for radio astronomy
arXiv:1703.02920 · doi:10.1016/j.ascom.2017.11.003
Abstract
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 a statistical basis, it has a better resolution and better image quality under certain conditions. This work presents a high performance GPU version of non-gridding MEM, which is tested using real 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 allows to reconstruct the HD142527 CO(6-5) short baseline data set in 2.1 minutes, instead of 2.5 days that takes a sequential version on CPU.
arXiv comment11 pages, 13 figures
References in corpus (15)
- LOFAR: The LOw-Frequency ARray
- First Results from High Angular Resolution ALMA Observations Toward the HL Tau Region
- A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry
- Dust and Gas in the disc of HL Tauri: Surface density, dust settling, and dust-to-gas ratio
- Accretion kinematics through the warped transition disk in HD142527 from resolved CO(6-5) observations
- Morphological analysis of the cm-wave continuum in the dark cloud LDN1622
- Dynamical Stability of Imaged Planetary Systems in Formation: Application to HL Tau
- Optimal Image Reconstruction in Radio Interferometry
- Cygrid: A fast Cython-powered convolution-based gridding module for Python
- Bayesian Inference for Radio Observations
- Montblanc: GPU accelerated Radio Interferometer Measurement Equations in support of Bayesian Inference for Radio Observations
- Multifrequency VLBA study of the blazar S5 0716+714 during the active state in 2004 II. Large-scale jet kinematics and the comparison of the different methods of VLBI data imaging as applied to kinematic studies of AGN
- The Complex North Transition Region of Centaurus A: Radio Structure
- Bayesian Image Reconstruction Based on Voronoi Diagrams
- Imaging VLBI polarimetry data from Active Galactic Nuclei using the Maximum Entropy Method
Cited by in corpus (19)
- Dust unveils the formation of a mini-Neptune planet in a protoplanetary ring
- An inner warp in the DoAr 44 T Tauri transition disk
- Kinematic detections of protoplanets: a Doppler-flip in the disk of HD100546
- Long Baseline Observations of the HD100546 Protoplanetary Disk with ALMA
- A dusty filament and turbulent CO spirals in HD135344B-SAO206462
- Sparse Modeling in Quantum Many-Body Problems
- Cm-wavelength observations of MWC758: resolved dust trapping in a vortex
- ALMA Super-resolution Imaging of T Tau: r = 12 au Gap in the Compact Dust Disk around T Tau N
- Super-resolution Imaging of the Protoplanetary Disk HD 142527 Using Sparse Modeling
- Spirals and clumps in V960 Mon: signs of planet formation via gravitational instability around an FU Ori star?
- Resolved spectral variations of the centimetre-wavelength continuum from the rho Oph W photo-dissociation-region
- High-resolution ALMA observations of V4046Sgr: a circumbinary disc with a thin ring
- Regularized Maximum Likelihood Image Synthesis and Validation for ALMA Continuum Observations of Protoplanetary Disks
- Resolved observations at 31 GHz of spinning dust emissivity variations in Oph
- A young gas giant and hidden substructures in a protoplanetary disk
- Multi-frequency observations of PDS 70c: Radio emission mechanisms in the circum-planetary environment
- Molecules with ALMA at Planet-forming Scales (MAPS) II: CLEAN Strategies for Synthesizing Images of Molecular Line Emission in Protoplanetary Disks
- The Doppler-flip in HD100546 as a disk eruption: the elephant in the room of kinematic protoplanet searches
- Variable structure in the PDS70 disc and uncertainties in radio-interferometric image restoration