Combining checkpointing and data compression for large scale seismic inversion
Authors:
Preprint
Abstract
Seismic inversion and imaging are adjoint-based optimization problems that processes up to terabytes of data, regularly exceeding the memory capacity of available computers. Data compression is an effective strategy to reduce this memory requirement by a certain factor, particularly if some loss in accuracy is acceptable. A popular alternative is checkpointing, where data is stored at selected points in time, and values at other times are recomputed as needed from the last stored state. This allows arbitrarily large adjoint computations with limited memory, at the cost of additional recomputations. In this paper we combine compression and checkpointing for the first time to compute a realistic seismic inversion. The combination of checkpointing and compression allows larger adjoint computations compared to using only compression, and reduces the recomputation overhead significantly compared to using only checkpointing.Keywords
All Related Versions
This publication has 13 references indexed in Scilit:
- Architecture and Performance of Devito, a System for Automated Stencil ComputationACM Transactions on Mathematical Software, 2020
- Efficient Lossy Compression for Scientific Data Based on Pointwise Relative Error BoundIEEE Transactions on Parallel and Distributed Systems, 2018
- Optimal Multistage Algorithm for Adjoint ComputationSIAM Journal on Scientific Computing, 2016
- Fixed-Rate Compressed Floating-Point ArraysIEEE Transactions on Visualization and Computer Graphics, 2014
- RTM using effective boundary saving: A staggered grid GPU implementationComputers & Geosciences, 2014
- Seismic wave modeling for seismic imagingThe Leading Edge, 2009
- Minimal Repetition Dynamic Checkpointing Algorithm for Unsteady Adjoint CalculationSIAM Journal on Scientific Computing, 2009
- MultiStage Approaches for Optimal Offline CheckpointingSIAM Journal on Scientific Computing, 2009
- A review of the adjoint-state method for computing the gradient of a functional with geophysical applicationsGeophysical Journal International, 2006
- The JPEG 2000 still image compression standardIEEE Signal Processing Magazine, 2001