irls - weighting functions
ierr=l1(res,weight)
ierr=cauchy(res,weight)
real(:)
Residual
real(:)
Resulting weighting function
Weighting function to apply to cgmethods (warning stability is not guaranteed when using these).
L1 - weight=1/abs(1+ res/rbar) CAUCHY- weight=1/sqrt(1+res/rbar**2)
geef90