occam: Occam inversion

In PEIP: Geophysical Inverse Theory and Optimization

Occam inversion

Description

Occam's inversion

Usage

occam(afun, ajac, L, d, m0, delta)

Arguments

afun	character, function handle that computes the forward problem
ajac	character, function handle that computes the Jacobian of the forward problem
L	regularization matrix
d	data that should be fit
m0	guess at the model
delta	cutoff to use for the discrepancy principle portion

Value

vector, model found

Note

This is a simple brute force way to do the line search. Much more sophisticated methods are available. Note: we've restricted the line search to the range from 1.0e-20 to 1. This seems to work well in practice, but might need to be adjusted for a particular problem.

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

bayes

PEIP documentation built on Aug. 21, 2023, 9:10 a.m.