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