use relative rather than absolute jitter at each step, for fitting p=1e-6 and smaller
use logging; fix --init=cov