For explink with deriv = 0,
  the exponential of theta, i.e.,
exp(theta) when inverse = FALSE.
  And if inverse = TRUE then
log(theta);
  if theta is not positive then it will return NaN.
For deriv = 1, then the function returns
d
eta / d
theta as a
  function of theta
if inverse = FALSE,
  else if inverse = TRUE then it returns the reciprocal.
Here, all logarithms are natural logarithms,
  i.e., to base e.