I usually return NAN if the uniform random variate U is outside of [0, 1] and thus does not meet any of the conditions above. However, `rk_double(state)` should always return a value within [0, 1], and I don't see any of the other distributions in `distributions.c` returning NAN. Should this final `else { return NAN;}` condition be removed for consistency with the other distributions?