improvement vs baseline scatter
fit failures
use tol=1E-3 in minimization
minimization now gets best P even when the procedure later derails
changed fit_loo() to do a single iteration for each held out point, starting from parameters of the overall fit