Hybrid-norm and Fortran 2003: Separating the physics from the solver

Next: Implementation in Fortran 2003 Up: Hybrid-norm and Fortran 2003: Previous: Solvers

# Supporting the hybrid norm

The L2-norm is often a non-optimal choice both because of the non-Gaussian nature of noise in data and it's tendency to create smooth models with traditional regularization schemes. To improve model estimation, authors have suggested using either direct L1-norm approaches, or more commonly, Iteratively Reweighed Least-Squares (IRLS) to approximate different norms. IRLS attempts to find the best model at a non-L2 norm by a series of linearization s. Each non-linear iteration updates the weighting of the residual. For example we can achieve L1 like behavior by updating the weighting operator through,

 (10)

IRLS methods tend to be difficult to use because the user must choose carefully the number of relinearizations and the numbers of steps between relinearizations carefully. (Claerbout, 2009) suggested an alternate approach that dynamically changes the weighting function every iteration and uses a Taylor expansion of the standard conjugate direction algorithm to update the solution. Further he suggests a different norm, the hybrid norm, that creates a smooth transition between the standard L2 problem and a L1 problem. Given an error function and a residual vector the hybrid norm is defined

 (11)

where is a user supplied bad-data percentile.

Creating an inversion framework that supports a Taylor expansion approach to conjugate directions requires adding two additional features to our vector class description. First, we need to associate a norm to each vector. Second we need to be able to multiply a vector by another vector, element by element. Adding support for the hybrid norm requires more changes. A vector must now have a bad-data percentage associated with it, it must be able to find its ith percentile value, and create a vector with this value.

 Hybrid-norm and Fortran 2003: Separating the physics from the solver

Next: Implementation in Fortran 2003 Up: Hybrid-norm and Fortran 2003: Previous: Solvers

2010-11-26