tfb_gumbel_cdf function

Compute Y = g(X) = exp(-exp(-(X - loc) / scale)), the Gumbel CDF.