Iterative Krylov solution methods for geophysical electromagnetic simulations on throughput-oriented processing units

Authors:
Michael Commer;Filipe Rnc Maia;Gregory A Newman
Affiliations:
Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, USA;National Energy Research Scientific Computing Center, Lawrence Berkeley National Laboratory, Berkeley, USA;Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, USA
Venue:
International Journal of High Performance Computing Applications
Year:
2012

Citing 5
Cited 1

Conjugate gradient-type methods for linear systems with complex symmetric coefficient matrices

SIAM Journal on Scientific and Statistical Computing - Special issue on iterative methods in numerical linear algebra
Improving performance of sparse matrix-vector multiplication

SC '99 Proceedings of the 1999 ACM/IEEE conference on Supercomputing
Massively parallel electrical-conductivity imaging of hydrocarbons using the IBM Blue Gene/L supercomputer

IBM Journal of Research and Development
Implementing sparse matrix-vector multiplication on throughput-oriented processors

Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis
Sparse Matrix-Vector Multiplication Optimizations based on Matrix Bandwidth Reduction using NVIDIA CUDA

DCABES '10 Proceedings of the 2010 Ninth International Symposium on Distributed Computing and Applications to Business, Engineering and Science

GPU parallelization of a three dimensional marine CSEM code

Computers & Geosciences

Quantified Score

Hi-index	0.00

Visualization

Abstract

Many geo-scientific applications involve boundary value problems arising in simulating electrostatic and electromagnetic fields for geophysical prospecting and subsurface imaging of electrical resistivity. Modeling complex geological media with three-dimensional finite-difference grids gives rise to large sparse linear systems of equations. For such systems, we have implemented three common iterative Krylov solution methods on graphics processing units and compared their performance with parallel host-based versions. The benchmarks show that the device efficiency improves with increasing grid sizes. Limitations are currently given by the device memory resources.