Up: R. Clapp: Regularization
Instead of approximating the Hessian
with a diagonal matrix we can attempt to
estimate the least squares inverse
using a conjugate gradient solver.
The model is preconditioned by using
polynomial division to apply the helical derivative
and the new preconditioned variable is estimated through