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.

GPU Implementation of Conjugate Gradient Method

Parallel Conjugate Gradient Method Based on Spline Difference Method for the One-Dimensional Heat Equation

Matrix Decomposition Graphics Processing Unit Solver for Poisson Image Editing

A Practical GPU-based Implementation of the Katsevich Algorithm

Parallel L-BFGS-B Algorithm on GPU.

CUDA-based PCG algorithm optimization for a large sparse matrix

3-D Surface Quality Evaluation Based on Graphics Processing Unit

GPU-based Volume Rendering

GPU-Acceleration of Parallel Unconditionally Stable Group Explicit Finite Difference Method

Multilevel Interior Penalty Methods on GPUs

Parallel optimization for sparse matrix-vector on GPU

Research and Implementation of Effective Jacobi Iteration Algorithms on GPU

High Performance Computing Via a GPU

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

Generating Approximate Inverse Preconditioners for Sparse Matrices Using CUDA and GPGPU

Acceleration of Approximate Matrix Multiplications on GPUs

Improving Dense Linear Equation Solver on Hybrid CPU-GPU System.

Enhancing data locality of the conjugate gradient method for high-order matrix-free finite-element implementations

A New Sparse Matrix Vector Multiplication GPU Algorithm Designed for Finite Element Problems

Optimizing sparse matrix-vector multiplication based on gpu

AmgX: A Library for GPU Accelerated Algebraic Multigrid and Preconditioned Iterative Methods