Imperial College London
Portrait Picture

ProfessorGerardGorman

Faculty of Engineering, Department of Earth Science & Engineering

Professor of Computational Science and Engineering
 
 
 
//

Contact

 

+44 (0)20 7594 9985g.gorman Website

 
 
//

Location

 

R4.92Royal School of MinesSouth Kensington Campus

//

Summary

 

Publications

Citation

BibTex format

@article{Witte:2019:10.1190/GEO2018-0174.1,
author = {Witte, PA and Louboutin, M and Kukreja, N and Luporini, F and Lange, M and Gorman, GJ and Herrmann, FJ},
doi = {10.1190/GEO2018-0174.1},
journal = {Geophysics},
pages = {F57--F71},
title = {A large-scale framework for symbolic implementations of seismic inversion algorithms in Julia},
url = {http://dx.doi.org/10.1190/GEO2018-0174.1},
volume = {84},
year = {2019}
}
Download

RIS format (EndNote, RefMan)

TY  - JOUR
AB  - 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.
AU  - Witte,PA
AU  - Louboutin,M
AU  - Kukreja,N
AU  - Luporini,F
AU  - Lange,M
AU  - Gorman,GJ
AU  - Herrmann,FJ
DO  - 10.1190/GEO2018-0174.1
EP  - 71
PY  - 2019///
SN  - 0016-8033
SP  - 57
TI  - A large-scale framework for symbolic implementations of seismic inversion algorithms in Julia
T2  - Geophysics
UR  - http://dx.doi.org/10.1190/GEO2018-0174.1
UR  - http://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcApp=PARTNER_APP&SrcAuth=LinksAMR&KeyUT=WOS:000468946700026&DestLinkType=FullRecord&DestApp=ALL_WOS&UsrCustomerID=1ba7043ffcc86c417c072aa74d649202
UR  - https://library.seg.org/doi/full/10.1190/geo2018-0174.1
UR  - http://hdl.handle.net/10044/1/110211
VL  - 84
ER  -
Download