docs: public static double getLowerBoundForBoverA(final long a, final long b, final double f) {
      checkInputs(a, b, f);
      if (a == 0) { return 0.0; }
      if (f == 1.0) { return (double) b / a; }
      return approximateLowerBoundOnP(a, b, NUM_STD_DEVS * hackyAdjuster(f));
    }