


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.4424902966070684e-39 or -1.802008489021528e-303 < y Initial program 1.4
rmApplied associate-/l*1.4
if -1.4424902966070684e-39 < y < -1.802008489021528e-303Initial program 4.1
rmApplied associate-/l*3.6
Taylor expanded around inf 3.7
Simplified2.2
rmApplied pow-neg2.2
Applied associate-*l/2.2
Applied associate-*l/2.2
Applied associate-/r/2.2
rmApplied associate-/r*3.3
Final simplification1.8
herbie shell --seed 2018365 +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: 30.2s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 1.0b
Found 4 expressions with local error:
| 1.1b | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 0.8b | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 0.2b | (* y (log z)) |
| 0.2b | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 21.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 16.0ms | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 4.0ms | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 161.0ms | (/ (* x (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b))) y) |
| 98.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 56.0ms | (* (- t 1.0) (log a)) |
| 51.0ms | (* y (log z)) |
34 calls. Slowest were:
| 324.0ms | (- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1)))) |
| 291.0ms | (* (* (* (- t 1.0) (- t 1.0)) (- t 1.0)) (* (* (log a) (log a)) (log a))) |
| 274.0ms | (+ (* t (* (log a) (exp (- (* 1.0 (log a)))))) (+ (exp (- (* 1.0 (log a)))) (* (log z) (* (exp (- (* 1.0 (log a)))) y)))) |
3 alts after pruning (3 fresh and 0 done)
Merged error: 0.1b
Found 4 expressions with local error:
| 1.1b | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 1.0b | (/ x (/ y (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
| 0.2b | (* y (log z)) |
| 0.2b | (* (- t 1.0) (log a)) |
4 calls. Slowest were:
| 41.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 4.0ms | (* (- t 1.0) (log a)) |
| 2.0ms | (/ x (/ y (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
4 calls. Slowest were:
| 185.0ms | (/ x (/ y (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)))) |
| 114.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 60.0ms | (* (- t 1.0) (log a)) |
| 34.0ms | (* y (log z)) |
39 calls. Slowest were:
| 371.0ms | (- (+ (* t (log -1)) (* 1.0 (log (/ -1 a)))) (+ (* t (log (/ -1 a))) (* 1.0 (log -1)))) |
| 353.0ms | (exp (- (+ (* y (log z)) (* (- t 1.0) (log a))) b)) |
| 325.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 (3 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 1.0b | (/ x (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))))) |
| 0.3b | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
| 0.2b | (fma t (log a) (- b)) |
| 0.0b | (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) |
4 calls. Slowest were:
| 24.0ms | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
| 9.0ms | (/ x (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))))) |
| 9.0ms | (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) |
4 calls. Slowest were:
| 468.0ms | (/ x (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))))) |
| 261.0ms | (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b)))) |
| 243.0ms | (/ y (* (* (pow a (- 1.0)) (pow z y)) (exp (fma t (log a) (- b))))) |
| 50.0ms | (fma t (log a) (- b)) |
37 calls. Slowest were:
| 429.0ms | (+ (+ (* (log a) (- 1.0)) (* (log z) y)) (fma t (log a) (- b))) |
| 422.0ms | (+ (+ (* (log a) (- 1.0)) (log (pow z y))) (fma t (log a) (- b))) |
| 412.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))))) |
6 alts after pruning (5 fresh and 1 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 1.0b | (/ x (* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0))) |
| 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:
| 13.0ms | (* (* 1 (pow z y)) (exp (fma t (log a) (- b)))) |
| 3.0ms | (exp (fma t (log a) (- b))) |
| 2.0ms | (/ x (* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0))) |
4 calls. Slowest were:
| 542.0ms | (/ x (* (/ y (* (* 1 (pow z y)) (exp (fma t (log a) (- b))))) (pow a 1.0))) |
| 174.0ms | (* (* 1 (pow z y)) (exp (fma t (log a) (- b)))) |
| 86.0ms | (fma t (log a) (- b)) |
| 63.0ms | (exp (fma t (log a) (- b))) |
26 calls. Slowest were:
| 351.0ms | (* (exp (* -1 (* (log (/ 1 z)) y))) (exp (- (+ b (* t (log (/ 1 a))))))) |
| 335.0ms | (/ (* (exp (* -1 (* (log (/ 1 z)) y))) (* x (exp (- (+ b (* t (log (/ 1 a)))))))) (* a y)) |
| 328.0ms | (exp (- (* t (log -1)) (+ (* t (log (/ -1 a))) b))) |
6 alts after pruning (4 fresh and 2 done)
Merged error: 0.0b
3.1% (1.8b remaining)
Error of 1.8b against oracle of 0.0b and baseline of 1.9b