Skip to content
Simple matlab code for testing optimization algorithms on seismic inverse problems
MATLAB Objective-C
Branch: master
Clone or download
Type Name Latest commit message Commit time
Failed to load latest commit information.
core update Jan 31, 2019
example Merge branch 'master' of Feb 11, 2019
operators added custom spot operators Jan 30, 2019
optimization update Jan 30, 2019
LICENSE Initial commit Oct 10, 2014 Update Oct 10, 2014
startup.m update Jan 31, 2019


This code provides the basic building blocks to test optimization algorithms on seismic inverse problems. The canonical seismic waveform inversion problem is given by

$$\min_{m} \sum_{i} \frac{1}{2}||P^TA(m)^{-1}q_i - d_i||_2^2 + \frac{\alpha}{2}||Lm||_2^2,$$

where $A(m)$ is the discretized Helmholtz operator $\omega^2 m + \nabla^2$ with absorbing boundary conditions, $P$ is the sampling operator, $q_i$ are the source vectors, $d_i$ are the observed data and $L$ is the discretized $\nabla$.

The main function is misfit.m, which takes as input the medium parameters $m$ and the data $d_i$ (and definitions of $P$, $q_i$, $\omega$ etc.) and returns the misfit value, its gradient and a function handle to compute the action of the Gauss-Newton hessian.

For an example, see /examples/marm.m, which uses a simple BB iteration to solve the optimization problem. Replace this with your favourite optimization algorithm and you're ready to roll!

Feel free to contact me with any questions, suggestions, etc.

Tristan van Leeuwen -

You can’t perform that action at this time.