return x - scale * (np.sqrt(c * c - 2 * logsf) - c)