log1mexp function

Compute f(a) = loglog(1 +/- expexp(-a)) Numerically Optimally