modify adjust_grid to not remove peaks
@@ -256,6 +256,7 @@ class Distribution(object):
dy = self.y-self.peak_val
prune = interp_err[::2] > rel_tol*(1+ (dy[1:-1:2]/yc)**4)
ind[1:-1:2] = prune
+ ind[self.peak_idx] = True
if np.mean(prune)<1.0:
self._func.y = self._func.y[ind]
self._func.x = self._func.x[ind]
