log_cdf = log(-expm1(logsf))