docs: public static double logpdf(double val, double loc, double scale, double shape) {
      val = (val - loc) / scale;
      double e = FastMath.exp(-val);
      return -(val + (shape + 1.0) * FastMath.log1p(e)) + FastMath.log(shape);
    }