log
The log-density computed directly, so a point far in the tail keeps its true value where pdf would underflow to zero: (shape - 1) ln x - x/scale - shape ln(scale) - ln Gamma(shape) for x > 0, and the ln(Double.MIN_VALUE) floor where pdf is zero.