



Bits error versus a




Bits error versus b




Bits error versus c
| Original | 33.2 |
|---|---|
| Target | 20.7 |
| Herbie | 6.7 |
if b < -8.276414657528413e+95Initial program 58.1
Simplified58.1
rmApplied div-inv58.1
Taylor expanded around -inf 2.4
Simplified2.4
if -8.276414657528413e+95 < b < -6.321237892731252e-282Initial program 32.3
Simplified32.2
rmApplied div-inv32.3
rmApplied flip--32.4
Applied associate-*l/32.4
Simplified15.1
Taylor expanded around -inf 8.3
rmApplied add-sqr-sqrt8.6
if -6.321237892731252e-282 < b < 3.7085352995569145e+106Initial program 9.6
Simplified9.6
rmApplied div-inv9.7
if 3.7085352995569145e+106 < b Initial program 46.2
Simplified46.1
Taylor expanded around inf 3.7
Final simplification6.7
herbie shell --seed 2019008 +o rules:numerics
(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: 30.4s)Debug log
1 calls. Slowest were:
| 109.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 33.9b
Found 4 expressions with local error:
| 23.7b | (sqrt (fma (* -4 a) c (* b b))) |
| 0.8b | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 0.1b | (/ (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (* 2 a)) |
| 0.0b | (fma (* -4 a) c (* b b)) |
4 calls. Slowest were:
| 12.0ms | (/ (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (* 2 a)) |
| 5.0ms | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 1.0ms | (sqrt (fma (* -4 a) c (* b b))) |
4 calls. Slowest were:
| 79.0ms | (sqrt (fma (* -4 a) c (* b b))) |
| 61.0ms | (/ (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (* 2 a)) |
| 57.0ms | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 44.0ms | (fma (* -4 a) c (* b b)) |
93 calls. Slowest were:
| 284.0ms | (- (* (- b) (- b)) (* (sqrt (fma (* -4 a) c (* b b))) (sqrt (fma (* -4 a) c (* b b))))) |
| 182.0ms | (- (pow (- b) 3) (pow (sqrt (fma (* -4 a) c (* b b))) 3)) |
| 124.0ms | (/ (* (* (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (- (- b) (sqrt (fma (* -4 a) c (* b b))))) (- (- b) (sqrt (fma (* -4 a) c (* b b))))) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 23.7b | (sqrt (fma (* -4 a) c (* b b))) |
| 0.8b | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 0.2b | (* (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (/ 1 (* 2 a))) |
| 0.0b | (fma (* -4 a) c (* b b)) |
4 calls. Slowest were:
| 13.0ms | (* (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (/ 1 (* 2 a))) |
| 5.0ms | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 1.0ms | (sqrt (fma (* -4 a) c (* b b))) |
4 calls. Slowest were:
| 64.0ms | (* (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (/ 1 (* 2 a))) |
| 47.0ms | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 46.0ms | (sqrt (fma (* -4 a) c (* b b))) |
| 36.0ms | (fma (* -4 a) c (* b b)) |
92 calls. Slowest were:
| 647.0ms | (* (- (* (- b) (- b)) (* (sqrt (fma (* -4 a) c (* b b))) (sqrt (fma (* -4 a) c (* b b))))) 1) |
| 590.0ms | (* (- (pow (- b) 3) (pow (sqrt (fma (* -4 a) c (* b b))) 3)) 1) |
| 369.0ms | (* (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (/ 1 (* 2 a))) |
10 alts after pruning (9 fresh and 1 done)
Merged error: 6.6b
Found 4 expressions with local error:
| 23.7b | (sqrt (fma (* -4 a) c (* b b))) |
| 15.4b | (* (- 0 (* (* c a) -4)) (/ 1/2 a)) |
| 0.9b | (+ (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 0.1b | (/ (* (- 0 (* (* c a) -4)) (/ 1/2 a)) (+ (- b) (sqrt (fma (* -4 a) c (* b b))))) |
4 calls. Slowest were:
| 20.0ms | (/ (* (- 0 (* (* c a) -4)) (/ 1/2 a)) (+ (- b) (sqrt (fma (* -4 a) c (* b b))))) |
| 8.0ms | (+ (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 7.0ms | (* (- 0 (* (* c a) -4)) (/ 1/2 a)) |
4 calls. Slowest were:
| 162.0ms | (/ (* (- 0 (* (* c a) -4)) (/ 1/2 a)) (+ (- b) (sqrt (fma (* -4 a) c (* b b))))) |
| 44.0ms | (+ (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 39.0ms | (sqrt (fma (* -4 a) c (* b b))) |
| 5.0ms | (* (- 0 (* (* c a) -4)) (/ 1/2 a)) |
61 calls. Slowest were:
| 678.0ms | (* (- 0 (* (* c a) -4)) (/ 1/2 a)) |
| 493.0ms | (* (- (pow 0 3) (pow (* (* c a) -4) 3)) (/ 1/2 a)) |
| 489.0ms | (* (- (pow 0 3) (pow (* (* c a) -4) 3)) (/ 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 (fma (* -4 a) c (* b b))) |
| 0.9b | (+ (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 0.1b | (/ (* 2 c) (+ (- b) (sqrt (fma (* -4 a) c (* b b))))) |
| 0.0b | (fma (* -4 a) c (* b b)) |
4 calls. Slowest were:
| 7.0ms | (/ (* 2 c) (+ (- b) (sqrt (fma (* -4 a) c (* b b))))) |
| 4.0ms | (+ (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 1.0ms | (sqrt (fma (* -4 a) c (* b b))) |
4 calls. Slowest were:
| 142.0ms | (/ (* 2 c) (+ (- b) (sqrt (fma (* -4 a) c (* b b))))) |
| 72.0ms | (+ (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 66.0ms | (sqrt (fma (* -4 a) c (* b b))) |
| 50.0ms | (fma (* -4 a) c (* b b)) |
36 calls. Slowest were:
| 530.0ms | (+ (* (- b) (- b)) (- (* (sqrt (fma (* -4 a) c (* b b))) (sqrt (fma (* -4 a) c (* b b)))) (* (- b) (sqrt (fma (* -4 a) c (* b b)))))) |
| 104.0ms | (/ (* (* (* 2 c) (* 2 c)) (* 2 c)) (* (* (+ (- b) (sqrt (fma (* -4 a) c (* b b)))) (+ (- b) (sqrt (fma (* -4 a) c (* b b))))) (+ (- b) (sqrt (fma (* -4 a) c (* b b)))))) |
| 73.0ms | (* (exp (- b)) (exp (sqrt (fma (* -4 a) c (* b b))))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 6.6b
95% (1.2b remaining)
Error of 6.7b against oracle of 5.5b and baseline of 29.7b