


Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b
if y < -1.9813937686979932e-95 or 0.00920090525632143 < y Initial program 0.4
rmApplied associate-/l*0.4
rmApplied *-un-lft-identity0.4
Applied exp-prod0.4
Simplified0.4
if -1.9813937686979932e-95 < y < 0.00920090525632143Initial program 3.7
Taylor expanded around inf 3.7
Simplified2.4
rmApplied add-sqr-sqrt2.4
Applied associate-*r*2.4
Final simplification1.3
herbie shell --seed 2019008 +o rules:numerics
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"
(/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y))
Time bar (total: 43.9s)Debug log
1 calls. Slowest were:
| 384.0ms | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 0.8b
Found 4 expressions with local error:
| 0.9b | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 0.9b | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 0.3b | (* y (log z)) |
| 0.1b | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 23.0ms | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 23.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 8.0ms | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 212.0ms | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 142.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 77.0ms | (* (- t 1.0) (log a)) |
| 35.0ms | (* y (log z)) |
34 calls. Slowest were:
| 340.0ms | (- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1)))) |
| 308.0ms | (+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y)))) |
| 201.0ms | (* (* (* (- t 1.0) (- t 1.0)) (- t 1.0)) (* (* (log a) (log a)) (log a))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 1.2b | (/ x (/ y (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
| 0.9b | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 0.3b | (* y (log z)) |
| 0.1b | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 42.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 5.0ms | (/ x (/ y (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
| 4.0ms | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 162.0ms | (/ x (/ y (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
| 105.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 72.0ms | (* (- t 1.0) (log a)) |
| 51.0ms | (* y (log z)) |
39 calls. Slowest were:
| 353.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 340.0ms | (- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1)))) |
| 323.0ms | (+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y)))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 1.2b | (/ x (/ y (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
| 0.9b | (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 0.3b | (* y (log z)) |
| 0.1b | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 7.0ms | (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 4.0ms | (* (- t 1.0) (log a)) |
| 2.0ms | (/ x (/ y (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
4 calls. Slowest were:
| 204.0ms | (/ x (/ y (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
| 125.0ms | (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 46.0ms | (* (- t 1.0) (log a)) |
| 40.0ms | (* y (log z)) |
38 calls. Slowest were:
| 303.0ms | (- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1)))) |
| 241.0ms | (sqrt (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 237.0ms | (/ 1 (pow E (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) |
4 alts after pruning (3 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.9b | (/ (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) y) |
| 0.2b | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
| 0.1b | (fma t (log a) (- b)) |
| 0.0b | (* (pow a (- 1.0)) (pow z y)) |
4 calls. Slowest were:
| 23.0ms | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
| 21.0ms | (/ (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) y) |
| 8.0ms | (* (pow a (- 1.0)) (pow z y)) |
4 calls. Slowest were:
| 447.0ms | (/ (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) y) |
| 279.0ms | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
| 174.0ms | (* (pow a (- 1.0)) (pow z y)) |
| 63.0ms | (fma t (log a) (- b)) |
41 calls. Slowest were:
| 763.0ms | (/ (* (exp (* -1 (* (log (/ 1 z)) y))) (* x (exp (- (+ b (* t (log (/ 1 a)))))))) (* a y)) |
| 387.0ms | (+ (+ (* (log a) (- 1.0)) (* (log z) y)) (fma t (log a) (- b))) |
| 369.0ms | (+ (+ (* (log a) (- 1.0)) (log (pow z y))) (fma t (log a) (- b))) |
5 alts after pruning (4 fresh and 1 done)
Merged error: 0b
32.6% (1.3b remaining)
Error of 1.3b against oracle of 0.0b and baseline of 1.9b