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