



Bits error versus a




Bits error versus b




Bits error versus c
| Original | 33.2 |
|---|---|
| Target | 20.4 |
| Herbie | 9.6 |
if b < -4.9529633933266e-08 or -1.432147794020809e-104 < b < -5.51185453096989e-127Initial program 53.5
Simplified53.5
Taylor expanded around -inf 7.8
Simplified7.8
if -4.9529633933266e-08 < b < -1.432147794020809e-104Initial program 34.9
Simplified34.9
rmApplied flip--34.9
Applied associate-/l/38.7
Simplified21.5
if -5.51185453096989e-127 < b < 8.046684077965657e+99Initial program 11.0
Simplified11.0
Taylor expanded around -inf 11.0
Simplified11.0
rmApplied *-un-lft-identity11.0
Applied *-un-lft-identity11.0
Applied distribute-lft-out--11.0
Applied associate-/l*11.1
if 8.046684077965657e+99 < b Initial program 44.6
Simplified44.6
Taylor expanded around inf 3.5
Final simplification9.6
herbie shell --seed 2019004 +o rules:numerics
(FPCore (a b c)
:name "The quadratic formula (r2)"
: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: 33.8s)Debug log
1 calls. Slowest were:
| 108.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 33.9b
Found 4 expressions with local error:
| 22.0b | (sqrt (fma (* -4 a) c (* b b))) |
| 0.9b | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 0.2b | (fma (* -4 a) c (* b b)) |
| 0.1b | (/ (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (* 2 a)) |
4 calls. Slowest were:
| 16.0ms | (/ (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (* 2 a)) |
| 7.0ms | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 2.0ms | (sqrt (fma (* -4 a) c (* b b))) |
4 calls. Slowest were:
| 60.0ms | (/ (- (- b) (sqrt (fma (* -4 a) c (* b b)))) (* 2 a)) |
| 55.0ms | (- (- b) (sqrt (fma (* -4 a) c (* b b)))) |
| 51.0ms | (sqrt (fma (* -4 a) c (* b b))) |
| 38.0ms | (fma (* -4 a) c (* b b)) |
93 calls. Slowest were:
| 341.0ms | (- (* (- b) (- b)) (* (sqrt (fma (* -4 a) c (* b b))) (sqrt (fma (* -4 a) c (* b b))))) |
| 165.0ms | (- (pow (- b) 3) (pow (sqrt (fma (* -4 a) c (* b b))) 3)) |
| 98.0ms | (- (* 2 (/ (* a c) b)) (* 2 b)) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 5.1b
Found 4 expressions with local error:
| 22.0b | (sqrt (fma (* a -4) c (* b b))) |
| 0.9b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.2b | (fma (* a -4) c (* b b)) |
| 0.1b | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
4 calls. Slowest were:
| 12.0ms | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
| 5.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 1.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 85.0ms | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
| 59.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 54.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 34.0ms | (fma (* a -4) c (* b b)) |
93 calls. Slowest were:
| 343.0ms | (- (* (- b) (- b)) (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b))))) |
| 215.0ms | (- (pow (- b) 3) (pow (sqrt (fma (* a -4) c (* b b))) 3)) |
| 95.0ms | (- (* 2 (/ (* a c) b)) (* 2 b)) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 5.1b
Found 4 expressions with local error:
| 22.0b | (sqrt (fma (* a -4) c (* b b))) |
| 0.9b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.2b | (fma (* a -4) c (* b b)) |
| 0.2b | (/ 1 (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
4 calls. Slowest were:
| 9.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 8.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 2.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 86.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 72.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 61.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 33.0ms | (fma (* a -4) c (* b b)) |
86 calls. Slowest were:
| 413.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (* (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 74.0ms | (- (* 2 (/ (* a c) b)) (* 2 b)) |
| 70.0ms | (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 5.1b
Found 4 expressions with local error:
| 22.0b | (sqrt (fma (* a -4) c (* b b))) |
| 0.9b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.3b | (/ (/ 1 (* 2 a)) (/ 1 (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 0.2b | (fma (* a -4) c (* b b)) |
4 calls. Slowest were:
| 32.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 10.0ms | (/ (/ 1 (* 2 a)) (/ 1 (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 2.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 72.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 65.0ms | (/ (/ 1 (* 2 a)) (/ 1 (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 48.0ms | (fma (* a -4) c (* b b)) |
| 44.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
134 calls. Slowest were:
| 620.0ms | (/ 1 (/ 1 (- (* (- b) (- b)) (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b))))))) |
| 570.0ms | (/ 1 (/ 1 (- (* (- b) (- b)) (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b))))))) |
| 560.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (* (- b) (sqrt (fma (* a -4) c (* b b)))))) |
10 alts after pruning (8 fresh and 2 done)
Merged error: 5.1b
85.3% (4.1b remaining)
Error of 9.6b against oracle of 5.5b and baseline of 33.2b