log_cdf = -log(-expm1(logsf))