A large-scale framework for symbolic implementations of seismic inversion algorithms in Julia
| Autor: | Gerard J. Gorman, Philipp Witte, Felix J. Herrmann, Mathias Louboutin, Navjot Kukreja, Fabio Luporini, Michael Lange |
|---|---|
| Rok vydání: | 2019 |
| Předmět: |
Scale (ratio)
Computer science business.industry Inversion (meteorology) 010103 numerical & computational mathematics 010502 geochemistry & geophysics 01 natural sciences Task (project management) Geophysics Software Geochemistry and Petrology Seismic inversion Time domain 0101 mathematics business Implementation Algorithm 0105 earth and related environmental sciences |
| Zdroj: | GEOPHYSICS. 84:F57-F71 |
| ISSN: | 1942-2156 0016-8033 |
| Popis: | Writing software packages for seismic inversion is a very challenging task because problems such as full-waveform inversion or least-squares imaging are algorithmically and computationally demanding due to the large number of unknown parameters and the fact that waves are propagated over many wavelengths. Therefore, software frameworks need to combine versatility and performance to provide geophysicists with the means and flexibility to implement complex algorithms that scale to exceedingly large 3D problems. Following these principles, we have developed the Julia Devito Inversion framework, an open-source software package in Julia for large-scale seismic modeling and inversion based on Devito, a domain-specific language compiler for automatic code generation. The framework consists of matrix-free linear operators for implementing seismic inversion algorithms that closely resemble the mathematical notation, a flexible resilient parallelization, and an interface to Devito for generating optimized stencil code to solve the underlying wave equations. In comparison with many manually optimized industry codes written in low-level languages, our software is built on the idea of independent layers of abstractions and user interfaces with symbolic operators. Through a series of numerical examples, we determined that this allows users to implement a series of increasingly complex algorithms for waveform inversion and imaging as simple Julia scripts that scale to large-scale 3D problems. This illustrates that software based on the paradigms of abstract user interfaces and automatic code generation and makes it possible to manage the complexity of the algorithms and performance optimizations, thus providing a high-performance research and production framework. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|