


Bits error versus a



Bits error versus b



Bits error versus c
Results
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 clear-num11.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 "Quadratic roots, full range"
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 20.0s)Debug log
1 calls. Slowest were:
| 77.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:
| 25.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 20.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 3.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 92.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 73.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 58.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 41.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 392.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 267.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 170.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:
| 10.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 6.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 2.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 116.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 58.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 48.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 27.0ms | (- (* b b) (* (* c a) 4)) |
30 calls. Slowest were:
| 402.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 180.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (+ (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 166.0ms | (- b (* 2 (/ (* a c) b))) |
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) |
| 20.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 74.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 70.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 59.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 48.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
24 calls. Slowest were:
| 425.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 372.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 153.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:
| 24.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 15.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 143.0ms | (/ (* (* c -4) a) (* (* 2 a) (+ (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 61.0ms | (+ (sqrt (- (* b b) (* (* c a) 4))) b) |
| 51.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 37.0ms | (- (* b b) (* (* c a) 4)) |
29 calls. Slowest were:
| 830.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)))) |
| 823.0ms | (/ (* (* (* (* 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)))) |
| 357.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