solver_reg - iteratively solve regularized system of equation
call solver_reg(oper, solv, reg, nreg, x, dat, niter, eps,
x0,nmem,nfreq,err,res,nloper,xmov,rmov,wght,verb,known,wt)
operator (standard interface)
step mecanism
regularization operator operator(standard interface)
number of output pts when reg is applied to the model
real(:)
data
number of iterations
real(:)
Initial model
Change behavior after nmem itterations use in conjunction with nfreq and/or wght
Restart step function every nfreq itterations after nmem steps
Non-linear operator to use for forward operator
func(res,m)
Weighting funtion operator
Wheter or not to be verbose
logical(:)
Specify (.true.) model values that are knonw
real(:)
Weighting operator to apply to residuals
real(:)
model
real(niter)
Error sum((d-d~)**2)
real(:)
Final residual (size of data)
real(:,niter)
Movie of model values as a function of iter
real(:,niter)
Movie of data residual
Generic regularized solver.
All operators must follow convetion
integer function (logical adj, logical add, real(:) model, real(:)data)
All solv routines must follow the convention
integer function(logical forget,real(:) x,real(:) g, real(:) rr,real(:) gg)
the solver_prec manpage,the solver manpage,the nonlin_solver manpage,the lsqr manpage
geef90