docs: public static double logcdf(double val, double k, double theta) {
      if(val < 0) {
        return Double.NEGATIVE_INFINITY;
      }
      double vt = val * theta;
      return (val == Double.POSITIVE_INFINITY) ? 0. : logregularizedGammaP(k, vt);
    }