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