


Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b
Results
Initial program 1.9
Taylor expanded around 0 1.9
Final simplification1.9
herbie shell --seed 2018362 +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: 26.4s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 0.9b
Found 4 expressions with local error:
| 0.8b | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 0.7b | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 0.3b | (* y (log z)) |
| 0.1b | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 26.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 12.0ms | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 4.0ms | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 204.0ms | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 159.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 55.0ms | (* (- t 1.0) (log a)) |
| 38.0ms | (* y (log z)) |
34 calls. Slowest were:
| 409.0ms | (- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1)))) |
| 344.0ms | (+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y)))) |
| 246.0ms | (* (* (* (- t 1.0) (- t 1.0)) (- t 1.0)) (* (* (log a) (log a)) (log a))) |
5 alts after pruning (5 fresh and 0 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 0.8b | (exp (- (+ (* (log z) y) (* (- t 1.0) (log a))) b)) |
| 0.7b | (/ (* x (exp (- (+ (* (log z) y) (* (- t 1.0) (log a))) b))) y) |
| 0.3b | (* (log z) y) |
| 0.1b | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 64.0ms | (exp (- (+ (* (log z) y) (* (- t 1.0) (log a))) b)) |
| 29.0ms | (/ (* x (exp (- (+ (* (log z) y) (* (- t 1.0) (log a))) b))) y) |
| 7.0ms | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 207.0ms | (/ (* x (exp (- (+ (* (log z) y) (* (- t 1.0) (log a))) b))) y) |
| 156.0ms | (exp (- (+ (* (log z) y) (* (- t 1.0) (log a))) b)) |
| 62.0ms | (* (log z) y) |
| 57.0ms | (* (- t 1.0) (log a)) |
33 calls. Slowest were:
| 293.0ms | (+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y)))) |
| 284.0ms | (- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1)))) |
| 235.0ms | (* (* (* (- t 1.0) (- t 1.0)) (- t 1.0)) (* (* (log a) (log a)) (log a))) |
5 alts after pruning (4 fresh and 1 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 0.7b | (/ (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) y) |
| 0.3b | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
| 0.2b | (fma t (log a) (- b)) |
| 0.1b | (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) |
4 calls. Slowest were:
| 34.0ms | (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) |
| 25.0ms | (/ (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) y) |
| 24.0ms | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
4 calls. Slowest were:
| 408.0ms | (/ (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) y) |
| 362.0ms | (* x (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) |
| 324.0ms | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
| 73.0ms | (fma t (log a) (- b)) |
32 calls. Slowest were:
| 473.0ms | (* (* (* (* (pow a (- 1.0)) (pow a (- 1.0))) (pow a (- 1.0))) (* (* (pow z y) (pow z y)) (pow z y))) (* (* (exp (fma t (log a) (- b))) (exp (fma t (log a) (- b)))) (exp (fma t (log a) (- b))))) |
| 454.0ms | (+ (+ (* (log a) (- 1.0)) (* (log z) y)) (fma t (log a) (- b))) |
| 389.0ms | (+ (+ (log (pow a (- 1.0))) (log (pow z y))) (fma t (log a) (- b))) |
4 alts after pruning (3 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.7b | (/ (/ (* x (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0)) y) |
| 0.2b | (* (* 1 (pow z y)) (exp (fma t (log a) (- b)))) |
| 0.2b | (fma t (log a) (- b)) |
| 0.0b | (exp (fma t (log a) (- b))) |
4 calls. Slowest were:
| 12.0ms | (/ (/ (* x (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0)) y) |
| 9.0ms | (* (* 1 (pow z y)) (exp (fma t (log a) (- b)))) |
| 1.0ms | (exp (fma t (log a) (- b))) |
4 calls. Slowest were:
| 525.0ms | (/ (/ (* x (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0)) y) |
| 125.0ms | (* (* 1 (pow z y)) (exp (fma t (log a) (- b)))) |
| 58.0ms | (exp (fma t (log a) (- b))) |
| 58.0ms | (fma t (log a) (- b)) |
31 calls. Slowest were:
| 351.0ms | (* (exp (* -1 (* (log (/ 1 z)) y))) (exp (- (+ b (* t (log (/ 1 a))))))) |
| 323.0ms | (exp (- (* t (log -1)) (+ (* t (log (/ -1 a))) b))) |
| 313.0ms | (/ (* (exp (* -1 (* (log (/ 1 z)) y))) (* x (exp (- (+ b (* t (log (/ 1 a)))))))) (* a y)) |
4 alts after pruning (2 fresh and 2 done)
Merged error: 0.0b
0% (1.9b remaining)
Error of 1.9b against oracle of 0.0b and baseline of 1.9b