If the initial p0 values are -np.inf due to the prior, it will iterate
through the other p0 values to replace the duff values. The replacement
takes another randomly selected p0 value and shifts it round by a
fractional gaussian with std=1e-10