docs: public static double logcdf(double x, double k, double theta, double shift) {
      x = (x - shift);
      return x <= 0. ? -Double.NEGATIVE_INFINITY : GammaDistribution.logregularizedGammaP(k, FastMath.log1p(x) * theta);
    }