Fix log-likelihood of `gev` for cases near zero

[lbelzile]

Interpolate using third-order series expansion in the vicinity of the Gumbel sub-case, with potentially log1p function. See Paul Northrop's code