return x - scale * np.exp(log1p(-np.exp(-log_cdf / par_d)) / par_c)