Exploring the use of SPIKE-based solvers on large electromagnetic modeling
MetadataShow full item record
Frequency-domain seismic and electromagnetic modeling requires solving the linear systems resulting from the discretization of the corresponding time-harmonic equations. Geophysical inversion is typically performed using several discrete frequencies and multiple (up to tens of thousands) source/receiver combinations. Limitations of classical direct and iterative sparse linear solvers have caused the development of the so-called hybrid methods that can be viewed as an intermediate approach between the direct and iterative methods. We present an efficient parallel solver based on the SPIKE algorithm. Several examples in frequency domain electromagnetic modeling illustrate the computational efficiency of the developed method in terms of memory demand and floating-point operations. Multiple sources can be efficiently handled by employing sparse direct solvers in the factorization of diagonal blocks of the system matrix. Based on the divide and conquer idea, this kind of algorithms exposes different parallelism levels, being suitable to take advantage of multiple accelerator devices. The SPIKE solver partially overcomes the fill-in problem of direct solvers, allowing to solve much larger domains on the same system.
Showing items related by title, author, creator and subject.
Paszynska, A.; Jopek, K.; Banas, K.; Paszynski, M.; Gurgul, P.; Lenerth, A.; Nguyen, D.; Pingali, K.; Dalcin, L.; Calo, Victor (2015)This paper describes a telescopic solver for two dimensional h adaptive grids with point singularities. The input for the telescopic solver is an h refined two dimensional computational mesh with rectangular finite elements. ...
Pardo, D.; Paszynski, M.; Collier, N.; Alvarez, J.; Dalcin, L.; Calo, Victor (2012)In this paper we describe the history, performance, and design concepts of direct solvers for algebraic systems resulting from Galerkin discretizations of partial differential equations. Popular direct solver implementations ...
Garcia, D.; Pardo, D.; Dalcin, L.; Calo, Victor (2018)Starting from a highly continuous Isogeometric Analysis (IGA) discretization, refined Isogeometric Analysis (rIGA) introduces C 0 hyperplanes that act as separators for the direct LU factorization solver. As a result, the ...