2014/10/14 by Gevorg Poghosyan, S. Matta, Sanchit Matta +5
Computer Science · Earth and Planetary Sciences · Physics and Astronomy · #Compiler #Computational Physics and Python Applications #Computational science #Computer science #Detector #Geophysics and Gravity Measurements #LIGO #Massively parallel #Message Passing Interface #Message passing #Operating system #Parallel computing #Pulsars and Gravitational Waves Research #Scalability #astro-ph.HE #cs.DC #cs.PF #cs.SE #gr-qc
paper · pdf · doi:10.1016/j.cpc.2014.10.025
published as Computer Physics Communications volume 188 pages 168 - 176 (2015) · 11 pages, 9 figures. Submitted to Computer Physics Communications
arxiv created 2014/10/14 · openalex publication_date 2014/11/17 · arxiv updated 2015/02/03 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/05
The parallelization, design and scalability of the \sky code to search for periodic gravitational waves from rotating neutron stars is discussed. The code is based on an efficient implementation of the F-statistic using the Fast Fourier Transform algorithm. To perform an analysis of data from the advanced LIGO and Virgo gravitational wave detectors' network, which will start operating in 2015, hundreds of millions of CPU hours will be required - the code utilizing the potential of massively parallel supercomputers is therefore mandatory. We have parallelized the code using the Message Passing Interface standard, implemented a mechanism for combining the searches at different sky-positions and frequency bands into one extremely scalable program. The parallel I/O interface is used to escape bottlenecks, when writing the generated data into file system. This allowed to develop a highly scalable computation code, which would enable the data analysis at large scales on acceptable time scales. Benchmarking of the code on a Cray XE6 system was performed to show efficiency of our parallelization concept and to demonstrate scaling up to 50 thousand cores in parallel.