Exploring the use of SPIKE-based solvers on large electromagnetic modeling
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/rece...
| Main Authors: | , , , |
|---|---|
| Format: | Conference Paper |
| Published: |
2016
|
| Online Access: | http://hdl.handle.net/20.500.11937/54484 |
| _version_ | 1848759382698885120 |
|---|---|
| author | Rodriguez Bernabeu, S. Puzyrev, Vladimir Hanzich, M. Fernández, S. |
| author_facet | Rodriguez Bernabeu, S. Puzyrev, Vladimir Hanzich, M. Fernández, S. |
| author_sort | Rodriguez Bernabeu, S. |
| building | Curtin Institutional Repository |
| collection | Online Access |
| description | 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. |
| first_indexed | 2025-11-14T09:59:00Z |
| format | Conference Paper |
| id | curtin-20.500.11937-54484 |
| institution | Curtin University Malaysia |
| institution_category | Local University |
| last_indexed | 2025-11-14T09:59:00Z |
| publishDate | 2016 |
| recordtype | eprints |
| repository_type | Digital Repository |
| spelling | curtin-20.500.11937-544842017-07-27T05:21:16Z Exploring the use of SPIKE-based solvers on large electromagnetic modeling Rodriguez Bernabeu, S. Puzyrev, Vladimir Hanzich, M. Fernández, S. 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. 2016 Conference Paper http://hdl.handle.net/20.500.11937/54484 restricted |
| spellingShingle | Rodriguez Bernabeu, S. Puzyrev, Vladimir Hanzich, M. Fernández, S. Exploring the use of SPIKE-based solvers on large electromagnetic modeling |
| title | Exploring the use of SPIKE-based solvers on large electromagnetic modeling |
| title_full | Exploring the use of SPIKE-based solvers on large electromagnetic modeling |
| title_fullStr | Exploring the use of SPIKE-based solvers on large electromagnetic modeling |
| title_full_unstemmed | Exploring the use of SPIKE-based solvers on large electromagnetic modeling |
| title_short | Exploring the use of SPIKE-based solvers on large electromagnetic modeling |
| title_sort | exploring the use of spike-based solvers on large electromagnetic modeling |
| url | http://hdl.handle.net/20.500.11937/54484 |