A Graphics Processing Unit Implementation and Optimization for Parallel Double-Difference Seismic Tomography

Abstract

Research Article| March 25, 2014 A Graphics Processing Unit Implementation and Optimization for Parallel Double‐Difference Seismic Tomography Pei‐Cheng Liao; Pei‐Cheng Liao aDepartment of Computer Science and Information Engineering, National Taiwan University of Science and Technology, No. 43, Sec. 4 Keelung Road, Taipei City, 106 Taiwan, ROCtom86046@gmail.comnucleargod823543@gmail.comyu-chi@mail.ntust.edu.tw Search for other works by this author on: GSW Google Scholar Cheng‐Chi Lii; Cheng‐Chi Lii aDepartment of Computer Science and Information Engineering, National Taiwan University of Science and Technology, No. 43, Sec. 4 Keelung Road, Taipei City, 106 Taiwan, ROCtom86046@gmail.comnucleargod823543@gmail.comyu-chi@mail.ntust.edu.tw Search for other works by this author on: GSW Google Scholar Yu‐Chi Lai; Yu‐Chi Lai aDepartment of Computer Science and Information Engineering, National Taiwan University of Science and Technology, No. 43, Sec. 4 Keelung Road, Taipei City, 106 Taiwan, ROCtom86046@gmail.comnucleargod823543@gmail.comyu-chi@mail.ntust.edu.tw Search for other works by this author on: GSW Google Scholar Ping‐Yu Chang; Ping‐Yu Chang bInstitute of Applied Geoscience, National Taiwan Ocean University, No. 2, Beining Road, Keelung, 20224 Taiwan, ROCpingyuc@mail.ntou.edu.tw Search for other works by this author on: GSW Google Scholar Haijiang Zhang; Haijiang Zhang cLaboratory of Seismology and Physics of Earth’s Interior, School of Earth and Space Sciences, University of Science and Technology of China, 96 Jinzhai Road, Hefei, Anhui 230026, P.R. Chinazhang11@ustc.edu.cn Search for other works by this author on: GSW Google Scholar Clifford Thurber Clifford Thurber dDepartment of Geoscience, University of Wisconsin–Madison, 1215 W. Dayton Street, Madison, Wisconsin 53706thurber@geology.wisc.edu Search for other works by this author on: GSW Google Scholar Author and Article Information Pei‐Cheng Liao aDepartment of Computer Science and Information Engineering, National Taiwan University of Science and Technology, No. 43, Sec. 4 Keelung Road, Taipei City, 106 Taiwan, ROCtom86046@gmail.comnucleargod823543@gmail.comyu-chi@mail.ntust.edu.tw Cheng‐Chi Lii aDepartment of Computer Science and Information Engineering, National Taiwan University of Science and Technology, No. 43, Sec. 4 Keelung Road, Taipei City, 106 Taiwan, ROCtom86046@gmail.comnucleargod823543@gmail.comyu-chi@mail.ntust.edu.tw Yu‐Chi Lai aDepartment of Computer Science and Information Engineering, National Taiwan University of Science and Technology, No. 43, Sec. 4 Keelung Road, Taipei City, 106 Taiwan, ROCtom86046@gmail.comnucleargod823543@gmail.comyu-chi@mail.ntust.edu.tw Ping‐Yu Chang bInstitute of Applied Geoscience, National Taiwan Ocean University, No. 2, Beining Road, Keelung, 20224 Taiwan, ROCpingyuc@mail.ntou.edu.tw Haijiang Zhang cLaboratory of Seismology and Physics of Earth’s Interior, School of Earth and Space Sciences, University of Science and Technology of China, 96 Jinzhai Road, Hefei, Anhui 230026, P.R. Chinazhang11@ustc.edu.cn Clifford Thurber dDepartment of Geoscience, University of Wisconsin–Madison, 1215 W. Dayton Street, Madison, Wisconsin 53706thurber@geology.wisc.edu Publisher: Seismological Society of America First Online: 14 Jul 2017 Online ISSN: 1943-3573 Print ISSN: 0037-1106 Bulletin of the Seismological Society of America (2014) 104 (2): 953–961. https://doi.org/10.1785/0120130236 Article history First Online: 14 Jul 2017 Cite View This Citation Add to Citation Manager Share Icon Share Facebook Twitter LinkedIn MailTo Tools Icon Tools Get Permissions Search Site Citation Pei‐Cheng Liao, Cheng‐Chi Lii, Yu‐Chi Lai, Ping‐Yu Chang, Haijiang Zhang, Clifford Thurber; A Graphics Processing Unit Implementation and Optimization for Parallel Double‐Difference Seismic Tomography. Bulletin of the Seismological Society of America 2014;; 104 (2): 953–961. doi: https://doi.org/10.1785/0120130236 Download citation file: Ris (Zotero) Refmanager EasyBib Bookends Mendeley Papers EndNote RefWorks BibTex toolbar search Search Dropdown Menu toolbar search search input Search input auto suggest filter your search All ContentBy SocietyBulletin of the Seismological Society of America Search Advanced Search Abstract Double‐difference seismic tomography can estimate velocity structure and event locations with high precision, but its high‐computation cost along with large memory usage prevents the use of a personal computer to process very large datasets and requires a long‐computation time. This work proposes graphics‐processing‐unit‐ (GPU)‐based acceleration schemes to run the algorithm on a personal computer for very large datasets more efficiently. Generally, the algorithm can be divided into five major steps: input, ray tracing, matrix construction, inversion, and output. This work focuses on accelerating the ray‐tracing and inversion steps, which take almost two‐thirds of the computation time. Before ray tracing, our algorithm preprocesses the data by sorting all recorded event–station paths according to their lengths. Therefore, those path estimation jobs assigned to GPU cores are suitable for the GPU architecture. Furthermore, our work also minimizes the usage of global and local memory to reduce the GPU computing time needed to handle a very large dataset. In addition to parallelizing the inversion computation, our work proposes a GPU‐based elimination method to reduce redundant computation in inversion for further acceleration. In our test, the proposed acceleration schemes can gain maximum speed‐up factors of 31.17 and 35.46 for ray tracing and inversion, respectively, in our test. Overall, the GPU‐based implementation can reach a maximum of 5.98 times faster than the central processing unit‐based implementation. You do not have access to this content, please speak to your institutional administrator if you feel you should have access.