d = np.sqrt(p + alpha * np.abs(p.min()))