1 - \log \left(1 - \frac{x - y}{1 - y}\right)double code(double x, double y) {
return ((double) (1.0 - ((double) log(((double) (1.0 - ((double) (((double) (x - y)) / ((double) (1.0 - y))))))))));
}
Please include this information when filing a bug report:
herbie shell --seed 2020169
(FPCore (x y)
:name "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, B"
:precision binary64
:herbie-target
(if (< y -81284752.61947241) (- 1.0 (log (- (/ x (* y y)) (- (/ 1.0 y) (/ x y))))) (if (< y 3.0094271212461764e+25) (log (/ (exp 1.0) (- 1.0 (/ (- x y) (- 1.0 y))))) (- 1.0 (log (- (/ x (* y y)) (- (/ 1.0 y) (/ x y)))))))
(- 1.0 (log (- 1.0 (/ (- x y) (- 1.0 y))))))
| match-define: no matching clause for '(== branch1333 -228783269951518.34) | L | C | |
|---|---|---|---|
| match:error | /opt/racket-7.5/collects/racket/match/runtime.rkt | 24 | 0 |
| program-body | /data/pavpan/nightlies/herbie/ival-find-range/src/programs.rkt | 25 | 0 |
| (unnamed) | /opt/racket-7.5/collects/racket/contract/private/arrow-higher-order.rkt | 379 | 33 |
| prepare-points | /data/pavpan/nightlies/herbie/ival-find-range/src/points.rkt | 130 | 0 |
| pred | /data/pavpan/nightlies/herbie/ival-find-range/src/core/regimes.rkt | 204 | 4 |
| .../more-scheme.rkt:261:28 | unknown | ||
| for-loop | /data/pavpan/nightlies/herbie/ival-find-range/src/core/regimes.rkt | 240 | 3 |
| sindices->spoints | /data/pavpan/nightlies/herbie/ival-find-range/src/core/regimes.rkt | 194 | 0 |
| f37 | unknown | ||
| get-final-combination | /data/pavpan/nightlies/herbie/ival-find-range/src/mainloop.rkt | 389 | 0 |
| .../more-scheme.rkt:261:28 | unknown | ||
| run | /opt/racket-7.5/share/pkgs/profile-lib/main.rkt | 39 | 2 |
| profile-thunk16 | /opt/racket-7.5/share/pkgs/profile-lib/main.rkt | 9 | 0 |
| (unnamed) | /opt/racket-7.5/collects/racket/private/more-scheme.rkt | 261 | 28 |
herbie shell --seed 2020169
(FPCore (x y)
:name "Numeric.SpecFunctions:invIncompleteGamma from math-functions-0.1.5.2, B"
:precision binary64
:herbie-target
(if (< y -81284752.61947241) (- 1.0 (log (- (/ x (* y y)) (- (/ 1.0 y) (/ x y))))) (if (< y 3.0094271212461764e+25) (log (/ (exp 1.0) (- 1.0 (/ (- x y) (- 1.0 y))))) (- 1.0 (log (- (/ x (* y y)) (- (/ 1.0 y) (/ x y)))))))
(- 1.0 (log (- 1.0 (/ (- x y) (- 1.0 y))))))