



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.2 |
|---|---|
| Target | 20.7 |
| Herbie | 6.6 |
if b < -1.1121654501400248e+139Initial program 61.2
Taylor expanded around -inf 1.6
Simplified1.6
if -1.1121654501400248e+139 < b < -6.179298479433663e-282Initial program 34.1
rmApplied div-inv34.2
rmApplied flip--34.3
Applied associate-*l/34.3
Simplified13.9
Taylor expanded around 0 7.7
if -6.179298479433663e-282 < b < 7.734953967976876e+108Initial program 9.5
rmApplied div-inv9.7
if 7.734953967976876e+108 < b Initial program 46.6
rmApplied div-inv46.6
rmApplied flip--61.4
Applied associate-*l/61.4
Simplified61.6
Taylor expanded around 0 61.5
Taylor expanded around 0 3.9
Simplified3.9
Final simplification6.6
herbie shell --seed 2019008
(FPCore (a b c)
:name "quadm (p42, negative)"
:herbie-target
(if (< b 0) (/ c (* a (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))) (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
(/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
Time bar (total: 22.8s)Debug log
1 calls. Slowest were:
| 448.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.0b
Found 4 expressions with local error:
| 23.7b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 1.0b | (- (* b b) (* 4 (* a c))) |
| 0.8b | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.1b | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
4 calls. Slowest were:
| 21.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 8.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 4.0ms | (- (* b b) (* 4 (* a c))) |
4 calls. Slowest were:
| 80.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 71.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 53.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 32.0ms | (- (* b b) (* 4 (* a c))) |
46 calls. Slowest were:
| 322.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 276.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 234.0ms | (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 23.7b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 1.0b | (- (* b b) (* 4 (* a c))) |
| 0.8b | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.2b | (* (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (/ 1 (* 2 a))) |
4 calls. Slowest were:
| 16.0ms | (* (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (/ 1 (* 2 a))) |
| 5.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 3.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
4 calls. Slowest were:
| 55.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 53.0ms | (* (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (/ 1 (* 2 a))) |
| 50.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 30.0ms | (- (* b b) (* 4 (* a c))) |
45 calls. Slowest were:
| 548.0ms | (* (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) 1) |
| 513.0ms | (* (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) 1) |
| 421.0ms | (* (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) (/ 1 (* 2 a))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 23.7b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 15.4b | (* (* (* 4 a) c) (/ 1/2 a)) |
| 1.0b | (- (* b b) (* 4 (* a c))) |
| 0.9b | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
4 calls. Slowest were:
| 6.0ms | (* (* (* 4 a) c) (/ 1/2 a)) |
| 5.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 2.0ms | (- (* b b) (* 4 (* a c))) |
4 calls. Slowest were:
| 57.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 52.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 43.0ms | (- (* b b) (* 4 (* a c))) |
| 11.0ms | (* (* (* 4 a) c) (/ 1/2 a)) |
29 calls. Slowest were:
| 339.0ms | (+ (log (* (* 4 a) c)) (log (/ 1/2 a))) |
| 314.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 284.0ms | (* (* (* 4 a) c) (/ 1/2 a)) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 23.7b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 1.0b | (- (* b b) (* 4 (* a c))) |
| 0.9b | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.1b | (/ (* 2 c) (+ (- b) (sqrt (- (* b b) (* 4 (* a c)))))) |
4 calls. Slowest were:
| 20.0ms | (/ (* 2 c) (+ (- b) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 11.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 4.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
4 calls. Slowest were:
| 151.0ms | (/ (* 2 c) (+ (- b) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 57.0ms | (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 50.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 26.0ms | (- (* b b) (* 4 (* a c))) |
37 calls. Slowest were:
| 359.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 348.0ms | (+ (* (- b) (- b)) (- (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c))))) (* (- b) (sqrt (- (* b b) (* 4 (* a c))))))) |
| 146.0ms | (sqrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 6.6b
95.4% (1.1b remaining)
Error of 6.6b against oracle of 5.5b and baseline of 29.4b