



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.1 |
|---|---|
| Target | 20.7 |
| Herbie | 10.3 |
if b < -2.8671071746705253e+127Initial program 51.2
Simplified51.2
Taylor expanded around -inf 2.8
if -2.8671071746705253e+127 < b < 1.354454736393541e-170Initial program 11.4
Simplified11.3
rmApplied *-un-lft-identity11.3
Applied associate-/l*11.5
if 1.354454736393541e-170 < b < 9.22274670563189e-75Initial program 24.9
Simplified24.9
rmApplied flip--25.0
Applied associate-/l/31.2
Simplified25.9
rmApplied times-frac18.2
Simplified18.2
if 9.22274670563189e-75 < b < 1.849572225033813e-22Initial program 37.0
Simplified37.0
rmApplied add-sqr-sqrt37.5
Applied times-frac37.5
if 1.849572225033813e-22 < b Initial program 54.3
Simplified54.3
Taylor expanded around inf 6.8
Simplified6.8
Final simplification10.3
herbie shell --seed 2019005
(FPCore (a b c)
:name "The quadratic formula (r1)"
:herbie-target
(if (< b 0) (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) (/ c (* a (/ (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 18.8s)Debug log
1 calls. Slowest were:
| 114.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 32.2b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* c a) 4))) |
| 1.1b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.4b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
4 calls. Slowest were:
| 20.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 19.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 69.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 67.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 64.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 24.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 377.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 366.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 251.0ms | (- (pow (sqrt (- (* b b) (* (* c a) 4))) 3) (pow b 3)) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 6.5b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* c a) 4))) |
| 1.1b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.4b | (- (* b b) (* (* c a) 4)) |
| 0.2b | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
4 calls. Slowest were:
| 20.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 9.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 99.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 55.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 51.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 29.0ms | (- (* b b) (* (* c a) 4)) |
30 calls. Slowest were:
| 394.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 228.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (+ (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 142.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 6.4b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* c a) 4))) |
| 23.2b | (sqrt (- (* b b) (* (* c a) 4))) |
| 1.1b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 1.1b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
4 calls. Slowest were:
| 21.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 16.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 76.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 73.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 54.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 52.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
24 calls. Slowest were:
| 404.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 382.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 147.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 6.4b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* c a) 4))) |
| 16.3b | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 1.1b | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 0.4b | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 12.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 11.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 137.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 70.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 68.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 43.0ms | (- (* b b) (* (* c a) 4)) |
29 calls. Slowest were:
| 1.1s | (/ (* (* (* (* c -4) a) (* (* c -4) a)) (* (* c -4) a)) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (+ (sqrt (- (* b b) (* (* c a) 4))) b) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) (+ (sqrt (- (* b b) (* (* c a) 4))) b)))) |
| 783.0ms | (/ (* (* (* (* c -4) a) (* (* c -4) a)) (* (* c -4) a)) (* (* (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b)))) |
| 323.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 6.4b
84.2% (4.3b remaining)
Error of 10.3b against oracle of 6.0b and baseline of 33.2b